Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,407,871 T→A 5.6% H406Q (CAT→CAA)  rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,407,8710TA 5.6% 224.8 / 3.0 72H406Q (CAT→CAA) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base A (2/2);  total (34/38)
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.21e-01

GACGGTTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGC  >  NC_000913/4407724‑4408017
                                                                                                                                                   |                                                                                                                                                  
gACGGTTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCAtc                                                                                                                                                   >  2:428359/1‑149 (MQ=255)
  cGGTTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCAtctc                                                                                                                                                 <  2:810698/149‑1 (MQ=255)
  cGGTTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGAGTTTATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCAtctc                                                                                                                                                 <  2:16079/149‑1 (MQ=255)
     tttCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGaa                                                                                                                                              >  1:397680/1‑149 (MQ=255)
     tttCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGaa                                                                                                                                              <  1:373738/149‑1 (MQ=255)
          tcgAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCAGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAagagat                                                                                                                                         >  1:693351/1‑147 (MQ=255)
            gAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGc                                                                                                                                       >  2:1041106/1‑149 (MQ=255)
            gAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGc                                                                                                                                       >  1:339031/1‑149 (MQ=255)
                ggCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATaa                                                                                                                                   <  1:509968/149‑1 (MQ=255)
                 gccgccTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAAc                                                                                                                                  >  2:198205/1‑149 (MQ=255)
                    gccTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTATGCAACTCTTCGAGATGCTTAAccag                                                                                                                               <  1:1098203/149‑3 (MQ=255)
                          cGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATaaa                                                                                                                         <  1:1067681/149‑1 (MQ=255)
                             gCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTg                                                                                                                      <  2:1590035/149‑1 (MQ=255)
                               tACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCt                                                                                                                    <  1:589622/149‑1 (MQ=255)
                               tACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCt                                                                                                                    <  1:65545/149‑1 (MQ=255)
                                aCAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTg                                                                                                                   <  2:144023/149‑1 (MQ=255)
                                 cAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTgg                                                                                                                  >  2:458669/1‑149 (MQ=255)
                                 cAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTgg                                                                                                                  >  2:1520339/1‑149 (MQ=255)
                                 cAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTACATAACCTCTATAAAGTGCTgg                                                                                                                  >  1:472210/1‑149 (MQ=255)
                                    aTTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTgggca                                                                                                               >  1:277246/1‑146 (MQ=255)
                                    aTTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATa                                                                                                               >  1:571051/1‑149 (MQ=255)
                                       cTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAg                                                                                                            <  2:397680/149‑1 (MQ=255)
                                       cTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAg                                                                                                            >  1:469183/1‑149 (MQ=255)
                                        tACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACATCTATAAAGTGCTGGATAAAGc                                                                                                           >  2:1237579/1‑149 (MQ=255)
                                          cGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGccc                                                                                                         <  1:772865/149‑1 (MQ=255)
                                            aaGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGt                                                                                                       <  1:656802/149‑1 (MQ=255)
                                                  tgatgaGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTgaagaa                                                                                                 <  2:469183/149‑1 (MQ=255)
                                                      gaGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAAcgcg                                                                                             <  1:1564220/149‑1 (MQ=255)
                                                       aGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCgg                                                                                            >  2:1219337/1‑149 (MQ=255)
                                                            tcACGCGCGTCTGACCTACACCAAAGTCTGGAATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTCGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATTAAGCCCGTGAAGATCGCGGTGCGa                                                                                       >  2:98360/1‑149 (MQ=255)
                                                             cACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAt                                                                                      <  1:198205/149‑1 (MQ=255)
                                                             cACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAt                                                                                      <  2:563666/149‑1 (MQ=255)
                                                              aCGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtc                                                                                     <  2:527913/149‑1 (MQ=255)
                                                              aCGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtc                                                                                     >  2:920151/1‑149 (MQ=255)
                                                                tcgcgTCTGACCTAAACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtctc                                                                                   <  2:1105021/148‑1 (MQ=255)
                                                                gcgcgTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtctc                                                                                   >  2:435777/1‑149 (MQ=255)
                                                                gcgcgTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtctc                                                                                   <  1:1369481/149‑1 (MQ=255)
                                                                  gcgTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCAt                                                                                 <  2:1105298/149‑1 (MQ=255)
                                                                   ggTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCAtt                                                                                <  2:725265/148‑1 (MQ=255)
                                                                     tCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTg                                                                              >  1:476317/1‑149 (MQ=255)
                                                                     gcTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTg                                                                              <  2:1476296/148‑1 (MQ=255)
                                                                      cTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTga                                                                             <  1:1041106/149‑1 (MQ=255)
                                                                      cTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTga                                                                             <  1:1524063/149‑1 (MQ=255)
                                                                          ccTACACCAAAGTCTGGCATATTCTTCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGc                                                                         <  1:1088673/149‑1 (MQ=255)
                                                                              cacCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCgaag                                                                     <  2:514081/149‑1 (MQ=255)
                                                                               acCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTTAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGATCGCGGGGGGATCTCATTTGAGAGCgaaga                                                                    >  1:1018876/1‑149 (MQ=255)
                                                                               acCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCgaaga                                                                    >  2:638278/1‑149 (MQ=255)
                                                                                 cAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                                                  >  2:1517878/1‑149 (MQ=255)
                                                                                   aaGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCg                                                                <  1:1520339/149‑1 (MQ=255)
                                                                                        tGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGtt                                                           >  1:639076/1‑149 (MQ=255)
                                                                                          gCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                                                         >  1:824378/1‑149 (MQ=255)
                                                                                            atatTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtt                                                       <  2:571051/149‑1 (MQ=255)
                                                                                                 cTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCaa                                                  <  2:339031/149‑1 (MQ=255)
                                                                                                  tGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAAc                                                 >  2:917564/1‑149 (MQ=255)
                                                                                                  tGCAGGGCGATCAGGAGCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAAc                                                 <  2:1300325/149‑1 (MQ=255)
                                                                                                           aTCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGc                                        >  2:657538/1‑149 (MQ=255)
                                                                                                           aTCAGGATCCGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGc                                        >  2:699647/1‑149 (MQ=255)
                                                                                                             cAGGATCTGCGCGCGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGGGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCTTTTTCATCGCTGAACGCCg                                      >  1:430353/1‑149 (MQ=255)
                                                                                                             cAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg                                      <  1:151232/149‑1 (MQ=255)
                                                                                                                gATCTGCGCGAGCAGTACGCCCCGCTGGTTAACCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGCAAGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTACATTTTCAACGCTGAACGCCgaaa                                   >  2:284117/1‑146 (MQ=255)
                                                                                                                 aTCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTAtt                                  <  1:489259/149‑1 (MQ=255)
                                                                                                                     gcgcGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAAc                              <  1:920151/149‑1 (MQ=255)
                                                                                                                             aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGaa                      >  2:451633/1‑149 (MQ=255)
                                                                                                                             aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGaa                      <  1:1083408/149‑1 (MQ=255)
                                                                                                                                aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAg                   <  1:659543/149‑1 (MQ=255)
                                                                                                                                aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAg                   <  1:699647/149‑1 (MQ=255)
                                                                                                                                  gCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGAc                 <  1:1242868/149‑1 (MQ=255)
                                                                                                                                    cccGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGAccc               >  1:1595863/1‑149 (MQ=255)
                                                                                                                                      tagaggTTATGCAACTCTTCGAGATGCTTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAg             >  2:1379716/5‑149 (MQ=255)
                                                                                                                                       gCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGc            <  2:824378/149‑1 (MQ=255)
                                                                                                                                       gCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGACGTATTGAACGTATCGAACAGACCCAGc            >  1:1272224/1‑149 (MQ=255)
                                                                                                                                        gaggTTATGCAACTCTTCGAGATGCTTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                                               >  1:1363529/3‑97 (MQ=25)
                                                                                                                                        gaggTTATGCAACTCTTCGAGATGCTTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                                               <  2:1363529/95‑1 (MQ=25)
                                                                                                                                               aGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATATCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAacgac    >  1:313608/1‑149 (MQ=255)
                                                                                                                                                 cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGAcgc  >  1:375619/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GACGGTTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGC  >  NC_000913/4407724‑4408017

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