Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 237 A→T 5.9% T16T (ACA→ACT) ‡ thrL → thr operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132370AT 5.9% 330.8 / 9.1 119T16T (ACA→ACT) ‡thrLthr operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (73/38);  new base T (4/3);  total (77/42)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 122 total aligned reads displayed.

GCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTT  >  NC_000913/88‑386
                                                                                                                                                     |                                                                                                                                                      
gCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACcaca                                                                                                                                                        <  1:188834/150‑1 (MQ=255)
  cGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGACATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAgg                                                                                                                                                      <  1:1681077/150‑1 (MQ=255)
        tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACgg                                                                                                                                                <  2:1320402/150‑1 (MQ=255)
        tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACACATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGCTAACgg                                                                                                                                                >  2:1729160/1‑150 (MQ=255)
           aTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCTGTGGTAATGGTga                                                                                                                                             <  1:1028031/150‑2 (MQ=255)
            taaagtatttAGTGACCTAAGTCAATAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCAAGAAAACCAATACCACCCCCATTACCACCACCATCACCATTAACACAGGTAACGGGGCg                                                                                                                                            >  2:478299/8‑150 (MQ=255)
               aaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTCCCACAGGTAACGGTGCGGGc                                                                                                                                         <  1:1993909/150‑1 (MQ=255)
                gtatttAGTGCCCTAAGTCAATAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCt                                                                                                                                        >  1:1472146/4‑150 (MQ=255)
                 aTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTg                                                                                                                                       >  1:632373/1‑150 (MQ=255)
                    ttATTGACTTAGGTCACTAAATACTTTATCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGAcg                                                                                                                                    <  2:1369790/150‑1 (MQ=255)
                      aTTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGAcgcg                                                                                                                                  >  1:34095/1‑150 (MQ=255)
                          aCTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACa                                                                                                                              >  1:228716/1‑150 (MQ=255)
                          aCTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACa                                                                                                                              <  2:905601/150‑1 (MQ=255)
                                gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                        <  2:1173387/150‑1 (MQ=255)
                                gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                        <  1:1108579/150‑1 (MQ=255)
                                 tCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAaca                                                                                                                       >  2:712170/1‑150 (MQ=255)
                                  cACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAacac                                                                                                                      >  2:1437508/1‑150 (MQ=255)
                                  cACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAacac                                                                                                                      >  2:503349/1‑150 (MQ=255)
                                      aaataaatTAACCAATATAGGCATAGCGCAAAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGAAATAGCACCACCAATACCACCAACATCACCATAACCAGAGGTAACGGTGCGGGCTGCCGGGTACTCGAACCACAGaa                                                                                                                  >  1:976637/1‑150 (MQ=255)
                                       aaTACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaaa                                                                                                                 <  2:929758/150‑1 (MQ=255)
                                           cTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACAACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAg                                                                                                             >  1:1296911/1‑150 (MQ=255)
                                             ttAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGcc                                                                                                           >  1:1336704/1‑150 (MQ=255)
                                                aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        >  2:1463295/1‑150 (MQ=255)
                                                aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        >  1:389417/1‑150 (MQ=255)
                                                aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        >  2:1373119/1‑150 (MQ=255)
                                                aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGAATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGTAAAAAGCCCGc                                                                                                        >  2:524535/1‑150 (MQ=255)
                                                    atatAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACGATTCCCACCACCATCACCATTACCCCAGGTTACGGTGCGGGCGGCCGCGGACCGGAAAGCCAGAAAAAAGCCCGCAccg                                                                                                    >  2:335522/1‑149 (MQ=255)
                                                       taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACAACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                 >  2:276055/1‑150 (MQ=255)
                                                         ggCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGAGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCACTGACac                                                                                               >  1:1506743/1‑149 (MQ=255)
                                                              aGCGCACAGACTGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTggcg                                                                                         >  2:911939/1‑147 (MQ=255)
                                                              aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                         >  1:6213/1‑150 (MQ=255)
                                                                cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAACAGCCCGCACCTGACAGTGCGGGc                                                                                       <  1:546117/150‑1 (MQ=255)
                                                                 gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACACGTAACGGTGCGGGCTGACGCGTACAGAAAACACAGAAAAAAGCCCGCCCCTGACCGTGCggtcg                                                                                      >  1:1240467/1‑147 (MQ=255)
                                                                 gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATAAGAACGACCATTACCACCACCATCACCATTAGCACAGGTAACGGTGCGGGCTTACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCt                                                                                      >  2:680629/1‑150 (MQ=255)
                                                                    cagacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttt                                                                                   >  2:1982060/1‑150 (MQ=255)
                                                                      gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttttt                                                                                 >  2:1477091/1‑150 (MQ=255)
                                                                        cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGTTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttttttt                                                                               <  2:40443/150‑1 (MQ=255)
                                                                        aagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttttttt                                                                               <  2:1827759/149‑1 (MQ=255)
                                                                          gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCttttttctgt                                                                             <  1:1300974/150‑5 (MQ=38)
                                                                           aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttttc                                                                            <  1:673779/150‑3 (MQ=255)
                                                                           aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCg                                                                            >  2:1819911/1‑150 (MQ=255)
                                                                              aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAcc                                                                         >  1:1639548/1‑150 (MQ=255)
                                                                              aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGAGCGGGCAATTTTTTACGAcc                                                                         >  1:155177/1‑150 (MQ=255)
                                                                                aaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCaa                                                                       <  1:672966/150‑1 (MQ=255)
                                                                                aaTTACAGAGTACACAACATCCAGGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCaa                                                                       <  1:201592/150‑1 (MQ=255)
                                                                                 aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCaaa                                                                      <  1:939270/150‑1 (MQ=255)
                                                                                  ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAg                                                                     >  2:1170421/1‑150 (MQ=255)
                                                                                   tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAgg                                                                    <  2:1172375/150‑1 (MQ=255)
                                                                                    aCAGAGTACACAACATCCATGAAACGCATTAGAACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCAGACGCGTACAGGAAACACAGAAAAAAGCCCGAACCTGACAGTGCGGGCTTTTTTGTTCGACCAAAGGt                                                                   >  2:1214648/1‑150 (MQ=255)
                                                                                     cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                                  <  2:534827/150‑1 (MQ=255)
                                                                                      agcgTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTCCGGGCTTTTTTTTTCGACCAAAGGTaa                                                                 >  2:625572/4‑150 (MQ=255)
                                                                                      agagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGCCCAAAGGTaa                                                                 >  2:1708209/1‑150 (MQ=255)
                                                                                      agagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                                 >  1:811434/1‑150 (MQ=255)
                                                                                           acacAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                            <  1:2011585/150‑1 (MQ=255)
                                                                                           acacAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCGGACGCGTACAGGAAACACAGAAAAACGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                            <  1:510222/150‑1 (MQ=255)
                                                                                           acacAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                            <  1:1554539/150‑1 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGTGCTTTTTTTTTCGACCAAAGGTAACGACGt                                                           >  1:631730/1‑150 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                           <  1:1437508/150‑1 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGt                                                           >  2:859199/1‑150 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGt                                                           >  1:1055564/1‑150 (MQ=255)
                                                                                                 caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                      >  2:399103/1‑150 (MQ=255)
                                                                                                  aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                     >  2:642865/1‑150 (MQ=255)
                                                                                                    ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCa                                                   >  1:578252/1‑150 (MQ=255)
                                                                                                    ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                   >  1:698304/1‑150 (MQ=255)
                                                                                                         aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                              >  2:862019/1‑150 (MQ=255)
                                                                                                          aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                             >  1:1344337/1‑150 (MQ=255)
                                                                                                           aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                            <  1:232044/150‑1 (MQ=255)
                                                                                                            cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTATTTTTAGACCAAAAGTAACGAGGTAACAACCAAGCGactt                                           >  1:969597/1‑147 (MQ=255)
                                                                                                            cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                           <  2:1146205/150‑1 (MQ=255)
                                                                                                                   gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                    >  2:423148/1‑150 (MQ=255)
                                                                                                                    caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGTTGAAGtt                                   >  2:1680777/1‑150 (MQ=255)
                                                                                                                     accaccATTACCACCACCATAACCATTACCACAGGTAACGGTGCGGGATGACGCGTACAGGCAACACATAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                  <  1:1823926/150‑1 (MQ=255)
                                                                                                                       caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                >  1:1204439/1‑150 (MQ=255)
                                                                                                                       caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                >  2:1623205/1‑150 (MQ=255)
                                                                                                                        accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                               <  2:269081/150‑1 (MQ=255)
                                                                                                                             tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                          <  2:632373/150‑1 (MQ=255)
                                                                                                                               ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                      <  1:495257/136‑1 (MQ=255)
                                                                                                                               ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                        >  1:386919/1‑150 (MQ=255)
                                                                                                                               ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAATGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                        >  1:792870/1‑150 (MQ=255)
                                                                                                                                 accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                            <  2:1877408/144‑1 (MQ=255)
                                                                                                                                 accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                            >  1:1877408/1‑144 (MQ=255)
                                                                                                                                    accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGCCCAAAGGTAACGAGGTACCAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                   >  2:982533/1‑150 (MQ=255)
                                                                                                                                    accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                   <  2:1199853/150‑1 (MQ=255)
                                                                                                                                    accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                   <  1:1819911/150‑1 (MQ=255)
                                                                                                                                     cccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                  >  2:729459/5‑150 (MQ=255)
                                                                                                                                     ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTTACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                  >  1:402029/1‑150 (MQ=255)
                                                                                                                                     ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                  >  1:518034/1‑150 (MQ=255)
                                                                                                                                      ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                 >  1:1304390/4‑150 (MQ=255)
                                                                                                                                      catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                 >  2:1682182/1‑150 (MQ=255)
                                                                                                                                      catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGc                 >  2:515966/1‑150 (MQ=255)
                                                                                                                                       cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                >  1:31068/3‑150 (MQ=255)
                                                                                                                                       cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                >  1:1484349/3‑150 (MQ=255)
                                                                                                                                       atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                <  2:1790669/150‑1 (MQ=255)
                                                                                                                                       atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                >  1:1701131/1‑150 (MQ=255)
                                                                                                                                          accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt             >  2:2012623/1‑150 (MQ=255)
                                                                                                                                          accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAACCCCGCACCTGACAGTGCGGGCTTTTTTTTTCCACCAAAGGTAACGAGGTAACAACCATGCGAGTGTAGAAGTTAGGCGGTACATCAGTGGCAAAt             >  2:640706/1‑150 (MQ=255)
                                                                                                                                                aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa       >  2:880980/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc      >  2:1365108/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc      >  2:1039319/1‑150 (MQ=255)
                                                                                                                                                  cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAACAACGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGTGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACg     <  2:1557100/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTT  >  NC_000913/88‑386

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: