Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,407,877 A→C 5.5% E408D (GAA→GAC) ‡ rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,407,8770AC 5.5% 225.0 / 3.6 73E408D (GAA→GAC) ‡rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base C (2/2);  total (35/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.53e-01

TTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAAT  >  NC_000913/4407729‑4408025
                                                                                                                                                    |                                                                                                                                                    
tttCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGaa                                                                                                                                                      <  1:373738/149‑1 (MQ=255)
tttCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGaa                                                                                                                                                      >  1:397680/1‑149 (MQ=255)
     tcgAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCAGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAagagat                                                                                                                                                 >  1:693351/1‑147 (MQ=255)
       gAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGc                                                                                                                                               >  1:339031/1‑149 (MQ=255)
       gAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGc                                                                                                                                               >  2:1041106/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                                                                                              <  2:563666/149‑1 (MQ=255)
                                                        cACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAt                                                                                              <  1:198205/149‑1 (MQ=255)
                                                         aCGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtc                                                                                             >  2:920151/1‑149 (MQ=255)
                                                         aCGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtc                                                                                             <  2:527913/149‑1 (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                              <  1:1083408/149‑1 (MQ=255)
                                                                                                                        aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGaa                              >  2:451633/1‑149 (MQ=255)
                                                                                                                           aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAg                           <  1:699647/149‑1 (MQ=255)
                                                                                                                           aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAg                           <  1:659543/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                                                                       <  2:1363529/95‑1 (MQ=25)
                                                                                                                                   gaggTTATGCAACTCTTCGAGATGCTTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                                                       >  1:1363529/3‑97 (MQ=25)
                                                                                                                                          aGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATATCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAacgac            >  1:313608/1‑149 (MQ=255)
                                                                                                                                            cgagtcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGAGTAACGAcgc          <  1:1064072/145‑1 (MQ=255)
                                                                                                                                            cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGAcgc          >  1:375619/1‑149 (MQ=255)
                                                                                                                                                tcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGcac      >  2:464936/1‑149 (MQ=255)
                                                                                                                                                 cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGcaca     >  1:679627/1‑149 (MQ=255)
                                                                                                                                                 cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGcaca     >  1:1105334/1‑149 (MQ=255)
                                                                                                                                                    agagTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaat  <  1:1237579/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCGTCGAAAGGCCGCCTGACGGGCTACAAATTCTACGAAGCGGTGATGAGCTCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAAT  >  NC_000913/4407729‑4408025

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: