Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 443,023 T→C 5.2% intergenic (+26/+28) thiI → / ← yajL tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913443,0230TC 5.2% 682.7 / 17.2 209intergenic (+26/+28)thiI/yajLtRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand):  ref base T (95/103);  new base C (6/5);  total (101/108)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-03

Only 100 of 210 total aligned reads displayed.

TTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGG  >  NC_000913/442875‑443170
                                                                                                                                                    |                                                                                                                                                   
    gCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGa                                                                                                                                                 <  2:1141549/149‑1 (MQ=255)
     cGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                <  3:769574/149‑1 (MQ=255)
     cGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                <  2:1192968/149‑1 (MQ=255)
         ctcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGa                                                                                                                                            <  3:1438408/149‑1 (MQ=255)
            gACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAc                                                                                                                                           >  1:1161495/1‑147 (MQ=255)
              ccAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGttt                                                                                                                                       >  1:1699338/1‑149 (MQ=255)
               cAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTg                                                                                                                                      >  3:701335/1‑149 (MQ=255)
                 gAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                                                                    <  1:1597123/149‑1 (MQ=255)
                  aaCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                >  4:966915/1‑136 (MQ=255)
                  aaCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                <  3:966915/136‑1 (MQ=255)
                  aaCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGc                                                                                                                                              >  1:35080/1‑138 (MQ=255)
                  aaCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGc                                                                                                                                              <  2:35080/138‑1 (MQ=255)
                   aCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGc                                                                                                                                  <  3:1081334/149‑1 (MQ=255)
                       aaCCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATcc                                                                                                                              >  1:793824/1‑149 (MQ=255)
                        aCCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCg                                                                                                                             >  3:624822/1‑149 (MQ=255)
                         ccTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGa                                                                                                                            <  4:880897/149‑1 (MQ=255)
                           tGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGaca                                                                                                                          <  1:1239309/149‑1 (MQ=255)
                                gctgTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                                                                     <  2:317419/133‑1 (MQ=255)
                                  tgtgGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCg                                                                                                                   >  2:1305136/1‑149 (MQ=255)
                                       tgtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGGAGGTATATCGCCCGTAATTCGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCTGCATCCGACACTAGTCGtaata                                                                                                              >  1:967008/1‑149 (MQ=255)
                                             cgcgGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtata                                                                                                        >  2:1129743/1‑149 (MQ=255)
                                                 gggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATc                                                                                                    <  1:55105/149‑1 (MQ=255)
                                                    tgatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCt                                                                                                 <  1:958668/149‑1 (MQ=255)
                                                     gatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGa                                                                                                                                            <  4:485971/105‑1 (MQ=255)
                                                     gatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGa                                                                                                                                            >  3:485971/1‑105 (MQ=255)
                                                      atgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgc                                                                                               <  2:610432/149‑1 (MQ=255)
                                                        gaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtat                                                                                                         <  3:35165/137‑1 (MQ=255)
                                                        gaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccg                                                                                             >  2:1313989/1‑149 (MQ=255)
                                                         aGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccgc                                                                                            >  2:281898/1‑149 (MQ=255)
                                                              tCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                >  2:1217004/1‑92 (MQ=255)
                                                                tGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAcc                                                                                     >  1:582830/1‑149 (MQ=255)
                                                                 gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCa                                                                                    <  1:734852/149‑1 (MQ=255)
                                                                   aGGCGCTCTATCTGAGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATCCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCGAACACCAGt                                                                                  <  2:1256388/149‑1 (MQ=255)
                                                                       gcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgtg                                                                              >  3:1306220/1‑149 (MQ=255)
                                                                       gcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgtg                                                                              >  2:835118/1‑149 (MQ=255)
                                                                         gctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTCGTTAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGAt                                                                            >  4:1462646/2‑149 (MQ=255)
                                                                           tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGc                                                                          <  2:803726/149‑1 (MQ=255)
                                                                           tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGc                                                                          <  3:1628582/149‑1 (MQ=255)
                                                                             tCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCa                                                                        >  1:50469/1‑149 (MQ=255)
                                                                                gcgcGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACtt                                                                     <  4:225882/149‑1 (MQ=255)
                                                                                 cgcgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTc                                                                    >  4:1546928/1‑149 (MQ=255)
                                                                                   cgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                  <  4:28730/149‑1 (MQ=255)
                                                                                   cgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                  <  1:281898/149‑1 (MQ=255)
                                                                                   cgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                  <  4:200365/149‑1 (MQ=255)
                                                                                         gggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtt                                                            >  3:1680652/1‑149 (MQ=255)
                                                                                          ggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttt                                                           <  2:1428778/149‑1 (MQ=255)
                                                                                          ggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttt                                                           <  1:814186/149‑1 (MQ=255)
                                                                                           gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtttt                                                          <  4:701335/149‑1 (MQ=255)
                                                                                             tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                        >  3:246995/1‑149 (MQ=255)
                                                                                             tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                        >  1:1538494/1‑149 (MQ=255)
                                                                                             tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                        >  1:151432/1‑149 (MQ=255)
                                                                                              ttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCa                                                       <  2:190729/149‑1 (MQ=255)
                                                                                               tAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCAc                                                      <  2:436280/149‑1 (MQ=255)
                                                                                                 aCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                       >  3:420016/1‑114 (MQ=255)
                                                                                                 aCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc                                                    >  3:306770/1‑149 (MQ=255)
                                                                                                   aacgagAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGccc                                                  >  2:74665/6‑149 (MQ=255)
                                                                                                   aaTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                                                                            <  2:1080768/91‑1 (MQ=255)
                                                                                                     tgtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCa                                                                                    <  4:1333583/113‑1 (MQ=255)
                                                                                                          aGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCcaaccaac                                           <  3:984026/149‑1 (MQ=255)
                                                                                                           ggTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACa                                          <  3:352821/149‑1 (MQ=255)
                                                                                                            gTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAg                                         <  4:171971/149‑1 (MQ=255)
                                                                                                            gTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAg                                         <  1:1058728/149‑1 (MQ=255)
                                                                                                               tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTc                                      <  4:624822/149‑1 (MQ=255)
                                                                                                               tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTc                                      <  2:214026/149‑1 (MQ=255)
                                                                                                                 tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGa                                    <  4:591446/149‑1 (MQ=255)
                                                                                                                 tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGa                                    <  4:361898/149‑1 (MQ=255)
                                                                                                                  cGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                        >  2:98721/1‑128 (MQ=255)
                                                                                                                  cGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                        <  1:98721/128‑1 (MQ=255)
                                                                                                                  cGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGAt                                   >  4:1141836/1‑149 (MQ=255)
                                                                                                                     ccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATTGGCTTTTTCACGCCCAACCAACAGGTCGATAAt                                >  3:372050/1‑149 (MQ=255)
                                                                                                                     ccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAt                                <  1:466287/149‑1 (MQ=255)
                                                                                                                      cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAtt                               >  1:259940/1‑149 (MQ=255)
                                                                                                                       gTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAttt                              >  3:1489184/1‑149 (MQ=255)
                                                                                                                            ttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGa                         <  1:1313989/149‑1 (MQ=255)
                                                                                                                              gtgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGAcc                       >  3:1658909/1‑149 (MQ=255)
                                                                                                                                ggTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAAATATAAATCCCTGCCGCCATCACCAGTTGTGATACCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaa                     >  4:1607748/1‑149 (MQ=255)
                                                                                                                                 gTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaaa                    >  1:1463159/1‑149 (MQ=255)
                                                                                                                                  tttttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATg                                                                 <  1:1435127/103‑1 (MQ=255)
                                                                                                                                   ttttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGt                  <  4:1658909/149‑1 (MQ=255)
                                                                                                                                       aCGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGAt              <  4:352168/149‑1 (MQ=255)
                                                                                                                                         gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaaccaa                                            >  2:963977/1‑117 (MQ=255)
                                                                                                                                         gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg            >  1:334039/1‑149 (MQ=255)
                                                                                                                                         gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg            >  2:1505322/1‑149 (MQ=255)
                                                                                                                                         gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg            <  3:1141836/149‑1 (MQ=255)
                                                                                                                                          tCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                       <  2:1519553/73‑1 (MQ=255)
                                                                                                                                          tCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc           >  1:712414/1‑149 (MQ=255)
                                                                                                                                          tCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTGGTGATGCCACTTCATGGGCtttt                                                          >  2:1590860/1‑102 (MQ=255)
                                                                                                                                          tCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAGTAATTATAAATCCCTGCCGCCATCACCAGTTGTGAt                                                                            >  4:1376031/1‑84 (MQ=39)
                                                                                                                                          tCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGAt                                                                            <  3:1376031/84‑1 (MQ=255)
                                                                                                                                          tCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc           >  1:406671/1‑149 (MQ=255)
                                                                                                                                           cGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCt          <  2:658919/149‑1 (MQ=255)
                                                                                                                                           cGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCt          <  4:443160/149‑1 (MQ=255)
                                                                                                                                            gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                             <  4:851920/81‑1 (MQ=255)
                                                                                                                                            gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                             >  3:851920/1‑81 (MQ=255)
                                                                                                                                            gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTg         >  4:1615029/1‑149 (MQ=255)
                                                                                                                                              aTCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                  >  2:120393/1‑90 (MQ=255)
                                                                                                                                              aTCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                  <  1:120393/90‑1 (MQ=255)
                                                                                                                                               tCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaaa                    >  3:924349/1‑135 (MQ=255)
                                                                                                                                               tCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaaa                    <  4:924349/135‑1 (MQ=255)
                                                                                                                                                   gTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg  >  2:372856/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGG  >  NC_000913/442875‑443170

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