Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,876 C→A 5.2% R163S (AGG→AGT ykgH ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,8760CA 5.2% 377.8 / 7.9 134R163S (AGG→AGTykgHputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (62/65);  new base A (4/3);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 137 total aligned reads displayed.

CCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGC  >  NC_000913/324728‑325022
                                                                                                                                                    |                                                                                                                                                  
ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                   <  2:686961/150‑1 (MQ=255)
ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                   >  1:803637/1‑150 (MQ=255)
accattttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                   >  2:1175364/5‑150 (MQ=255)
 cATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCt                                                                                                                                                  >  2:1494582/1‑150 (MQ=255)
   tttttttGACAAGGCAAAACATTACAGAAATTACAATTCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTgg                                                                                                                                                <  2:1451730/150‑1 (MQ=255)
     tttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTATGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGta                                                                                                                                              >  1:203263/1‑149 (MQ=255)
      ttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCa                                                                                                                                             >  1:744133/1‑150 (MQ=255)
      ttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCa                                                                                                                                             >  2:1622163/1‑150 (MQ=255)
       tttGACAAGTCAAAACATTACAGAAATTACAATTCTTAGAAATAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAAAGACCAGTTCCTGGTCAg                                                                                                                                            <  1:327852/150‑1 (MQ=255)
        ttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGAACTGGTCatt                                                                                                                                           <  1:1133392/150‑3 (MQ=255)
        ttGACAAGGCAAAACATTACAGAAATTACAAAGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGAACTGGTCatt                                                                                                                                           <  1:1532426/150‑3 (MQ=255)
         tGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGt                                                                                                                                          <  2:616784/150‑1 (MQ=255)
          gACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGta                                                                                                                                         >  1:459532/1‑150 (MQ=255)
          gACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGta                                                                                                                                         <  2:803637/150‑1 (MQ=255)
           aCAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGCTATTTATTTTTTATATCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGtat                                                                                                                                        >  1:1749625/1‑150 (MQ=255)
            cAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGtata                                                                                                                                       >  2:143023/1‑150 (MQ=255)
                gCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAt                                                                                                                                   >  1:849206/1‑150 (MQ=255)
                 cAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAtt                                                                                                                                  >  2:1469171/1‑150 (MQ=255)
                     aaaTTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTg                                                                                                                              <  2:1678306/148‑1 (MQ=255)
                      cATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGt                                                                                                                             <  1:1269894/150‑1 (MQ=255)
                       aTTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGta                                                                                                                            <  1:1156328/150‑1 (MQ=255)
                       aTTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGta                                                                                                                            <  1:1703271/150‑1 (MQ=255)
                       aTTACAGAAATAACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGTGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGta                                                                                                                            <  2:937217/150‑1 (MQ=255)
                          aCAGAAATTAGAATGCGTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaa                                                                                                                         >  2:47085/1‑150 (MQ=255)
                          aCAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaa                                                                                                                         >  2:1561511/1‑150 (MQ=255)
                          aCAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGCGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGTCCTGGTCAGGTATATCATTTTTGTATaa                                                                                                                         >  2:1603375/1‑150 (MQ=255)
                           cAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaaa                                                                                                                        >  1:1316286/1‑150 (MQ=255)
                           cAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaaa                                                                                                                        >  2:297773/1‑150 (MQ=255)
                             gAAATTACAATGCTTAGATAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGc                                                                                                                      >  2:542011/1‑150 (MQ=255)
                               aaTTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTc                                                                                                                    >  2:1302501/1‑150 (MQ=255)
                                     aaTGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTTAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGtt                                                                                                              <  2:1207800/150‑1 (MQ=255)
                                     aaTGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGtt                                                                                                              <  2:513769/150‑1 (MQ=255)
                                     aaTGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTAATTTTTATTTCTGGTGAGCTTAAAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGATCACTGtt                                                                                                              >  1:1576917/1‑150 (MQ=255)
                                      aTGATTAGAAAAAATACCAGAGCCAGTCAAATGTCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTa                                                                                                             <  1:1123908/150‑1 (MQ=255)
                                            aGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCt                                                                                                       <  1:28344/150‑1 (MQ=255)
                                            aGAAAAAATCCCAGAGCCAGGCAAATGGCAGCACGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTAAATCATTTTTGTATAAAGCTCACTGTTAATTGCt                                                                                                       >  2:853361/1‑150 (MQ=255)
                                             gAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTg                                                                                                      >  2:448487/1‑150 (MQ=255)
                                             gAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGATCGGATACAATGACCAGTTCCTGGTCAGGTATATCAATTATGTATAAAGCTCACTGTTAATTGCTg                                                                                                      >  1:726271/1‑150 (MQ=255)
                                                 aaaTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTTAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAgg                                                                                                  <  1:571449/150‑1 (MQ=255)
                                                   aTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCg                                                                                                <  1:1797888/150‑1 (MQ=255)
                                                   aTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACACCCATCACCCATCGGATACAAAGACCACTTCCTGGTCAGGTATATCATTTTAGTATAAAGCTCACTGTAAATTTCTGAAGGCg                                                                                                >  1:1131629/1‑150 (MQ=255)
                                                     cccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                              >  1:1361421/1‑150 (MQ=255)
                                                     cccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                              >  2:161616/1‑150 (MQ=255)
                                                     cagagagCCAGGCAAATGGGAGCAGGGTATTTATTTTTTATTTCTGGTGAGATTATAGTTACAGAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                              <  2:806605/147‑1 (MQ=255)
                                                      ccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGttt                                                                                             >  2:1560337/1‑150 (MQ=255)
                                                       cAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTg                                                                                            <  1:1004535/150‑1 (MQ=255)
                                                       cAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTg                                                                                            >  2:1644283/1‑150 (MQ=255)
                                                          atccAGGCAAATTGCAGCAGTGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAt                                                                                         <  1:1521357/148‑1 (MQ=255)
                                                           gCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAtc                                                                                        <  2:1654636/150‑1 (MQ=255)
                                                            ccAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAtct                                                                                       >  1:685489/1‑150 (MQ=255)
                                                            ccAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAtct                                                                                       <  2:481072/150‑1 (MQ=255)
                                                             cAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGGTCGGATACAATGACCAGTGCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAtctc                                                                                      >  2:714519/1‑150 (MQ=255)
                                                               ggCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCAt                                                                                    >  1:371797/1‑150 (MQ=255)
                                                                gCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATc                                                                                   <  2:567168/150‑1 (MQ=255)
                                                                      ggCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAAc                                                                             <  2:654518/150‑1 (MQ=255)
                                                                        cagcagGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACtt                                                                           >  1:462209/1‑150 (MQ=255)
                                                                         agcagGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtg                                                                          <  1:1302501/150‑1 (MQ=255)
                                                                             gggTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgcctg                                                                      <  2:700227/150‑1 (MQ=255)
                                                                                 atttatttTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCAt                                                                  >  1:1172554/1‑150 (MQ=255)
                                                                                 atttatttTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCAt                                                                  >  2:981274/1‑150 (MQ=255)
                                                                                 atttatttTTTATTTCTGCTGAGCTTATAGTTACACAGCCATCACCGATCGGATACAATGACAAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCAt                                                                  <  1:1051593/150‑1 (MQ=255)
                                                                                      ttttttatttCTGGTGAGCTTATAGTGACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGa                                                             <  1:931053/150‑1 (MQ=255)
                                                                                       tttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGat                                                            >  1:1410876/1‑150 (MQ=255)
                                                                                        ttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGata                                                           >  2:134508/1‑150 (MQ=255)
                                                                                         tttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCAATTATGatat                                                          >  2:1747986/1‑150 (MQ=255)
                                                                                           tatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTc                                                        >  1:717857/1‑150 (MQ=255)
                                                                                                cTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACg                                                   <  1:251529/150‑1 (MQ=255)
                                                                                                     gAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGTTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATATCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATAc                                              >  2:671083/1‑150 (MQ=255)
                                                                                                     gAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATAc                                              <  1:1486731/150‑1 (MQ=255)
                                                                                                             aGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGat                                      >  2:949961/1‑150 (MQ=255)
                                                                                                             aGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGat                                      <  1:195504/150‑1 (MQ=255)
                                                                                                              gTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGCTATTCTGCCGCATACCATTTGata                                     >  1:134067/1‑150 (MQ=255)
                                                                                                              gTTACACAGCCATCACCGTTCGGATACAAAGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGATAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGata                                     <  2:914409/150‑1 (MQ=255)
                                                                                                                tACACCTCCATCACCGTTCGGATACAAAGACCAGTTCCTGGGCAGGTATATCATTTTTGTATAAAGCTCACTTTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGataat                                   <  2:1105828/150‑1 (MQ=255)
                                                                                                                tACACAGCCATCACCTTTCGGATACAATGACCAGTTCCTGGGCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGataat                                   <  1:1197028/150‑1 (MQ=255)
                                                                                                                  cacaGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTAACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTATTACTCATACCATTTTATAATag                                 >  2:312371/1‑150 (MQ=255)
                                                                                                                   acaGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATaga                                <  2:1160655/150‑1 (MQ=255)
                                                                                                                       ccATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGcc                            <  1:1279324/150‑1 (MQ=255)
                                                                                                                         aTCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATATAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTg                          <  1:1053611/150‑1 (MQ=255)
                                                                                                                           cACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCa                        <  1:427254/150‑1 (MQ=255)
                                                                                                                            aCCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAGCTTGCCTGACATTATGATTTTCTGACGCATACCATTTGATAATAGAGGCCTGCAc                       >  1:347849/1‑150 (MQ=255)
                                                                                                                              cGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCAcgc                     >  1:1792732/1‑150 (MQ=255)
                                                                                                                               gTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCAcgcg                    <  2:444681/150‑1 (MQ=255)
                                                                                                                                 tCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTTTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCgag                  <  1:143023/150‑1 (MQ=255)
                                                                                                                                 tCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCgag                  <  2:744133/150‑1 (MQ=255)
                                                                                                                                 tCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCgag                  <  2:1423934/150‑1 (MQ=255)
                                                                                                                                    gATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCACACCATTTGATAATAGAGGCCTGCACGCGAGata               <  2:202758/150‑1 (MQ=255)
                                                                                                                                      tACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATAtt             >  1:614375/1‑150 (MQ=255)
                                                                                                                                        cAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTc           >  2:1101965/1‑150 (MQ=255)
                                                                                                                                        cAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTc           >  2:1743641/1‑150 (MQ=255)
                                                                                                                                        cAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTc           <  1:134888/150‑1 (MQ=255)
                                                                                                                                         cctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          >  2:994782/3‑150 (MQ=255)
                                                                                                                                         caTGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGATCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          >  2:1359619/2‑150 (MQ=255)
                                                                                                                                         aaTTACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACTCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          <  1:714872/150‑1 (MQ=255)
                                                                                                                                         aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAACGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          <  1:85380/150‑1 (MQ=255)
                                                                                                                                         aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTTCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          >  1:892481/1‑150 (MQ=255)
                                                                                                                                             aCCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCtt      <  2:459532/150‑1 (MQ=255)
                                                                                                                                               cAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATTATAGAGGCCTGCACGCGAGATATTTCCTCTTca    >  1:321051/1‑150 (MQ=255)
                                                                                                                                                 gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagc  <  2:685489/150‑1 (MQ=255)
                                                                                                                                                 gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagc  >  1:36874/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGC  >  NC_000913/324728‑325022

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