Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,161,067 C→T 6.3% R162H (CGT→CAT)  fhuE ← ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,161,0670CT 6.3% 204.1 / 2.8 63R162H (CGT→CAT) fhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

TGTAAATCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGATC  >  NC_000913/1160920‑1161206
                                                                                                                                                   |                                                                                                                                           
tGTAAATCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACTACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCAc                                                                                                                                             >  5:36400/1‑148 (MQ=255)
tGTAAATCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCAc                                                                                                                                             >  3:24368/1‑148 (MQ=255)
  tAATTCCGCCACATACCGTTCTTTGTTCCAGATACCTTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATTTATTACCCTTCCCGGTAATGAGTACTGTCGCGCCACGc                                                                                                                                           <  6:233766/148‑1 (MQ=255)
  tAAATCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGc                                                                                                                                           <  4:24368/148‑1 (MQ=255)
  tAAATCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGc                                                                                                                                           <  8:226144/148‑1 (MQ=255)
     aTCCGCCACATACCGTGCTTTGTTCCAGCTACCGTATTCCGCCGCGACATCGCCTTTAAATTCACGACTGGTCGCGGGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCact                                                                                                                                        <  5:76357/148‑1 (MQ=255)
      tCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCacta                                                                                                                                       <  5:310388/148‑1 (MQ=255)
           cacaTACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCt                                                                                                                                  <  6:36400/148‑1 (MQ=255)
           cacaTACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCt                                                                                                                                  <  8:132976/148‑1 (MQ=255)
            acaTACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTa                                                                                                                                 <  5:317055/148‑1 (MQ=255)
               tACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTAcgc                                                                                                                              >  7:246063/1‑148 (MQ=255)
               tACCGTTCTTTGTTCCAGCTACCATATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTAcgc                                                                                                                              <  7:326116/148‑1 (MQ=255)
                aCCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTAcgcg                                                                                                                             >  1:330597/1‑148 (MQ=255)
                      cTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCaa                                                                                                                        <  8:246063/147‑1 (MQ=255)
                      cTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCaa                                                                                                                        >  8:141789/1‑147 (MQ=255)
                       tttGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAag                                                                                                                      <  2:193725/148‑1 (MQ=255)
                         tGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAag                                                                                                                      <  2:330597/146‑1 (MQ=255)
                         tGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAag                                                                                                                      <  2:344782/146‑1 (MQ=255)
                              cAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGcc                                                                                                                <  1:116536/147‑1 (MQ=255)
                               aGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCa                                                                                                               >  3:162140/1‑147 (MQ=255)
                               aGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCa                                                                                                               >  3:156825/1‑147 (MQ=255)
                                      gTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGa                                                                                                        >  3:141031/1‑147 (MQ=255)
                                        aTTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCAc                                                                                                                                             >  1:227215/1‑108 (MQ=255)
                                        aTTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCAc                                                                                                                                             <  2:227215/108‑1 (MQ=255)
                                        aTTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGc                                                                                                                                           <  8:253920/110‑1 (MQ=255)
                                        aTTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGc                                                                                                                                           >  7:253920/1‑110 (MQ=255)
                                                       cGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCATGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGAtt                                                                                      >  2:241386/1‑148 (MQ=255)
                                                       cGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCATGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGAtt                                                                                      >  2:135437/1‑148 (MQ=255)
                                                        gCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTc                                                                                     >  1:302413/1‑148 (MQ=255)
                                                             tAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGc                                                                                                                                           <  2:362462/89‑1 (MQ=255)
                                                             tAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGc                                                                                                                                           >  1:362462/1‑89 (MQ=255)
                                                                aTTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCaa                                                                                                                        >  7:377749/1‑105 (MQ=255)
                                                                aTTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCaa                                                                                                                        <  8:377749/105‑1 (MQ=255)
                                                                  tCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCactact                                                                                                                                     >  2:85464/1‑90 (MQ=255)
                                                                  tCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCactact                                                                                                                                     <  1:85464/90‑1 (MQ=255)
                                                                    aCGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACGCGTCCCGGTCATGAGTCCTGTCGCGCCTCGAAATACTTCTACGCGTTCAAAGAGTGCCTTATCATAAAGTGCGTCGCCCTGATTCCAGCGCGATTCa                                                                         <  2:375463/148‑1 (MQ=255)
                                                                        cTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAAt                                                                     <  2:302413/148‑1 (MQ=255)
                                                                           gTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCTCCCAGATTCCAGCTCGATTCAAAATAg                                                                   <  4:156825/147‑1 (MQ=255)
                                                                           gTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAg                                                                   <  4:162140/147‑1 (MQ=255)
                                                                           gTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAg                                                                   <  4:141031/147‑1 (MQ=255)
                                                                                               atTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAAt                                                                     <  5:343129/125‑1 (MQ=255)
                                                                                               atTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAAt                                                                     >  6:343129/1‑125 (MQ=255)
                                                                                                ttAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCata                                              >  2:390/1‑147 (MQ=255)
                                                                                                   aTTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAg                                           <  7:141789/147‑1 (MQ=255)
                                                                                                           aGATGGATTACCAGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATATGTGGGGATACCATCAACCATAAAGTTATCGCTc                                  >  6:197/1‑148 (MQ=255)
                                                                                                             aTGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCat                                               >  1:101210/1‑133 (MQ=255)
                                                                                                             aTGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCat                                               <  2:101210/133‑1 (MQ=255)
                                                                                                               ggATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATGCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTgg                               >  2:129737/1‑147 (MQ=255)
                                                                                                                   tACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATcc                          >  3:36574/1‑148 (MQ=255)
                                                                                                                               tCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCaa                                                                        <  8:33259/90‑1 (MQ=255)
                                                                                                                               tCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCaa                                                                        >  7:33259/1‑90 (MQ=255)
                                                                                                                               tCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGa                              <  1:393468/132‑1 (MQ=255)
                                                                                                                               tCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGa                              >  2:393468/1‑132 (MQ=255)
                                                                                                                                cATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAataata              >  3:304043/1‑147 (MQ=255)
                                                                                                                                 aTGAGTCCTGTCGCGCCATGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAAtaataa             <  1:135437/147‑1 (MQ=255)
                                                                                                                                 aTGAGTCCTGTCGCGCCATGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAAtaataa             <  1:241386/147‑1 (MQ=255)
                                                                                                                                  tGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAataata              >  5:13906/1‑145 (MQ=255)
                                                                                                                                   gAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAATAATaaa            >  4:89630/1‑146 (MQ=255)
                                                                                                                                    aGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGaa                             <  8:7512/128‑1 (MQ=255)
                                                                                                                                    aGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGaa                             <  8:13356/128‑1 (MQ=255)
                                                                                                                                    aGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGaa                             >  7:7512/1‑128 (MQ=255)
                                                                                                                                    aGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGaa                             >  7:13356/1‑128 (MQ=255)
                                                                                                                                         tGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGa    >  6:364176/1‑148 (MQ=255)
                                                                                                                                          gTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGa                                    <  7:264069/115‑1 (MQ=255)
                                                                                                                                          gTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGa                                    >  8:264069/1‑115 (MQ=255)
                                                                                                                                          gTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCATCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCACCAACCATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGAt   >  2:62650/1‑148 (MQ=255)
                                                                                                                                           tCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGATc  >  4:88866/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TGTAAATCCGCCACATACCGTTCTTTGTTCCAGCTACCGTATTCCGCCGAGACATCGCCTTTAAATTCACGACTGGTCGCGTGTTTTCGAACCATATTAATTGCCGCAGATGGATTACCCGTCCCGGTCATGAGTCCTGTCGCGCCACGCACTACTTCTACGCGTTCAAAGAGTGCCATATCAGAAAGTGCGTCGCCCAGATTCCAGCGCGATTCAAAATAGGTGGGGATACCATCAACCATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGATC  >  NC_000913/1160920‑1161206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: