Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,310,602 A→G 5.3% V189V (GTT→GTC pnp ← polynucleotide phosphorylase/polyadenylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,310,6020AG 5.3% 227.0 / 2.1 76V189V (GTT→GTCpnppolynucleotide phosphorylase/polyadenylase
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base G (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CGGCTCCGGCTGCCAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACA  >  NC_000913/3310458‑3310750
                                                                                                                                                |                                                                                                                                                    
cggctccggctGCCAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCa                                                                                                                                                  <  1:1599399/149‑1 (MQ=255)
  gctccggctGCCAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATc                                                                                                                                                >  2:690177/1‑149 (MQ=255)
   ctccggctGCCAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCACAGCAGTTTCGCTTCAGATTCAACCATCa                                                                                                                                               >  2:286706/1‑149 (MQ=255)
          tGCCAGTCCCAACGCGGTTTACCGGCTTCTTTCTCCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGc                                                                                                                                        <  1:697601/149‑1 (MQ=255)
            ccAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCgg                                                                                                                                      >  2:484372/1‑149 (MQ=255)
            ccAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCgg                                                                                                                                      <  2:238282/149‑1 (MQ=255)
             ctgTCCCAAAGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACAGAACACTACTGCGCCCAGCATCTGGGCTACGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGc                                                                                                                                     <  2:171104/147‑1 (MQ=255)
              aGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCGACCATCAGTACAGCGGCt                                                                                                                                    >  1:1036947/1‑149 (MQ=255)
              aGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCGACCATCAGTACAGCGGCt                                                                                                                                    >  1:1528931/1‑149 (MQ=255)
                tCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTc                                                                                                                                  <  1:1193674/149‑1 (MQ=255)
                 cccAACGCGGTTTACCGGCTTCTTTCTCCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCa                                                                                                                                 <  2:799629/149‑1 (MQ=255)
                   cAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACCCTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGt                                                                                                                               >  2:237884/1‑149 (MQ=255)
                   cAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGt                                                                                                                               <  2:1145105/149‑1 (MQ=255)
                    aaCGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTa                                                                                                                              <  2:84134/149‑1 (MQ=255)
                      cgcgGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTAcc                                                                                                                            <  2:837761/149‑1 (MQ=255)
                           tttACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGcaa                                                                                                                       <  2:949581/149‑1 (MQ=255)
                           tttACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGcaa                                                                                                                       <  2:590098/149‑1 (MQ=255)
                              aCCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGcaacaa                                                                                                                    <  2:1109929/149‑1 (MQ=255)
                               ccGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCATCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCaacaac                                                                                                                   <  2:296759/149‑1 (MQ=255)
                                cGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAAcc                                                                                                                  >  2:430717/1‑149 (MQ=255)
                                cGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAAcc                                                                                                                  >  1:1218818/1‑149 (MQ=255)
                                 ggCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCa                                                                                                                 >  1:428499/1‑149 (MQ=255)
                                    ttcttTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGAt                                                                                                              <  1:1264677/149‑1 (MQ=255)
                                       tttCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCa                                                                                                           >  2:1255240/1‑149 (MQ=255)
                                       tttCACCAGTTCATTGATGTTCTGAATAACAACCTCATGTTGTTAATGACCGAGCACTACTGCGCCCAGCATCAGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCACTACACCGGCTTCAGTCACGGCAACAACCAGATGCa                                                                                                           >  1:796052/1‑149 (MQ=255)
                                           aCCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGttt                                                                                                       <  1:1349665/149‑1 (MQ=255)
                                                  cATTGATGTTCTGTATAACTACCTGATGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGctctct                                                                                                <  2:1008042/149‑1 (MQ=255)
                                                   aTTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCtt                                                                                               <  1:237884/149‑1 (MQ=255)
                                                   aTTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGCCAACAACCAGATCCAGTTTGCTCTCtt                                                                                               <  2:1159938/149‑1 (MQ=255)
                                                   aTTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAACAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCtt                                                                                               >  2:750574/1‑149 (MQ=255)
                                                       aTGTTCTGAATAACAACCTCCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCGCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCGCTCTTTCAg                                                                                           <  1:268850/149‑1 (MQ=255)
                                                        tGTTCTGAATAAGAACCTGCTGTTGTTCATGAGCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACTGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCTGc                                                                                          <  2:612303/149‑1 (MQ=255)
                                                          ttCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTc                                                                                        >  1:1040623/1‑149 (MQ=255)
                                                          ttCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTc                                                                                        >  1:10101/1‑149 (MQ=255)
                                                             tGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                                                     >  1:344120/1‑149 (MQ=255)
                                                              gAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTcc                                                                                    >  1:673548/1‑149 (MQ=255)
                                                                 tAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGa                                                                                 >  2:1486478/1‑149 (MQ=255)
                                                                 tAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGa                                                                                 <  1:855736/149‑1 (MQ=255)
                                                                   acaaaCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGt                                                                               <  1:228364/149‑1 (MQ=255)
                                                                    caacCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCGACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTc                                                                              <  2:1528931/149‑1 (MQ=255)
                                                                    caacCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCGACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTc                                                                              <  2:1036947/149‑1 (MQ=255)
                                                                     aacCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCg                                                                             >  1:695926/1‑149 (MQ=255)
                                                                      acCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCgg                                                                            <  2:59916/149‑1 (MQ=255)
                                                                      acCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCgg                                                                            <  1:656036/149‑1 (MQ=255)
                                                                       ccTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCggg                                                                           >  1:651722/1‑149 (MQ=255)
                                                                       cagcaggttgttCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCggg                                                                           >  2:309401/7‑149 (MQ=255)
                                                                           ctgTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCa                                                                       <  2:1040623/149‑1 (MQ=255)
                                                                           ctgTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCa                                                                       <  1:1509758/149‑1 (MQ=255)
                                                                           ctgTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCa                                                                       <  2:135148/149‑1 (MQ=255)
                                                                            tgttgtTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAg                                                                      <  1:430717/149‑1 (MQ=255)
                                                                              ttgttCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAgta                                                                    >  2:1327860/1‑149 (MQ=255)
                                                                                    aTGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACt                                                              <  1:286706/149‑1 (MQ=255)
                                                                                           aaCACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCAtt                                                       >  2:615806/1‑149 (MQ=255)
                                                                                            acacTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTg                                                      <  2:458318/149‑1 (MQ=255)
                                                                                            acacTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCAGCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTg                                                      <  1:966344/149‑1 (MQ=255)
                                                                                             cacTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGGTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGGACTGGTCATTGa                                                     >  2:187951/1‑149 (MQ=255)
                                                                                                tactGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGt                                                  >  2:1111718/1‑149 (MQ=255)
                                                                                                tactGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCGGAGTCGGGTTCAGTACGTACTGGTCATTGATGt                                                  >  1:1399991/1‑149 (MQ=255)
                                                                                                   tGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAAc                                               >  2:1142382/1‑149 (MQ=255)
                                                                                                   tGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAAc                                               >  1:311029/1‑149 (MQ=255)
                                                                                                        ccAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTAcg                                          <  1:484372/149‑1 (MQ=255)
                                                                                                        ccAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTAcg                                          >  2:360193/1‑149 (MQ=255)
                                                                                                         cAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTAcgc                                         <  2:651722/149‑1 (MQ=255)
                                                                                                         cAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTAAGTACGTACTGGTCATTGATGTAACCTAcgc                                         >  1:1018573/1‑149 (MQ=255)
                                                                                                                     ttCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTGCCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCaa                             >  2:764974/1‑149 (MQ=255)
                                                                                                                     ttCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCaa                             >  1:1087332/1‑149 (MQ=255)
                                                                                                                          tCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCggg                        >  2:1180464/1‑149 (MQ=255)
                                                                                                                               agTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCAtt                   >  2:158155/1‑149 (MQ=255)
                                                                                                                                      cTTCAGATTCCACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCACCTCGTCCTGAGTCGGCTTCAGTACGTACTGGTCATTGATGTAACCAACGCGGGACGCACCAATAGGGCCATTGAACGGa            >  2:184945/1‑149 (MQ=255)
                                                                                                                                      cTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGa            >  2:1013330/1‑149 (MQ=255)
                                                                                                                                       ttCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCGATCGGGCCATTGAACGGaa           <  1:1327860/149‑1 (MQ=255)
                                                                                                                                       ttCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGaa           <  1:690177/149‑1 (MQ=255)
                                                                                                                                          aGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATAc        <  1:309401/149‑1 (MQ=255)
                                                                                                                                            aTTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACca      >  2:305477/1‑149 (MQ=255)
                                                                                                                                            aTTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACca      <  2:1218818/149‑1 (MQ=255)
                                                                                                                                            aTTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCAGTGAACGGAATACca      >  1:13957/1‑149 (MQ=255)
                                                                                                                                                aaCCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACcagaca  >  1:403171/1‑149 (MQ=255)
                                                                                                                                                aaCCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACcagaca  >  1:798465/1‑149 (MQ=255)
                                                                                                                                                aaCCATCAGTACAGCGGCTTCAGTACAGGCAACAACCAGCTCCAGTATGCTCTCTTTCAGCTCGTCGTGAGTCGGGTTCAGTACGGACTGGTCATTGATGTCACCTACGCGGGCAGCACCCAGCGGGCCAGTCCGCGGCATACcagaca  >  1:251580/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGGCTCCGGCTGCCAGTCCCAACGCGGTTTACCGGCTTCTTTCACCAGTTCATTGATGTTCTGAATAACAACCTGCTGTTGTTCATGACCGAACACTACTGCGCCCAGCATCTGGTCTTCGCTCAGCAGTTGAGCTTCAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACA  >  NC_000913/3310458‑3310750

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