Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,453,100 G→T 7.7% V57V (GTC→GTA rpsJ ← 30S ribosomal subunit protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,453,1000GT 7.7% 264.4 / 5.3 77V57V (GTC→GTArpsJ30S ribosomal subunit protein S10
Reads supporting (aligned to +/- strand):  ref base G (44/27);  new base T (3/3);  total (47/30)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGCTGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATC  >  NC_000913/3452967‑3453231
                                                                                                                                     |                                                                                                                                   
ggCTGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGc                                                                                                                               >  1:1028967/1‑140 (MQ=255)
 gCTGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACTGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCg                                                                                                                              <  2:955935/140‑1 (MQ=255)
 gCTGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACTGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCg                                                                                                                              <  2:955856/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  1:1145890/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  2:923359/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  2:910634/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  2:434682/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  1:1483076/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  2:1589517/140‑1 (MQ=255)
   tGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCggg                                                                                                                            <  1:486244/140‑1 (MQ=255)
     aTCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGga                                                                                                                          >  2:516305/1‑140 (MQ=255)
     aTCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGga                                                                                                                          >  2:1250045/1‑140 (MQ=255)
     aTCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGga                                                                                                                          >  2:1250044/1‑140 (MQ=255)
        tGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCACCGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGAt                                                                                                                       <  1:981944/140‑1 (MQ=255)
        tGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGTGAGAt                                                                                                                       <  2:1131385/140‑1 (MQ=255)
        tGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGTGAGAt                                                                                                                       <  2:1131280/140‑1 (MQ=255)
        tGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGAt                                                                                                                       <  1:146059/140‑1 (MQ=255)
        tGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGAt                                                                                                                       <  1:146027/140‑1 (MQ=255)
          cACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTTATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCa                                                                                                                     >  1:925772/1‑140 (MQ=255)
          cACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCa                                                                                                                     >  2:163508/1‑140 (MQ=255)
               cTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACa                                                                                                                >  2:1381755/1‑140 (MQ=255)
               cTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACa                                                                                                                >  2:1381739/1‑140 (MQ=255)
               cTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACa                                                                                                                >  2:1381738/1‑140 (MQ=255)
               cTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACa                                                                                                                >  2:641176/1‑140 (MQ=255)
                     cGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAg                                                                                                          <  2:1596214/140‑1 (MQ=255)
                       gCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTCCGGAGTTCGTACTGATCGCGCGCGTCTGGGTTGACGTGCGGGGAGAGCAGACCAGTGAcgcg                                                                                                        >  2:1689608/1‑140 (MQ=255)
                       gCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAgcg                                                                                                        >  2:955870/1‑140 (MQ=255)
                       gCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAgcg                                                                                                        >  2:665960/1‑140 (MQ=255)
                       gCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAgcg                                                                                                        >  1:1364641/1‑140 (MQ=255)
                       gCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAgcg                                                                                                        >  2:1123831/1‑140 (MQ=255)
                       gCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCACAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTGTGTTGACGTGCGGGGAGATCAGAACAGTGAAgcg                                                                                                        >  1:1364633/1‑140 (MQ=255)
                          gCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGctc                                                                                                     >  2:44506/1‑140 (MQ=255)
                                     aCGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCg                                                                                          >  1:68381/1‑140 (MQ=255)
                                     aCGCATCAGAGCATCAACGGTTTTCTCGGGTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCg                                                                                          >  1:68448/1‑140 (MQ=255)
                                         aTCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAg                                                                                      >  2:253761/1‑140 (MQ=255)
                                         aTCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTGTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAg                                                                                      >  2:462875/1‑140 (MQ=255)
                                               gCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGAt                                                                                <  2:132377/140‑1 (MQ=255)
                                                cATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTTACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATc                                                                               >  1:871309/1‑140 (MQ=255)
                                                cATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTTACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATc                                                                               >  1:872199/1‑140 (MQ=255)
                                                cATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTTACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATc                                                                               >  1:871993/1‑140 (MQ=255)
                                                   cAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGa                                                                            >  1:1652864/1‑140 (MQ=255)
                                                   cAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGa                                                                            >  2:1048103/1‑140 (MQ=255)
                                                         ttttCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACgg                                                                      <  1:163508/140‑1 (MQ=255)
                                                          tttCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGa                                                                     >  1:820947/1‑140 (MQ=255)
                                                           ttCTCTGTTGGATCAAATATGTCAAACAGAAGAAAGTTAGTACGGATTTCGTACTGATCGAGCGAGTCTTTGTTGACGTGAGGGGAGATCAGAACAGTGAATCTCTCTTTTCGTGTCGGCAGCGGGATCGGACCACGGAc                                                                    <  2:185057/140‑1 (MQ=255)
                                                           ttCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGAc                                                                    <  2:185009/140‑1 (MQ=255)
                                                           ttCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGAc                                                                    <  1:1048103/140‑1 (MQ=255)
                                                                     gCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTTACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCa                                                          <  2:871993/140‑1 (MQ=255)
                                                                     gCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTTACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCa                                                          <  2:871309/140‑1 (MQ=255)
                                                                     gCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTGTACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCa                                                          <  2:872199/140‑1 (MQ=255)
                                                                       tCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGt                                                        <  2:299896/140‑1 (MQ=255)
                                                                                tCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                               >  1:1496003/1‑140 (MQ=255)
                                                                                tCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                               >  1:537545/1‑140 (MQ=255)
                                                                                tCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                               >  1:709384/1‑140 (MQ=255)
                                                                                tCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                               >  1:1561735/1‑140 (MQ=255)
                                                                                tCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                               >  1:1561822/1‑140 (MQ=255)
                                                                                tCAACCAGACGCAAGTGAGTACGGATTTAGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                               >  2:783854/1‑140 (MQ=255)
                                                                                            aaGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGAt                                   >  1:693219/1‑140 (MQ=255)
                                                                                            aaGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGAt                                   >  1:693214/1‑140 (MQ=255)
                                                                                               tGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTc                                >  2:575811/1‑140 (MQ=255)
                                                                                                    aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCgg                           >  1:891483/1‑140 (MQ=255)
                                                                                                    aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCgg                           >  1:888313/1‑140 (MQ=255)
                                                                                                    aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCgg                           >  1:888186/1‑140 (MQ=255)
                                                                                                    aCGGATGTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCgg                           <  2:1438008/140‑1 (MQ=255)
                                                                                                        aTTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                       >  1:817156/1‑140 (MQ=255)
                                                                                                             gTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTAGCTTgat                  >  2:1051145/1‑140 (MQ=255)
                                                                                                                 tGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgatcgat              <  2:252306/140‑1 (MQ=255)
                                                                                                                 tGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgatcgat              <  1:138608/140‑1 (MQ=255)
                                                                                                                     cgcgcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGa          <  2:679435/140‑1 (MQ=255)
                                                                                                                      gcgcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGAc         >  2:650475/1‑140 (MQ=255)
                                                                                                                        gcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACga       <  1:768166/140‑1 (MQ=255)
                                                                                                                        gcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACga       >  2:188389/1‑140 (MQ=255)
                                                                                                                             tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCa                                              <  2:854477/96‑1 (MQ=255)
                                                                                                                             tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCa                                              >  1:854477/1‑96 (MQ=255)
                                                                                                                             tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATc  >  2:1589534/1‑140 (MQ=255)
                                                                                                                             tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATc  >  2:411874/1‑140 (MQ=255)
                                                                                                                             tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATc  >  1:480580/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                   
GGCTGATCTGCACGTCTACACCGGCAGCCAGATCCAGACGCATCAGAGCATCAACGGTTTTCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATC  >  NC_000913/3452967‑3453231

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