Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 769,791 A→T 5.2% D605V (GAT→GTT)  mngB → alpha‑mannosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913769,7910AT 5.2% 249.0 / 2.9 76D605V (GAT→GTT) mngBalpha‑mannosidase
Reads supporting (aligned to +/- strand):  ref base A (30/42);  new base T (2/2);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAA  >  NC_000913/769643‑769939
                                                                                                                                                    |                                                                                                                                                    
cATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGat                                                                                                                                                     >  1:379928/1‑150 (MQ=255)
 aTGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGata                                                                                                                                                    <  2:648419/150‑1 (MQ=255)
  tGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGatat                                                                                                                                                   >  2:1316202/1‑150 (MQ=255)
  tGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGataa                                                                                                                                                   >  2:1423571/1‑149 (MQ=255)
     aGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCaa                                                                                                                                                >  1:972953/1‑150 (MQ=255)
        cTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCt                                                                                                                                             <  2:628463/150‑1 (MQ=255)
             agagCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATa                                                                                                                                        <  1:769793/150‑1 (MQ=255)
                  aGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGTTATCAAGATTGATATCCac                                                                                                                                   <  2:31245/150‑2 (MQ=255)
                  aGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAg                                                                                                                                   >  1:32327/1‑150 (MQ=255)
                   gggCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAAt                                                                                                                                               >  2:1203250/1‑137 (MQ=255)
                   gggCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAAt                                                                                                                                               <  1:1203250/137‑1 (MQ=255)
                   gggCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAgg                                                                                                                                  >  1:607816/1‑150 (MQ=255)
                    ggCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGc                                                                                                                                 <  2:1219981/150‑1 (MQ=255)
                     gCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCt                                                                                                                                >  1:664152/1‑150 (MQ=255)
                     gCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCt                                                                                                                                >  1:331860/1‑150 (MQ=255)
                        gTTATCCGCTATCACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACGGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGAGATCAATCTTGATAACCAGGCGGAc                                                                                                                             >  2:499115/1‑150 (MQ=255)
                           aTCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGat                                                                                                                          >  1:1039770/1‑150 (MQ=255)
                             ccGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATACCAATCTTGATAACCAGGCTGACGatca                                                                                                                        <  1:1423571/150‑1 (MQ=255)
                                 aaTGACATGGCAGTGCCGCTCATTTTGTATGACTGCAGCGCCCGGCAATCCCCTGTCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCCCATCAGCATGCGTATTGATGTGGATAGCAATCTTGATAACCAGGCTGACGATCATCGc                                                                                                                    <  2:1002886/149‑1 (MQ=255)
                                     aCATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTc                                                                                                                >  1:666041/1‑150 (MQ=255)
                                         ggCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgtgt                                                                                                            <  1:379037/150‑1 (MQ=255)
                                         ggCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgtgt                                                                                                            <  2:607816/150‑1 (MQ=255)
                                          gCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTc                                                                                                           >  2:607330/1‑150 (MQ=255)
                                            aGTGGCGCTCATTTTGTGAGAACGCAGCGCCCGGCAATCCACTTGCATAGTAGTGGTGGTTTTGGTTGTCACTCTTAGTCATAACATCAGGCGTTTTGATGTGTATATCAATATTGATAACCTGGCTGACTATCATCGCCTTCGTGTCCt                                                                                                         <  2:1167469/150‑1 (MQ=255)
                                            aGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGGGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCt                                                                                                         <  2:233836/150‑1 (MQ=255)
                                             gTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGTCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTg                                                                                                        <  2:920602/150‑1 (MQ=255)
                                                   cTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCt                                                                                                  <  2:1325980/150‑1 (MQ=255)
                                                        tttGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACAcc                                                                                             <  1:1387782/150‑1 (MQ=255)
                                                         ttGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCt                                                                                            <  2:1101062/150‑1 (MQ=255)
                                                            tCAGAACGCAGCTCCCGGCAATCCACTGGGAGAGTAGGGGTTGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCTTATTGATGTTGATATAATTGTTGATAACCAGGCTGAAGATGATCGTCTTCGTGTCTTGGTCTATACAACtttt                                                                                         <  2:1214885/150‑1 (MQ=255)
                                                            tCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCtttt                                                                                         <  2:331860/150‑1 (MQ=255)
                                                            tCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGATGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCtttt                                                                                         <  2:1341217/150‑1 (MQ=255)
                                                            gcAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGGGGTGGTGGTTGTCACTCTTAGTCATAACAGCAGGCGGAGTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCtttt                                                                                         <  2:934911/149‑1 (MQ=255)
                                                             cAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTa                                                                                        <  1:1026438/150‑1 (MQ=255)
                                                               caaCGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAac                                                                                      >  1:1012995/2‑150 (MQ=255)
                                                                  cGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACAcc                                                                                   <  2:664152/150‑1 (MQ=255)
                                                                    cAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGa                                                                                 <  2:819320/150‑1 (MQ=255)
                                                                     aGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGAc                                                                                <  2:666041/150‑1 (MQ=255)
                                                                      gcgcCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAACAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACa                                                                               <  1:986558/150‑1 (MQ=255)
                                                                      gcgcCCGGCAATCCACTGGCAGAGTAGGGGTGGTGGTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACa                                                                               <  2:593917/150‑1 (MQ=255)
                                                                      gcgcCCGGCAATCCACTGGCAGAGTAGGGGGGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACa                                                                               <  2:1305831/150‑1 (MQ=255)
                                                                           cGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGtt                                                                          >  2:580144/1‑150 (MQ=255)
                                                                                  ccACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCGCTACACCTTTTATCACCGACAGTGTTCTCGCAg                                                                   >  1:945450/1‑150 (MQ=255)
                                                                                      tGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATAc                                                               >  1:1281836/1‑150 (MQ=255)
                                                                                       ggCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACg                                                              <  2:1082970/150‑1 (MQ=255)
                                                                                       ggCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACg                                                              >  2:463152/1‑150 (MQ=255)
                                                                                       ggCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACg                                                              >  1:1421089/1‑150 (MQ=255)
                                                                                       ggCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACg                                                              >  1:1088587/1‑150 (MQ=255)
                                                                                       ggCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACAACTTTTAACACCGACAGTGTTCTGGCAGATACg                                                              >  2:881035/1‑150 (MQ=255)
                                                                                               aGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTgg                                                      <  2:1088587/150‑1 (MQ=255)
                                                                                                ggggTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGt                                                     >  2:1050463/1‑150 (MQ=255)
                                                                                                  tgtggtgTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGtta                                                   >  1:1086789/2‑149 (MQ=255)
                                                                                                   gtggtgTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCg                                                  >  2:483803/1‑150 (MQ=255)
                                                                                                    tggtgTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGc                                                 <  1:470141/150‑1 (MQ=255)
                                                                                                      gtgtTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTa                                               <  1:1282553/150‑1 (MQ=255)
                                                                                                        gttggttgTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAAc                                             >  1:1348200/1‑150 (MQ=255)
                                                                                                                tCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCg                                     >  2:1281967/1‑150 (MQ=255)
                                                                                                                     ctTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAAc                                <  2:874150/150‑1 (MQ=255)
                                                                                                                     ctTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAAc                                <  1:126452/150‑1 (MQ=255)
                                                                                                                     ctTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGAGCATCGCCTTCGTGTCCTGGTCCCTTCACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAAc                                <  1:1316202/150‑1 (MQ=255)
                                                                                                                      ttAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACg                               <  1:607330/150‑1 (MQ=255)
                                                                                                                        aGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGATGACGATCATCGCCTTCGTGTCCTg                                                                                                        >  1:729121/1‑75 (MQ=255)
                                                                                                                        aGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGATGACGATCATCGCCTTCGTGTCCTg                                                                                                        <  2:729121/75‑1 (MQ=255)
                                                                                                                            aTAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTg                         <  1:1281967/150‑1 (MQ=255)
                                                                                                                             tAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCCGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGc                        <  2:1299062/150‑1 (MQ=255)
                                                                                                                             tAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGc                        <  1:719179/150‑1 (MQ=255)
                                                                                                                             tAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGc                        <  2:1039770/150‑1 (MQ=255)
                                                                                                                                    agGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGaac                 <  2:972953/150‑1 (MQ=255)
                                                                                                                                          aTTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTgg           <  1:106624/150‑1 (MQ=255)
                                                                                                                                           ttGATGTGGATATCAATCTTGATAACCAGGCTGACGCTCATCGCCATCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgc          >  1:985517/1‑150 (MQ=255)
                                                                                                                                           ttGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgc          <  2:379928/150‑1 (MQ=255)
                                                                                                                                           ttGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgc          >  2:1194233/1‑150 (MQ=255)
                                                                                                                                           ttGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTCGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgc          <  2:291765/150‑1 (MQ=255)
                                                                                                                                            cagcctgGTTATCAAGATTGATATCCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgca         >  1:480887/6‑150 (MQ=255)
                                                                                                                                            cagcctgGTTATCAAGATTGATATCCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgca         <  2:480887/145‑1 (MQ=255)
                                                                                                                                              gcctgGTTATCAAGATTGATATCCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGgcagc       >  1:984496/4‑150 (MQ=255)
                                                                                                                                                 tgGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCaag    >  1:60853/1‑150 (MQ=255)
                                                                                                                                                  ggATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCaaga   >  1:1368877/1‑150 (MQ=255)
                                                                                                                                                  ggATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCaaga   >  2:701485/1‑150 (MQ=255)
                                                                                                                                                   gATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCaagaa  >  2:614137/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAA  >  NC_000913/769643‑769939

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