Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 243 C→T 5.6% N18N (AAC→AAT)  thrL → thr operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132430CT 5.6% 400.3 / 7.8 125N18N (AAC→AAT) thrLthr operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (75/43);  new base T (4/3);  total (79/46)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 127 total aligned reads displayed.

AGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCG  >  NC_000913/94‑392
                                                                                                                                                     |                                                                                                                                                       
aGTAAATTAAAATTTCATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAAc                                                                                                                                                         >  2:629657/1‑150 (MQ=255)
  tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACgg                                                                                                                                                       <  2:1320402/150‑1 (MQ=255)
  tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACACATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGCTAACgg                                                                                                                                                       >  2:1729160/1‑150 (MQ=255)
    aaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTg                                                                                                                                                     <  1:1407657/150‑1 (MQ=255)
     aTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCTGTGGTAATGGTga                                                                                                                                                    <  1:1028031/150‑2 (MQ=255)
         aaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATTGCGCACAGACAGATAAAAATTACAGAGTACACAACAACCATGTAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGTgatggt                                                                                                                                                <  2:104628/150‑6 (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)
                    aCTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACa                                                                                                                                     <  1:928798/150‑1 (MQ=255)
                     cTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAg                                                                                                                                    >  1:1289116/1‑150 (MQ=255)
                          gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                               <  1:1108579/150‑1 (MQ=255)
                          gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                               <  2:1173387/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)
                                 aaTACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaaa                                                                                                                        <  2:929758/150‑1 (MQ=255)
                                     cTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACAACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAg                                                                                                                    >  1:1296911/1‑150 (MQ=255)
                                       ttAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGcc                                                                                                                  >  1:1336704/1‑150 (MQ=255)
                                          aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAATGCCCGc                                                                                                               >  2:693107/1‑150 (MQ=255)
                                          aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                               >  2:1373119/1‑150 (MQ=255)
                                          aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                               >  2:1463295/1‑150 (MQ=255)
                                          aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                               >  1:389417/1‑150 (MQ=255)
                                          aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACAACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                               >  2:894142/1‑150 (MQ=255)
                                          aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGAATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGTAAAAAGCCCGc                                                                                                               >  2:524535/1‑150 (MQ=255)
                                                 taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACAACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                        >  2:276055/1‑150 (MQ=255)
                                                  aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                       >  2:1758409/1‑150 (MQ=255)
                                                    gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAGTAGCACCCCCATTACCACCACCATCACCATTACCACAGGTAACGGAGCGGGCTGCCGCGTACAGGACACACAGAAAAAAGCCCGCCccagcca‑gg                                                                                                    >  1:394141/1‑149 (MQ=255)
                                                        aGCGCACAGACTGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTggcg                                                                                                >  2:911939/1‑147 (MQ=255)
                                                        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                >  1:6213/1‑150 (MQ=255)
                                                         gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                               <  2:1209034/150‑1 (MQ=255)
                                                          cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAACAGCCCGCACCTGACA‑GTGCGGGc                                                                                              <  1:546117/150‑1 (MQ=255)
                                                           gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATAAGAACGACCATTACCACCACCATCACCATTAGCACAGGTAACGGTGCGGGCTTACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                            >  2:680629/1‑150 (MQ=255)
                                                              cagacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttt                                                                                         >  2:1982060/1‑150 (MQ=255)
                                                                gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttttt                                                                                       >  2:1477091/1‑150 (MQ=255)
                                                                  cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGTTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttttttt                                                                                     <  2:40443/150‑1 (MQ=255)
                                                                  aagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttttttt                                                                                     <  2:1827759/149‑1 (MQ=255)
                                                                    gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttttttctgt                                                                                   <  1:1300974/150‑5 (MQ=38)
                                                                     aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttttc                                                                                  <  1:673779/150‑3 (MQ=255)
                                                                       ataaaTTACAGAGTACACTACATCCATGAAACGCATTAGCACCACCATTACCACCTCCATCTCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAACCACAGAAAACAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAc                                                                                <  2:375420/148‑1 (MQ=255)
                                                                        aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAcc                                                                               >  1:1639548/1‑150 (MQ=255)
                                                                        aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GAGCGGGC‑AATTTTTTACGAcc                                                                               >  1:155177/1‑150 (MQ=255)
                                                                         aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCTGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCa                                                                              >  1:242959/1‑150 (MQ=255)
                                                                          aaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaa                                                                             <  1:672966/150‑1 (MQ=255)
                                                                          aaTTACAGAGTACACAACATCCAGGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaa                                                                             <  1:201592/150‑1 (MQ=255)
                                                                           aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaaa                                                                            <  1:939270/150‑1 (MQ=255)
                                                                             tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAgg                                                                          <  2:1172375/150‑1 (MQ=255)
                                                                                agcgTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTCCGGGC‑TTTTTTTTTCGACCAAAGGTaa                                                                       >  2:625572/4‑150 (MQ=255)
                                                                                agagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGCCCAAAGGTaa                                                                       >  2:1708209/1‑150 (MQ=255)
                                                                                 cagaaaacaacaTCCATGAAACGCATTACCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAc                                                                      <  1:1752826/144‑1 (MQ=255)
                                                                                     acacAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAgg                                                                  <  1:2011585/150‑1 (MQ=255)
                                                                                     acacAACATCCATGAAACGCATTCGCACCACCATTACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCGGACGCGTACAGGAAACACAGAAAAACGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAgg                                                                  <  1:510222/150‑1 (MQ=255)
                                                                                     acacAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAgg                                                                  <  1:1554539/150‑1 (MQ=255)
                                                                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGTGC‑TTTTTTTTTCGACCAAAGGTAACGACGt                                                                 >  1:631730/1‑150 (MQ=255)
                                                                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                 <  1:1437508/150‑1 (MQ=255)
                                                                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGt                                                                 >  2:859199/1‑150 (MQ=255)
                                                                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGt                                                                 >  1:1055564/1‑150 (MQ=255)
                                                                                           caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                            >  2:399103/1‑150 (MQ=255)
                                                                                            aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                           >  2:642865/1‑150 (MQ=255)
                                                                                              ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                         >  1:698304/1‑150 (MQ=255)
                                                                                                   aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                    >  2:862019/1‑150 (MQ=255)
                                                                                                    aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                   >  1:1344337/1‑150 (MQ=255)
                                                                                                     aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                  <  1:232044/150‑1 (MQ=255)
                                                                                                      cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTATTTTTAGACCAAAAGTAACGAGGTAACAACCAAGCGactt                                                 >  1:969597/1‑147 (MQ=255)
                                                                                                      cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                                 <  2:1146205/150‑1 (MQ=255)
                                                                                                             gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                          >  2:423148/1‑150 (MQ=255)
                                                                                                              caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGTTGAAGtt                                         >  2:1680777/1‑150 (MQ=255)
                                                                                                               accaccATTACCACCACCATAACCATTACCACAGGTAACGGTGCGGGATGACGCGTACAGGCAACACATAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                        <  1:1823926/150‑1 (MQ=255)
                                                                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                      >  2:1623205/1‑150 (MQ=255)
                                                                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                      >  1:1204439/1‑150 (MQ=255)
                                                                                                                  accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                     <  2:269081/150‑1 (MQ=255)
                                                                                                                  accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                     >  2:1561131/1‑150 (MQ=255)
                                                                                                                       tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                                <  2:632373/150‑1 (MQ=255)
                                                                                                                        accaccaccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                               >  1:1538310/1‑150 (MQ=255)
                                                                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                            <  1:495257/136‑1 (MQ=255)
                                                                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                            >  2:495257/1‑136 (MQ=255)
                                                                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑ATGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                              >  1:792870/1‑150 (MQ=255)
                                                                                                                           accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                                  <  2:1877408/144‑1 (MQ=255)
                                                                                                                           accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                                  >  1:1877408/1‑144 (MQ=255)
                                                                                                                              accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGCCCAAAGGTAACGAGGTACCAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                         >  2:982533/1‑150 (MQ=255)
                                                                                                                              accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                         <  2:1199853/150‑1 (MQ=255)
                                                                                                                              accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                         <  1:1819911/150‑1 (MQ=255)
                                                                                                                               ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTTACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                        >  1:402029/1‑150 (MQ=255)
                                                                                                                               ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                        >  1:518034/1‑150 (MQ=255)
                                                                                                                                ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                       >  1:1304390/4‑150 (MQ=255)
                                                                                                                                catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGc                       >  2:515966/1‑150 (MQ=255)
                                                                                                                                 cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                      >  1:31068/3‑150 (MQ=255)
                                                                                                                                 cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                      >  1:1484349/3‑150 (MQ=255)
                                                                                                                                 atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                      <  2:1790669/150‑1 (MQ=255)
                                                                                                                                    accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCCACCAAAGGTAACGAGGTAACCACCATTCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt                   >  1:55633/1‑150 (MQ=255)
                                                                                                                                          aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa             >  2:880980/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc            >  2:1365108/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc            >  2:1039319/1‑150 (MQ=255)
                                                                                                                                               aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt        >  2:192709/1‑150 (MQ=255)
                                                                                                                                                ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGtttt       >  1:1491754/1‑150 (MQ=255)
                                                                                                                                                  tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc      <  1:894142/150‑1 (MQ=255)
                                                                                                                                                   aaCGGTGCGGGCTGACGCGTACAGGAAACTCAGAAAAACGCCCGCACCTGACA‑TTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg    <  2:1667031/150‑1 (MQ=255)
                                                                                                                                                     cGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg  <  2:285446/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCG  >  NC_000913/94‑392

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: ‑