Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,302 C→T R296C (CGT→TGT)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,3020CT99.3% 505.3 / NA 152R296C (CGT→TGT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (73/78);  minor base A (0/1);  total (73/79)
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 152 total aligned reads displayed.

ATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCAC  >  NC_000913/3967155‑3967448
                                                                                                                                                   |                                                                                                                                                  
aTCGATAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACTTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGt                                                                                                                                                  <  2:1464542/150‑1 (MQ=255)
  ggagaAGGCCAAAAGCCTGGTTGAGCACAAGTAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTAGAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGACAACGCCCTGCATTGTcc                                                                                                                                                <  2:1353049/149‑1 (MQ=255)
   gagaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCg                                                                                                                                               >  1:253066/1‑150 (MQ=255)
       aGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTTTGATTCTGCTCGAAGCCATCCCTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCAGTGTCCGTAAc                                                                                                                                           >  2:106840/1‑150 (MQ=255)
          ccAAACGCCTGGTGGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGct                                                                                                                                        <  1:1759551/150‑1 (MQ=255)
            aaaCGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGcttc                                                                                                                                      <  2:1041825/150‑1 (MQ=255)
            aaaCGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGcttc                                                                                                                                      >  2:1936732/1‑150 (MQ=255)
             aaCGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAAAGCCCTGCATTGTCCGAAACGcttct                                                                                                                                     >  1:39001/1‑150 (MQ=255)
              aCGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGcttctt                                                                                                                                    >  2:224381/1‑150 (MQ=255)
               cGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGGCCGAAACGCTTCttt                                                                                                                                   <  1:121392/150‑1 (MQ=255)
                    ggTTGAGCACAAGAAAGACGTTATCATTCCGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGGGGTGTGGATGCCAACGCCCTGCATTGTCCGACACGCTTCTTTGGTgc                                                                                                                              <  1:1181671/150‑1 (MQ=255)
                      ttGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTgcgg                                                                                                                            <  1:814551/150‑1 (MQ=255)
                       tGAGCACAAGAATGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCACTGCATTGTCCGAAACGGTTCTTTGGTgcggc                                                                                                                           <  2:1601587/150‑1 (MQ=255)
                        gAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTGTGGTgcggcg                                                                                                                          >  1:362348/1‑150 (MQ=255)
                        gAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGGTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTgcggcg                                                                                                                          <  2:243726/150‑1 (MQ=255)
                          gCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGcgcg                                                                                                                        >  1:545426/1‑150 (MQ=255)
                          gCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGcgcg                                                                                                                        <  1:864827/150‑1 (MQ=255)
                          gCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGGTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCACCGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGcgcg                                                                                                                        >  2:1028462/1‑150 (MQ=255)
                           cacaAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGt                                                                                                                       >  2:472242/1‑150 (MQ=255)
                             caAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTaa                                                                                                                     >  2:1005099/1‑150 (MQ=255)
                               agaaagaCGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACg                                                                                                                   >  1:1923096/1‑150 (MQ=255)
                                gaaagaCGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGt                                                                                                                  >  2:1589889/1‑150 (MQ=255)
                                     aCGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCACTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAag                                                                                                             <  1:681763/150‑1 (MQ=255)
                                      cGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCACTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAaga                                                                                                            >  1:453597/1‑150 (MQ=255)
                                        ttATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAgg                                                                                                          >  1:1685183/1‑150 (MQ=255)
                                            cATTATGCGCGACTCTATAACTCGGCTGGCGCGCGCGTACAACACCGTTGTCCCGGCGTCATGTAATTTGTTGACCGTTTGTGTGGATGCCAACGCCCTGGATTGTGCGTAACGCTTGTTTGGTTCGGCGCGTAACGTGGAAGAGggcgg                                                                                                      <  2:121225/150‑1 (MQ=255)
                                             aTTCAGCCCGACTCCATCACTCGTCTGGCGCGCGCGTACACCACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTTGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGggcggc                                                                                                     <  2:674577/150‑1 (MQ=255)
                                               tCTGCTCGCCTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGGCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAg                                                                                                   <  2:1302818/150‑1 (MQ=255)
                                               tCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAg                                                                                                   <  2:273584/150‑1 (MQ=255)
                                               tCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTAGTCCGGTGGTGTGGATGCCGACGCCCTGCATTGTCCGAAACGATTCTTTGGTGCGGCGCGTATCGTGGAAGAGGGCGGCAg                                                                                                   >  1:276589/1‑150 (MQ=255)
                                                 tGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGcc                                                                                                 <  1:75371/150‑1 (MQ=255)
                                                 tGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGcc                                                                                                 >  2:872652/1‑150 (MQ=255)
                                                   cTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTg                                                                                               >  1:55248/1‑150 (MQ=255)
                                                   cTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTg                                                                                               >  2:1192358/1‑150 (MQ=255)
                                                   cTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGAAACGTGGAAGAGGGCGGCAGCCTg                                                                                               >  1:251366/1‑150 (MQ=255)
                                                       aCTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCa                                                                                           <  2:27514/150‑1 (MQ=255)
                                                       aCTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCa                                                                                           <  2:396440/150‑1 (MQ=255)
                                                        cTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAACCGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGCCCAt                                                                                          >  1:59068/1‑150 (MQ=255)
                                                           cATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTAt                                                                                       >  1:825460/1‑150 (MQ=255)
                                                              cACTCGTCTGGCGCGCGCTTACACCACCGTTGTTCCGGCGTCAGGTAAAGTGGTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCGTCTTGGGTGCGGCGCGTAACGTGGATGAGGGCGGCAGCCTGACCATTATcgc                                                                                    >  1:382057/1‑150 (MQ=255)
                                                                  cGTCTGGCGCGCGCTTACAACACCGTTGTTACGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACg                                                                                >  2:419492/1‑150 (MQ=255)
                                                                   gTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACgg                                                                               >  2:1833789/1‑150 (MQ=255)
                                                                     cTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGCCGgcg                                                                             >  1:16271/1‑150 (MQ=255)
                                                                     cTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGgcg                                                                             <  1:1961172/150‑1 (MQ=255)
                                                                      tGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGCCGgcgg                                                                            >  1:416417/1‑149 (MQ=255)
                                                                       ggCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCt                                                                           <  2:1360778/150‑1 (MQ=255)
                                                                       ggCGCGCGCTTACAACACCGATGTTCCGGCGTCAGGTAAAGGGTTGACCGGTGGGGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCt                                                                           <  1:1536366/150‑1 (MQ=255)
                                                                        gcgcgtgcTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCAAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACAATTATCGCGACGGCGCtt                                                                          >  2:364339/1‑150 (MQ=255)
                                                                         cgcgcgcTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTa                                                                         <  2:716991/150‑1 (MQ=255)
                                                                            gcgcTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCg                                                                      <  2:1836442/150‑1 (MQ=255)
                                                                            gcgcTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCGGACCATTATCGCGACGGCGCTTATCg                                                                      <  2:124385/150‑1 (MQ=255)
                                                                                 tACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCCAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGCCAGCCTGACCATTATCGCGACGGCGCTTATCGATAcc                                                                 >  2:1173251/1‑150 (MQ=255)
                                                                                 aacaacaCCGTTGTTCCGGCGCCAGGTAAAGTGTTGACCGGTGGAGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATAcc                                                                 <  1:676576/149‑1 (MQ=255)
                                                                                      cacCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTc                                                            >  1:722462/1‑150 (MQ=255)
                                                                                       acCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCt                                                           >  2:1430477/1‑150 (MQ=255)
                                                                                        ccGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTa                                                          >  1:1311022/1‑150 (MQ=255)
                                                                                         cGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTaa                                                         <  2:870809/150‑1 (MQ=255)
                                                                                              ttCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATgg                                                    <  2:262248/150‑1 (MQ=255)
                                                                                               tCCGGCGGCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGa                                                   <  2:589545/150‑1 (MQ=255)
                                                                                                ccGGCGTCAGGTAAAGTGTTTAACTGTGGTGTGGAAGCCAAAGCCCTACATTGACCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGA‑CATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACg                                                 >  2:223285/1‑150 (MQ=255)
                                                                                                ccGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGTGCTTATCGATACCGGTTCTAAAATGGAc                                                  <  2:545426/150‑1 (MQ=255)
                                                                                                 cGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACg                                                 <  1:1820453/150‑1 (MQ=255)
                                                                                                  ggCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGa                                                >  1:159569/1‑150 (MQ=255)
                                                                                                    cGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAg                                              >  2:1147318/1‑150 (MQ=255)
                                                                                                     gTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGt                                             <  1:137739/150‑1 (MQ=255)
                                                                                                     gTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGt                                             <  2:1678763/150‑1 (MQ=255)
                                                                                                           tAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACTCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTa                                       <  1:1443351/150‑1 (MQ=255)
                                                                                                            aaaGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTAc                                      <  2:55248/150‑1 (MQ=255)
                                                                                                                tgtTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAag                                  <  1:643362/150‑1 (MQ=255)
                                                                                                                tgtTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAag                                  <  2:825460/150‑1 (MQ=255)
                                                                                                                  ttGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAagag                                >  1:1163516/1‑150 (MQ=255)
                                                                                                                   tGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGt                               >  1:1514496/1‑150 (MQ=255)
                                                                                                                    gACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGtt                              >  2:232692/1‑150 (MQ=255)
                                                                                                                         gtggtgTGGATGCCAACGCCCTGCATTGTCCGAAACGCTGCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAg                         <  1:858066/150‑1 (MQ=255)
                                                                                                                         gtggtgTGGATGCCAACGCCCTGCATTCTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAg                         <  1:1147357/150‑1 (MQ=255)
                                                                                                                          tggtgTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAgg                        <  1:1192358/150‑1 (MQ=255)
                                                                                                                           ggtgTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGt                       >  1:768013/1‑150 (MQ=255)
                                                                                                                              gtgGATGCCAACGCCCTGCGTTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACa                    <  2:775788/150‑1 (MQ=255)
                                                                                                                              gtgGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACa                    >  1:1516490/1‑150 (MQ=255)
                                                                                                                               tgGATGCCAACGCCCTGCATTGTCCGAACAGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTTTCGCGACGGCGCTTATAGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAg                   <  2:261100/150‑1 (MQ=255)
                                                                                                                                 gATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGc                 <  1:1807454/150‑1 (MQ=255)
                                                                                                                                 gATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGc                 <  2:1098681/150‑1 (MQ=255)
                                                                                                                                 gATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGc                 >  1:798984/1‑150 (MQ=255)
                                                                                                                                  aTGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCa                <  1:1079808/150‑1 (MQ=255)
                                                                                                                                    gCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAAc              <  1:1005099/150‑1 (MQ=255)
                                                                                                                                     ccAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGATGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACTAAGAGTTTAAAGCTACAGGCAACa             >  1:1227609/1‑150 (MQ=255)
                                                                                                                                     ccAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGCCCATTATCGCGCCGGCGCTTATCGATACCGGTTCTAAAAGGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACa             >  1:1690239/1‑150 (MQ=255)
                                                                                                                                     ccAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACa             >  1:749009/1‑150 (MQ=255)
                                                                                                                                     ccAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACa             >  1:1722229/1‑150 (MQ=255)
                                                                                                                                     ccAAAGCCCTGCATTGTCCGAACCGCGTCTTTGGTGCGGCGCGGAACGGGGAAGCGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGATGTTATCTACGAAGAGTTTAAAGGTACAGGCAACa             <  2:1035131/150‑1 (MQ=255)
                                                                                                                                      cAACGCCCAGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTTTCGCGACGGCGCGTTTCGATACCGGTTATAAAATGGACGCAGTTTTCGTCGAAGAGTTTAAAGGTACAGGCAACAt            >  2:657451/1‑150 (MQ=255)
                                                                                                                                        acccccTGCATTGTCCGAAACGCTACTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATAACGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATgg          >  2:984742/4‑150 (MQ=255)
                                                                                                                                         cGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGa         >  2:1512926/1‑150 (MQ=255)
                                                                                                                                          gCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGaa        >  2:740119/1‑150 (MQ=255)
                                                                                                                                           gccTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAAc       <  2:53258/149‑1 (MQ=255)
                                                                                                                                            ccTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACt      >  1:128848/1‑150 (MQ=255)
                                                                                                                                            ccTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAAAGTGGCAGAGGGCGGCAGCATGACCATTATCGCTACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGATATCTTCTAAGAGTTTAAAGGTACATGCAAAATGGAACt      >  1:326129/1‑150 (MQ=255)
                                                                                                                                              tGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGc    <  1:1941278/150‑1 (MQ=255)
                                                                                                                                               gCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCa   <  1:1613991/150‑1 (MQ=255)
                                                                                                                                                cATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCAc  <  1:271338/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCAC  >  NC_000913/3967155‑3967448

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