Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,357 C→T 5.6% intergenic (‑28/+43) tufA ← / ← fusA translation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,471,3570CT 5.6% 388.2 / 7.4 125intergenic (‑28/+43)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (53/65);  new base T (3/4);  total (56/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

Only 100 of 130 total aligned reads displayed.

CGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCA  >  NC_000913/3471209‑3471506
                                                                                                                                                     |                                                                                                                                                     
cgTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACAGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTccc                                                                                                                                                      <  1:269376/150‑1 (MQ=255)
 gTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcccc                                                                                                                                                     >  2:63349/1‑150 (MQ=255)
 gTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcccc                                                                                                                                                     >  2:843753/1‑150 (MQ=255)
    ggTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                                                                  >  1:1107516/1‑150 (MQ=255)
    ggTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATACCTTACTATAGTGCTCTCCCCTTc                                                                                                                                                  >  1:929312/1‑150 (MQ=255)
       tttAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGAGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAgg                                                                                                                                               >  1:988618/1‑150 (MQ=255)
           gCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCATTACTATAGTGCTCTCCCCTTCAGgagag                                                                                                                                           >  1:1242729/1‑150 (MQ=255)
           gCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTCTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                                                                           >  1:740529/1‑150 (MQ=255)
            ccAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                                                                          >  1:1383493/1‑150 (MQ=255)
             cAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                                                                         >  2:102900/1‑150 (MQ=255)
             cAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTACagaagacgc                                                                                                                                         >  2:1580033/1‑147 (MQ=255)
              aGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGc                                                                                                                                        <  1:779217/150‑1 (MQ=255)
                tACGGGGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                                                                      <  1:1350478/150‑1 (MQ=21)
                 ccGGAGGTGCTTGCAGCGGTCAGAGGAGTTTTACCGTGGTCAACGTGGCCGAACGTACCAACGCTAACGTGCGGTTTGGTACGTTCAACTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCACTCCTGAAGGGGAGAGCAct                                                                                                                                     <  1:42876/149‑2 (MQ=255)
                 aCGGGGGTGATTGCAGCGGGCAGAGTAGGTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACg                                                                                                                                     <  2:1525676/150‑1 (MQ=255)
                  cGGTGGTGATTGCAGCGGTCAGTGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCATCGTTAACGTGCGGTTTTGTTCGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                                                                    <  2:367674/150‑1 (MQ=255)
                  cGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                                                                    >  1:554313/1‑150 (MQ=255)
                   ggtggtGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                                                                   <  2:1630627/150‑1 (MQ=255)
                   ggtggtGATTGCAGCGGTCAGAGTAGATTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAAATTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                                                                   <  1:1227232/150‑1 (MQ=255)
                        tgATTGCAGCGGTCAGAGGAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTATCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                                                                              <  1:1010278/150‑1 (MQ=255)
                        ggcttGCAGCCGTCAGAGTCGTTTTACCGTGGTCCACGTGGCCGATCGTACCCACGTTAACCTGCGGGTTTGTACGTTCAAATTTTTCTTTAGACCCGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                                                                              <  2:19838/147‑1 (MQ=255)
                         tttgtcAGCGGTAAGAGTCGTTTTCCGGGGGTAAACGTGGCCGATAGTACCAAGTTTATCGTGGGGTTTTTTACTTTAAAATTTCTCTTTAGACACGGCTATATTCCTTAATATTGTGCTCTCCCCTTCTGGAGAGAGCACGGGACTTTg                                                                                                                             <  2:740529/145‑1 (MQ=255)
                          aTTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                                                                            >  2:1361614/1‑150 (MQ=255)
                              cAGCGGTCAGAGATGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGCGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACAATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAtt                                                                                                                        <  1:193397/150‑1 (MQ=255)
                               aGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTa                                                                                                                       >  1:623370/1‑150 (MQ=255)
                                  gggcAGAGTAGTTTTACCGGGGACAACGTGGCCGATAGAACCAACGTTAACGTGCGGTTTTGTACGTACAAATTTTTCTTTAGACACGGCTATATTCCTAACTATAGTGCTCCCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                    <  1:371916/147‑1 (MQ=37)
                                     cAGAGTAGTTTTACCGTGGCCAACGTGGCCGATAGTACAAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTAAGACACGGCTATATTCCTTACTATAGTGCGCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCtt                                                                                                                 <  1:733542/150‑1 (MQ=255)
                                      agagTAGGTTTACCGTGGGCAACTTGGCCGATAGTACCAACGTAAACGTGCGGTTTTGTACGTTTAAATTATTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTACTAACCCTTa                                                                                                                <  1:902610/150‑1 (MQ=255)
                                       gagTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAg                                                                                                               >  2:388075/1‑150 (MQ=255)
                                           agtTTTACCGTGGTCAACTTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTAAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCtt                                                                                                           <  2:1663988/150‑1 (MQ=255)
                                            gtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTa                                                                                                          >  2:762560/1‑150 (MQ=255)
                                              tttACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAtt                                                                                                        <  2:1127424/150‑1 (MQ=255)
                                               ttACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAttt                                                                                                       <  2:263925/150‑1 (MQ=255)
                                                 aCCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                                                                     >  1:156143/1‑150 (MQ=255)
                                                 aCCGTGGGCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTAGTTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCGCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                                                                     <  1:495817/150‑1 (MQ=255)
                                                  ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAcc                                                                                                    >  1:275525/1‑150 (MQ=255)
                                                  ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCGTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGTTATTAACCCTTAGGCTTATTTAcc                                                                                                    >  2:119825/1‑150 (MQ=255)
                                                       gTCAACGTGTCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                                                                               <  1:54530/150‑1 (MQ=255)
                                                       gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                                                                               >  1:813428/1‑150 (MQ=255)
                                                        tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                              >  2:1448900/1‑150 (MQ=255)
                                                        tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                              <  1:899379/150‑1 (MQ=255)
                                                        tCAACATGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                              <  2:914934/150‑1 (MQ=255)
                                                         cAACGTGGCCGATAGTACCAACGTTAACGTGCGGTATTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCATACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGCCTTATTTACCACGGGCt                                                                                             >  1:18524/1‑150 (MQ=255)
                                                          atcGTGGCCTAAACTACCACCTTTAACGTGCCGTTTTGAACGTCCAAATTTTTCTTAACACAGGCCAAAATACCTTACTATTGTCCTCACCCCTTCCGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCtt                                                                                            <  2:165318/148‑1 (MQ=255)
                                                            cGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCa                                                                                          <  2:1204019/150‑1 (MQ=255)
                                                             gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCaa                                                                                         >  1:436997/1‑150 (MQ=255)
                                                              tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCaa                                                                                         >  2:659408/1‑150 (MQ=255)
                                                                   gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                                                                   <  1:297970/150‑1 (MQ=255)
                                                                   gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                                                                   <  1:195439/150‑1 (MQ=255)
                                                                    aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                                                                  >  2:1624542/1‑150 (MQ=255)
                                                                     tAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACTCGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGcc                                                                                 <  2:929868/150‑1 (MQ=255)
                                                                     tAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGcc                                                                                 <  2:504131/150‑1 (MQ=255)
                                                                     tAGTAAAAACGTTAACGTGCGGTTTTGTACGTACAAATTTTTATTTAGACACGGCTTTATTACTTACTATATTGCTCTCCCCTTCATGATAGATCACGTGACTTTGGAATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTActtcc                                                                                 >  2:592802/1‑146 (MQ=255)
                                                                        tACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACAACGGGCTTCAATTACGGCCTGa                                                                              >  1:1647304/1‑150 (MQ=255)
                                                                        tACCAACGTTAACGTGCGGATTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTATGCTTATTTACCACGGGCTTCAATTACGGCCTGa                                                                              >  2:1218310/1‑150 (MQ=255)
                                                                         aCCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCGATATTCCTTACAATAGAGCTCTCCCCATAATGAGAGAGCACTGGACTTTGGTATTAACCCATAGGCTTATTTACCAAGGGATTCAATTACGGCCTgcg                                                                             >  2:1047281/1‑148 (MQ=37)
                                                                          ccAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                                                                            <  1:1626745/150‑1 (MQ=255)
                                                                              cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                                                                        >  1:1223228/1‑150 (MQ=255)
                                                                               tttAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                                                                       <  1:530837/149‑1 (MQ=255)
                                                                                  aaCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTAc                                                                    <  2:1616171/150‑1 (MQ=255)
                                                                                  aaCGGGTGGTTTTGAACGTTAAAATTTTTCTTTACACCCGGCAATATTCCTTACTATTGTGCTCTCCCCTTCTGTATATAGAACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTAc                                                                    <  1:363115/150‑1 (MQ=255)
                                                                                        cGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg                                                              >  2:1021735/1‑150 (MQ=255)
                                                                                          gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt                                                            >  2:1233529/1‑150 (MQ=255)
                                                                                          gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt                                                            >  2:186918/1‑150 (MQ=255)
                                                                                            tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc                                                          >  1:1275413/1‑150 (MQ=255)
                                                                                             ttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca                                                         <  1:1000737/150‑1 (MQ=255)
                                                                                              tttaccgtCAAATTTTTTATTAGCCACGGCTAATTTCCTTAAAAAAGGGCCCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat                                                        <  1:1580033/143‑1 (MQ=255)
                                                                                               gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGTGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc                                                       <  1:843753/150‑1 (MQ=255)
                                                                                               gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCTCGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc                                                       <  2:1242729/150‑1 (MQ=255)
                                                                                               gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc                                                       <  1:63349/150‑1 (MQ=255)
                                                                                                   gTTCAAATTTTGCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATAc                                                   <  2:481635/150‑1 (MQ=255)
                                                                                                        aaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCACTTCTGGAGAGAGCACGGGACTTTGGTATTAACCATTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAg                                              >  1:115408/1‑150 (MQ=255)
                                                                                                         aTTTTCCTTTCGACACGCCTGTATACCTTCCTCTCGGGGTCACCCCTTCTTGTCAGAGCACGGGACGTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAgg                                             <  2:185463/150‑1 (MQ=255)
                                                                                                             ttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAAtt                                         >  1:1009502/1‑150 (MQ=255)
                                                                                                                tttAGACACGGCTATATTCCTTACTATAGCGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCa                                      <  1:1690069/150‑1 (MQ=255)
                                                                                                                 ttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCAt                                     <  1:431155/150‑1 (MQ=255)
                                                                                                                 ttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCAt                                     >  1:530742/1‑150 (MQ=255)
                                                                                                                   aGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAg                                   >  2:145107/1‑150 (MQ=255)
                                                                                                                     acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAgtg                                 <  1:102900/150‑1 (MQ=255)
                                                                                                                      cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAgtgt                                >  2:902126/1‑150 (MQ=255)
                                                                                                                       acGGCTATATTCCTTACTATAGTGCGCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGTAATTCCATAGTGTa                               <  2:851381/150‑1 (MQ=255)
                                                                                                                          gCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTatga                            <  2:1107516/150‑1 (MQ=255)
                                                                                                                          gCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTatga                            >  1:1565996/1‑150 (MQ=255)
                                                                                                                          gCTATATTCCTTACAATAGTGCTCTCCCCTGCAGGAGAGAGCCCGGGACTGTGGTATTAACCCTTAGGCTTATTTACCACGGGATTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCCTAGTGTatga                            >  2:1478482/1‑150 (MQ=255)
                                                                                                                            tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTatgatg                          <  2:554313/150‑1 (MQ=255)
                                                                                                                                ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACg                      <  2:275525/150‑1 (MQ=255)
                                                                                                                                  ccTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGAc                    >  2:1694200/1‑150 (MQ=255)
                                                                                                                                  ccTTACTATAGTGCTCTCCCCTTCAGGAGAGAGAACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAAATACGGCCTGAGCCACGTTACTCGGCGCTTCATCATACGTCAGGCATTACATAGTGTAAGATTCACGAc                    >  1:315308/1‑150 (MQ=255)
                                                                                                                                    ttACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCt                  <  1:1336019/150‑1 (MQ=255)
                                                                                                                                        tataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTgg              >  1:718227/1‑150 (MQ=255)
                                                                                                                                          taGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGGTGCACGACCTTTGGTc            <  2:1147014/150‑1 (MQ=255)
                                                                                                                                          taGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTc            <  2:618231/150‑1 (MQ=255)
                                                                                                                                          ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTc            >  1:1109267/3‑150 (MQ=255)
                                                                                                                                            ttGCTCACCCCTTCAGGAGAGAGCACGGGACTTTGGTAATAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCag          <  2:854963/149‑1 (MQ=255)
                                                                                                                                            gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCag          <  2:1580514/150‑1 (MQ=255)
                                                                                                                                              gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCagag        <  2:662426/150‑1 (MQ=255)
                                                                                                                                              gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCagag        >  1:857980/1‑150 (MQ=255)
                                                                                                                                                  tcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACg    <  2:1613717/150‑1 (MQ=255)
                                                                                                                                                    ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGc   >  1:983584/1‑150 (MQ=255)
                                                                                                                                                     cccTTCAGGAGAGAGCACGGGACTTTGGTATTACCCCTTAGGCTAATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTCCTCGGCGCTTCATCATACTTCAGGAATTACATAGTGTATGAAGCACGACCTGTGGTCAGAGAACGCa  >  1:209743/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCA  >  NC_000913/3471209‑3471506

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