Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,302 C→T 19.0% 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,3020CT19.0% 173.4 / 32.4 79R296C (CGT→TGT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (8/7);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

ATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCA  >  NC_000913/3967155‑3967447
                                                                                                                                                   |                                                                                                                                                 
aTCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCg                                                                                                                                                  >  1:1030884/1‑149 (MQ=255)
aTCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCg                                                                                                                                                  >  2:1330959/1‑149 (MQ=255)
aTCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACAACGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGGGGTTTGGATGCCAACGCCATGCATCg                                                                                                                                                  >  1:286824/1‑149 (MQ=255)
  cGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTc                                                                                                                                                <  1:1058066/149‑1 (MQ=255)
   gagaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGTGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTcc                                                                                                                                               >  2:390812/1‑149 (MQ=255)
   gagaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCACCGCCCTGCATCGTcc                                                                                                                                               >  2:1349931/1‑149 (MQ=255)
   gagaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTcc                                                                                                                                               >  1:1030847/1‑149 (MQ=255)
   gagaAGGCCAAACGCCTGGTAGAGCACAAGAAAGACGATATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGAGCGCTGGTGAGGATGCAAACCCCCTGCAACGTcc                                                                                                                                               >  1:1036180/1‑149 (MQ=255)
    agaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCCCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTCACCGGAGGTGTGGATGCCAACGCACAGCATCGAcca                                                                                                                                              >  2:1437281/1‑148 (MQ=255)
       aGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGaaa                                                                                                                                           <  1:585971/149‑1 (MQ=255)
          ccAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGc                                                                                                                                        >  2:244322/1‑149 (MQ=255)
          ccAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGc                                                                                                                                        <  1:382104/149‑1 (MQ=255)
             aaCGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGcttc                                                                                                                                     <  2:45985/149‑1 (MQ=255)
                 ccTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTg                                                                                                                                 >  1:1490303/1‑149 (MQ=255)
                 ccTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTg                                                                                                                                 <  1:312961/149‑1 (MQ=255)
                    ggTTGAGCACAAGAAAGACGTAATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTGGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTg                                                                                                                              <  2:793494/149‑1 (MQ=255)
                          gCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGgcgc                                                                                                                        <  2:1500157/149‑1 (MQ=255)
                              aagaaagaCGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTaa                                                                                                                    >  1:472214/1‑149 (MQ=255)
                                 aaagaCGTTCTCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCGTACACCACCGTTGTTCGGGCGTCAGGGACAGTGTGGACCGGGGGTGTGGAGGCCAACGCACTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGt                                                                                                                 <  1:42365/149‑1 (MQ=255)
                                 aaagaCGTTATCATTCTGCTCGACTCCATGACTCGTCTGGCGCGCGCTTCCACGTCCGTTGTTCCGGCGTCAGGTAAATTGTTGACCGGTGGTGTGGATGGCAACGCCCTGCATCGTCCGAAAGGCTTCTTTGGTGCGGCTCGTAACGt                                                                                                                 <  2:1036180/149‑1 (MQ=255)
                                 aaagaCGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTTCGGCGCGTAACGt                                                                                                                 <  2:1030847/149‑1 (MQ=255)
                                     aCGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGaa                                                                                                             >  2:320691/1‑149 (MQ=255)
                                            cATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGggcg                                                                                                      <  1:1330959/149‑1 (MQ=255)
                                             aTTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTCTTCCGGCGTCAGGTAACGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGggcgg                                                                                                     <  1:242802/149‑1 (MQ=255)
                                             aTTCTGCTCGACTCCATCACTCCTCTCGCGCGCGCTTACCACACCGGTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGggcgg                                                                                                     <  2:87460/149‑1 (MQ=255)
                                                    tCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTg                                                                                              >  1:1371127/1‑149 (MQ=255)
                                                    tCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAggg                                                                                                        >  1:723870/1‑139 (MQ=255)
                                                    tCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAggg                                                                                                        <  2:723870/139‑1 (MQ=255)
                                                     cGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGGTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGa                                                                                             <  2:987519/149‑1 (MQ=255)
                                                      gACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGAc                                                                                            <  1:426677/149‑1 (MQ=255)
                                                      gACTCCATCACTCGTCTGGCGCGCGCGTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGAc                                                                                            <  1:733135/149‑1 (MQ=255)
                                                       acaccATCACTCGTCTGGCGCGAGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGCATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGAcc                                                                                           <  1:139125/146‑1 (MQ=255)
                                                          ccATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCAtt                                                                                        >  1:1279430/1‑149 (MQ=255)
                                                           cATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCGGACCATTa                                                                                       >  2:977638/1‑149 (MQ=255)
                                                            aTCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTAt                                                                                      >  1:895978/1‑149 (MQ=255)
                                                            aTCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGGCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTAt                                                                                      <  2:1030884/149‑1 (MQ=255)
                                                             tCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATc                                                                                     <  1:1442776/149‑1 (MQ=255)
                                                              cACTCGTCTGGTGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATcg                                                                                    <  2:668092/149‑1 (MQ=255)
                                                                cTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTAGAACTCTGGTGTCGAAGCAAACGAAATACAATGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATcgcg                                                                                  >  2:1222197/1‑149 (MQ=255)
                                                                 tCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGa                                                                                 >  1:1227596/1‑149 (MQ=255)
                                                                       ggCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGgcgc                                                                           <  2:745085/149‑1 (MQ=255)
                                                                           cgcgcTTACAACACCGTTGTTCCGGCGTCCGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCGTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTAt                                                                       <  1:566771/149‑1 (MQ=255)
                                                                                    aacaCCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCTGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACAgg                                                              >  2:759370/1‑149 (MQ=255)
                                                                                    aacaCCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCgg                                                              >  1:251652/1‑149 (MQ=255)
                                                                                    aacaCCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCgg                                                              <  2:688378/149‑1 (MQ=255)
                                                                                       acCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTTTCGAGTCCGGtta                                                           >  1:271127/1‑148 (MQ=255)
                                                                                       acCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTc                                                           <  2:187955/149‑1 (MQ=255)
                                                                                       acCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGATGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGAAGGCGCTTATCGATACCGGTTc                                                           >  1:297896/1‑149 (MQ=255)
                                                                                          gttgttCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTaa                                                        >  1:66079/1‑149 (MQ=255)
                                                                                          gttgttCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGGGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATAGATACCGGTTCtta                                                        >  2:300108/1‑147 (MQ=255)
                                                                                            tgttCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTaaaa                                                      <  2:1490303/149‑1 (MQ=255)
                                                                                             gttCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAAt                                                     <  2:1101694/149‑1 (MQ=255)
                                                                                                ccGGCTTCAGTTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGTAAGGCTTCTTTGGTGCGGCGCGTAACGTGGAAGATGGCGGCATCCTGACCATTATCGCGACGGCGCTTTTCGATACCGGTTCTAAAATGGa                                                  <  1:747454/149‑1 (MQ=255)
                                                                                                 cGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCACCGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGAc                                                 >  1:986405/1‑149 (MQ=255)
                                                                                                 cGGCGTCAGGTAAAGTGTTGACAGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACTCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATAGCGACGGCGCTGATCGATACCGGTTCTAAAATGGAc                                                 >  2:416796/1‑149 (MQ=255)
                                                                                                    cGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGaa                                              <  2:760593/149‑1 (MQ=255)
                                                                                                          gTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATc                                        >  2:681890/1‑149 (MQ=255)
                                                                                                          gTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATc                                        >  2:781091/1‑149 (MQ=255)
                                                                                                          gTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATc                                        >  2:1071364/1‑149 (MQ=255)
                                                                                                          gTAAAGTGGTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATc                                        <  1:1331385/149‑1 (MQ=255)
                                                                                                           tAAAGTGTTGACCGGTGGTGTGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGCCCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTa                                      >  1:228249/1‑149 (MQ=255)
                                                                                                              aGTGTTGACCGGTGGTGTGTATGCCAACGCCCTGAATTGTCCGAAACGCTTCTTTGGTGCTGCGCGTAACGTGGAAGAGTGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACg                                    <  1:1499181/149‑1 (MQ=255)
                                                                                                                tgtTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGaa                                  >  2:1155704/1‑149 (MQ=255)
                                                                                                                  ttGACCGGTGGTGTGGATGCCAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAaga                                <  2:66079/149‑1 (MQ=255)
                                                                                                                         gtggtgTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTCCGGAGCGTAACGTGGAAGAGGGCGGCAGCCTGCCCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTaaa                         >  2:470790/1‑149 (MQ=255)
                                                                                                                          tggtgTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAg                        <  2:986405/149‑1 (MQ=255)
                                                                                                                          gggggcggACGCCAACGCCCTGCAGCGTCCGAAACGCTTCTTTGGGGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTCTCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAg                        <  2:1498069/143‑1 (MQ=255)
                                                                                                                            gtgtGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGt                      >  2:26981/1‑149 (MQ=255)
                                                                                                                               tgGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGCCCATTATCGCGACGGCGCTTATCGCTACCGGGTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACa                   >  1:395797/1‑149 (MQ=255)
                                                                                                                               gggATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACa                   <  1:609727/148‑1 (MQ=255)
                                                                                                                                 ttaagcAACGCTCTGCATTGTCCGAAACGATTCTTTGGTGCGGCGCGTTACGTGGAAGAGGGCGGCAGCCTGACCATTTTCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGTAGTTATCTACGAAGAGTTTAAAGGTACAgg                 <  2:143061/144‑1 (MQ=255)
                                                                                                                                    gCCAACGCCCTGCATCGTCCGAAACTCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCaa              <  2:251652/149‑1 (MQ=255)
                                                                                                                                    gCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCTCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCaa              <  1:62511/149‑1 (MQ=255)
                                                                                                                                      cAACGCCCTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACa            >  2:1535583/1‑149 (MQ=255)
                                                                                                                                           cccTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGaa       <  1:1349931/149‑1 (MQ=255)
                                                                                                                                           ccatGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCTACGGCGCTTATCGATACAGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGaa       <  2:668267/146‑1 (MQ=255)
                                                                                                                                             cTGCATTGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACt     >  2:573368/1‑149 (MQ=255)
                                                                                                                                                cATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCa  >  1:457816/1‑149 (MQ=255)
                                                                                                                                                cATCGTCCGAAACGCTTATTTGGTGCTGCGCGTAACGTGGAAGAGGGCGGCAGCCTTACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCa  <  1:1437281/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCA  >  NC_000913/3967155‑3967447

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: