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.

GCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTT  >  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)
                                                              aGCGCACAGACTGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTggcg                                                                                         >  2:911939/1‑147 (MQ=255)
                                                              aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                         >  1:6213/1‑150 (MQ=255)
                                                                cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAACAGCCCGCACCTGACA‑GTGCGGGc                                                                                       <  1:546117/150‑1 (MQ=255)
                                                                 gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACACGTAACGGTGCGGGCTGACGCGTACAGAAAACACAGAAAAAAGCCCGCCCCTGACC‑GTGCggtcg                                                                                      >  1:1240467/1‑147 (MQ=255)
                                                                 gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATAAGAACGACCATTACCACCACCATCACCATTAGCACAGGTAACGGTGCGGGCTTACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCt                                                                                      >  2:680629/1‑150 (MQ=255)
                                                                    cagacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCtttt                                                                                   >  2:1982060/1‑150 (MQ=255)
                                                                      gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCtttttt                                                                                 >  2:1477091/1‑150 (MQ=255)
                                                                        cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGTTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCtttttttt                                                                               <  2:40443/150‑1 (MQ=255)
                                                                        aagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCtttttttt                                                                               <  2:1827759/149‑1 (MQ=255)
                                                                          gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGCttttttctgt                                                                             <  1:1300974/150‑5 (MQ=38)
                                                                           aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCttttttttttc                                                                            <  1:673779/150‑3 (MQ=255)
                                                                           aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCg                                                                            >  2:1819911/1‑150 (MQ=255)
                                                                              aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGAcc                                                                         >  1:1639548/1‑150 (MQ=255)
                                                                              aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GAGCGGGCAATTTTTTACGAcc                                                                         >  1:155177/1‑150 (MQ=255)
                                                                                aaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCaa                                                                       <  1:672966/150‑1 (MQ=255)
                                                                                aaTTACAGAGTACACAACATCCAGGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCaa                                                                       <  1:201592/150‑1 (MQ=255)
                                                                                 aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCaaa                                                                      <  1:939270/150‑1 (MQ=255)
                                                                                  ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAg                                                                     >  2:1170421/1‑150 (MQ=255)
                                                                                   tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAgg                                                                    <  2:1172375/150‑1 (MQ=255)
                                                                                    aCAGAGTACACAACATCCATGAAACGCATTAGAACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCAGACGCGTACAGGAAACACAGAAAAAAGCCCGAACCTGACA‑GTGCGGGCTTTTTTGTTCGACCAAAGGt                                                                   >  2:1214648/1‑150 (MQ=255)
                                                                                     cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                                  <  2:534827/150‑1 (MQ=255)
                                                                                      agcgTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTCCGGGCTTTTTTTTTCGACCAAAGGTaa                                                                 >  2:625572/4‑150 (MQ=255)
                                                                                      agagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGCCCAAAGGTaa                                                                 >  2:1708209/1‑150 (MQ=255)
                                                                                      agagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                                 >  1:811434/1‑150 (MQ=255)
                                                                                           acacAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                            <  1:2011585/150‑1 (MQ=255)
                                                                                           acacAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCGGACGCGTACAGGAAACACAGAAAAACGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                            <  1:510222/150‑1 (MQ=255)
                                                                                           acacAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                            <  1:1554539/150‑1 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGTGCTTTTTTTTTCGACCAAAGGTAACGACGt                                                           >  1:631730/1‑150 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                           <  1:1437508/150‑1 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGt                                                           >  2:859199/1‑150 (MQ=255)
                                                                                            cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGt                                                           >  1:1055564/1‑150 (MQ=255)
                                                                                                 caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                      >  2:399103/1‑150 (MQ=255)
                                                                                                  aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                     >  2:642865/1‑150 (MQ=255)
                                                                                                    ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCa                                                   >  1:578252/1‑150 (MQ=255)
                                                                                                    ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                   >  1:698304/1‑150 (MQ=255)
                                                                                                         aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                              >  2:862019/1‑150 (MQ=255)
                                                                                                          aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                             >  1:1344337/1‑150 (MQ=255)
                                                                                                           aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                            <  1:232044/150‑1 (MQ=255)
                                                                                                            cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACA‑GTGCGGGCTTTATTTTTAGACCAAAAGTAACGAGGTAACAACCAAGCGactt                                           >  1:969597/1‑147 (MQ=255)
                                                                                                            cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                           <  2:1146205/150‑1 (MQ=255)
                                                                                                                   gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                    >  2:423148/1‑150 (MQ=255)
                                                                                                                    caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGTTGAAGtt                                   >  2:1680777/1‑150 (MQ=255)
                                                                                                                     accaccATTACCACCACCATAACCATTACCACAGGTAACGGTGCGGGATGACGCGTACAGGCAACACATAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                  <  1:1823926/150‑1 (MQ=255)
                                                                                                                       caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                >  1:1204439/1‑150 (MQ=255)
                                                                                                                       caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                >  2:1623205/1‑150 (MQ=255)
                                                                                                                        accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                               <  2:269081/150‑1 (MQ=255)
                                                                                                                             tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                          <  2:632373/150‑1 (MQ=255)
                                                                                                                               ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                      <  1:495257/136‑1 (MQ=255)
                                                                                                                               ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                        >  1:386919/1‑150 (MQ=255)
                                                                                                                               ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑ATGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                        >  1:792870/1‑150 (MQ=255)
                                                                                                                                 accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                            <  2:1877408/144‑1 (MQ=255)
                                                                                                                                 accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                            >  1:1877408/1‑144 (MQ=255)
                                                                                                                                    accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGCCCAAAGGTAACGAGGTACCAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                   >  2:982533/1‑150 (MQ=255)
                                                                                                                                    accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                   <  2:1199853/150‑1 (MQ=255)
                                                                                                                                    accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                   <  1:1819911/150‑1 (MQ=255)
                                                                                                                                     cccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                  >  2:729459/5‑150 (MQ=255)
                                                                                                                                     ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTTACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                  >  1:402029/1‑150 (MQ=255)
                                                                                                                                     ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                  >  1:518034/1‑150 (MQ=255)
                                                                                                                                      ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                 >  1:1304390/4‑150 (MQ=255)
                                                                                                                                      catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                 >  2:1682182/1‑150 (MQ=255)
                                                                                                                                      catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGc                 >  2:515966/1‑150 (MQ=255)
                                                                                                                                       cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                >  1:31068/3‑150 (MQ=255)
                                                                                                                                       cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                >  1:1484349/3‑150 (MQ=255)
                                                                                                                                       atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                <  2:1790669/150‑1 (MQ=255)
                                                                                                                                       atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                >  1:1701131/1‑150 (MQ=255)
                                                                                                                                          accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt             >  2:2012623/1‑150 (MQ=255)
                                                                                                                                          accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAACCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCCACCAAAGGTAACGAGGTAACAACCATGCGAGTGTAGAAGTTAGGCGGTACATCAGTGGCAAAt             >  2:640706/1‑150 (MQ=255)
                                                                                                                                                aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa       >  2:880980/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc      >  2:1365108/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc      >  2:1039319/1‑150 (MQ=255)
                                                                                                                                                  cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAACAACGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGTGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACg     <  2:1557100/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTT  >  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: ‑