Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,140,400 G→A 5.0% R405W (CGG→TGG)  asmA ← putative assembly protein AsmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,140,4000GA 5.0% 319.1 / 2.1 98R405W (CGG→TGG) asmAputative assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base G (40/53);  new base A (2/3);  total (42/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTCA  >  NC_000913/2140251‑2140549
                                                                                                                                                     |                                                                                                                                                     
gACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTATACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCg                                                                                                                                                       <  2:661830/150‑1 (MQ=255)
gACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCg                                                                                                                                                       >  2:895921/1‑150 (MQ=255)
   aTGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTgg                                                                                                                                                    <  1:1567822/150‑1 (MQ=255)
     gtgCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCt                                                                                                                                                  <  2:1216138/150‑1 (MQ=255)
      tgCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTg                                                                                                                                                 >  2:1288282/1‑150 (MQ=255)
       gCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTgg                                                                                                                                                <  2:649833/150‑1 (MQ=255)
       gCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTgg                                                                                                                                                >  2:1610683/1‑150 (MQ=255)
       gCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTgg                                                                                                                                                <  2:1545422/150‑1 (MQ=255)
        cgtGTCCTTGCGAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACGTGATAAATCACCAGCCAGTGACATTTTTCAGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGTCTGGa                                                                                                                                               <  2:817446/148‑1 (MQ=255)
        ccTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGa                                                                                                                                               <  2:1283356/150‑1 (MQ=255)
        ccTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGa                                                                                                                                               <  1:1242815/150‑1 (MQ=255)
          tGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAg                                                                                                                                             <  1:949736/150‑1 (MQ=255)
          tGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAg                                                                                                                                             >  1:474232/1‑150 (MQ=255)
           gTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCt                                                                                                                                                  <  1:906226/144‑1 (MQ=255)
           gTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCt                                                                                                                                                  >  2:906226/1‑144 (MQ=255)
             ccTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGttt                                                                                                                                          <  1:1554639/150‑1 (MQ=255)
              cTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTa                                                                                                                                         <  2:1019310/150‑1 (MQ=255)
                tGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATc                                                                                                                                       >  1:919850/1‑150 (MQ=255)
                tGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacgttg                                                                                                                                       <  1:146374/150‑8 (MQ=255)
                      ttGTGGCGGATTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGTAATCGGATAGTTAAAGGCCTTCAGAATGGTACCAATCTCATCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGa                                                                                                                                 <  1:907818/150‑1 (MQ=255)
                        gtgGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGAtt                                                                                                                               <  1:422160/150‑1 (MQ=255)
                        gtgGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGAtt                                                                                                                               >  1:225918/1‑150 (MQ=255)
                           gCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAt                                                                                                                            >  1:508079/1‑150 (MQ=255)
                            cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAtt                                                                                                                           >  1:681939/1‑150 (MQ=255)
                            cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAtt                                                                                                                           <  2:1143362/150‑1 (MQ=255)
                             ggAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTg                                                                                                                          <  2:797670/150‑1 (MQ=255)
                             ggAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTg                                                                                                                          <  1:446614/150‑1 (MQ=255)
                              gAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGa                                                                                                                         <  1:1160283/150‑1 (MQ=255)
                                 tGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCAGCGGATTTATTGATGt                                                                                                                      >  1:1128275/1‑150 (MQ=255)
                                   cGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcg                                                                                                                    <  2:1351149/150‑1 (MQ=255)
                                   cGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcg                                                                                                                    <  2:225918/150‑1 (MQ=255)
                                   cGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcg                                                                                                                    <  1:985058/150‑1 (MQ=255)
                                   cGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcg                                                                                                                    >  2:1088836/1‑150 (MQ=255)
                                   cGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcg                                                                                                                    >  2:1394779/1‑150 (MQ=255)
                                       ggCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTc                                                                                                                >  1:1344459/1‑150 (MQ=255)
                                        gCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTcc                                                                                                               >  1:162150/1‑150 (MQ=255)
                                         cATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCa                                                                                                              <  1:611921/150‑1 (MQ=255)
                                          aTCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAg                                                                                                             <  2:1529910/150‑1 (MQ=255)
                                            cTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCg                                                                                                           <  2:1344459/150‑1 (MQ=255)
                                              aTGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTg                                                                                                         >  1:1440254/1‑150 (MQ=255)
                                               tGTCAGCACCGGAGAAATTACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGc                                                                                                        >  1:1096681/1‑150 (MQ=255)
                                                  cAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCg                                                                                                     <  2:919850/150‑1 (MQ=255)
                                                     cACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCa                                                                                                  <  1:1394779/150‑1 (MQ=255)
                                                      aCCGGAGAAATTACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAg                                                                                                 <  2:1096681/150‑1 (MQ=255)
                                                      aCCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAg                                                                                                 >  2:263627/1‑150 (MQ=255)
                                                       ccGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTg                                                                                                         <  2:494133/141‑1 (MQ=255)
                                                       ccGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTg                                                                                                         >  1:494133/1‑141 (MQ=255)
                                                       ccGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGt                                                                                                <  2:319229/150‑1 (MQ=255)
                                                         ggAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGa                                                                                              <  1:77292/150‑1 (MQ=255)
                                                         ggAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGa                                                                                              <  2:151405/150‑1 (MQ=255)
                                                         ggAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGa                                                                                              >  2:1252764/1‑150 (MQ=255)
                                                            gaAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAc                                                                                           >  2:1061200/1‑150 (MQ=255)
                                                             aaaTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAcc                                                                                          <  1:653882/150‑1 (MQ=255)
                                                             aaaTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAcc                                                                                          <  1:1061200/150‑1 (MQ=255)
                                                               aTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTg                                                                                        >  2:292370/1‑150 (MQ=255)
                                                                    cagccacTTACATTTTTCAGGTCAACGAAATAGGATAGTTAAACGCCTTCAGAATGTTACAAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAc                                                                                   >  2:376798/1‑150 (MQ=255)
                                                                     agccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAcc                                                                                  >  2:527282/1‑150 (MQ=255)
                                                                       ccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGt                                                                                <  2:681939/150‑1 (MQ=255)
                                                                          gTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTa                                                                                                                                         <  1:1257421/90‑1 (MQ=255)
                                                                          gTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTa                                                                                                                                         >  2:1257421/1‑90 (MQ=255)
                                                                                    tCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTg                                                                   <  1:1095213/150‑1 (MQ=255)
                                                                                     ccGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGc                                                                  >  2:440464/1‑150 (MQ=255)
                                                                                     ccGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGc                                                                  <  1:263627/150‑1 (MQ=255)
                                                                                      cGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCa                                                                 >  1:430475/1‑150 (MQ=255)
                                                                                       ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAg                                                                >  2:1319176/1‑150 (MQ=255)
                                                                                       ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCCTTAAGTTTGCCCTGCAg                                                                >  2:275930/1‑150 (MQ=255)
                                                                                          cAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTg                                                             <  1:418481/150‑1 (MQ=255)
                                                                                           aaCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                            >  2:344859/1‑150 (MQ=255)
                                                                                            aCGAAATCGGATAGTTAAACTCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTTAg                                                           <  1:860464/150‑1 (MQ=255)
                                                                                            aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGc                                                                  >  2:180996/1‑143 (MQ=255)
                                                                                            aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGc                                                                  <  1:180996/143‑1 (MQ=255)
                                                                                            aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAg                                                           <  2:1539123/150‑1 (MQ=255)
                                                                                                    ggATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTcc                                                   <  1:244935/150‑1 (MQ=255)
                                                                                                      aTAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCGGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACGTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAg                                                 <  2:1128275/150‑1 (MQ=255)
                                                                                                       tAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGc                                                >  1:947154/1‑150 (MQ=255)
                                                                                                           tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAc                                            >  2:636315/1‑150 (MQ=255)
                                                                                                           tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAc                                            >  2:637204/1‑150 (MQ=255)
                                                                                                             aaCGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAccc                                          <  2:474232/150‑1 (MQ=255)
                                                                                                               cGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGa                                        <  1:1288282/150‑1 (MQ=255)
                                                                                                                 ccTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACt                                      >  1:1451683/1‑150 (MQ=255)
                                                                                                                  cTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACtt                                     >  1:330664/1‑150 (MQ=255)
                                                                                                                     cAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                  <  1:625064/150‑1 (MQ=255)
                                                                                                                      aGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgttg                                 <  1:245314/150‑1 (MQ=255)
                                                                                                                       gAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTgg                                <  1:895921/150‑1 (MQ=255)
                                                                                                                            gTACCAATCTCAACGGTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAtt                           <  1:325655/150‑1 (MQ=255)
                                                                                                                            ctACCAATCTCAACGTTTTCCAGCCGGGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCTCCGTTAAGTTTGCCCTCCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAtt                           <  1:319734/149‑1 (MQ=255)
                                                                                                                              aCCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTg                         <  1:1436073/150‑1 (MQ=255)
                                                                                                                              aCCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTg                         <  1:440464/150‑1 (MQ=255)
                                                                                                                                 aaTCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGt                      >  1:1204838/1‑150 (MQ=255)
                                                                                                                                     tcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                  >  1:72417/1‑150 (MQ=255)
                                                                                                                                      cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCaa                 <  1:62559/150‑1 (MQ=255)
                                                                                                                                       aaCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAAc                >  1:807500/1‑150 (MQ=255)
                                                                                                                                        aCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACa               <  1:945021/150‑1 (MQ=255)
                                                                                                                                         cGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACAt              >  2:610366/1‑150 (MQ=255)
                                                                                                                                          aacttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc             >  1:466665/4‑150 (MQ=255)
                                                                                                                                            cttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTg           >  1:1271801/2‑150 (MQ=255)
                                                                                                                                              tCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTa         >  1:1414606/1‑150 (MQ=255)
                                                                                                                                               ccAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTaa        >  1:1623974/1‑150 (MQ=255)
                                                                                                                                                cAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTaaa       <  1:1252764/150‑1 (MQ=255)
                                                                                                                                                     gTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAattca  >  1:398501/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTCA  >  NC_000913/2140251‑2140549

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