Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,073 A→T 19.1% noncoding (85/1542 nt) rrsG ← 16S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,0730AT19.1% 371.4 / 74.9 158noncoding (85/1542 nt)rrsG16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base A (60/67);  new base T (17/13);  total (77/81)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 172 total aligned reads displayed.

ATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTC  >  NC_000913/2730975‑2731172
                                                                                                  |                                                                                                   
attattCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTCCTCCCCCGTCCGCCACTCGTCAGCGAAGCAGCaa                                                                                                    <  1:1122877/100‑1 (MQ=255)
   aTGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTGCTCACACGTCCGCCTCTCGTCAGCGAAGCAGCAAgct                                                                                                 >  2:2823225/1‑100 (MQ=255)
     gCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgc                                                                                               <  1:2094421/100‑1 (MQ=38)
     gCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCAACTCGTCAGCGAAGCAGCAAgctgc                                                                                               >  1:606127/1‑100 (MQ=255)
      cGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCGGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgct                                                                                              >  2:2169556/1‑100 (MQ=255)
        gTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTc                                                                                            >  2:1196395/1‑100 (MQ=40)
          aTTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgc                                                                                          <  1:203273/100‑3 (MQ=255)
          aTTAGCTACCGGTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTcc                                                                                           >  1:1447378/1‑99 (MQ=255)
           ttAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgct                                                                                         <  1:331355/100‑4 (MQ=14)
           ttAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgct                                                                                         <  2:892695/100‑4 (MQ=14)
           ttAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTg                                                                                         <  1:3137874/100‑1 (MQ=255)
           ttACCTACCGTTTCCAGTAGTTATCCCCCTTCATCAGGCGGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgct                                                                                         <  2:2915032/100‑4 (MQ=14)
            tAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgctg                                                                                        <  2:2794259/100‑5 (MQ=14)
            tAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGt                                                                                        <  2:3817506/100‑1 (MQ=255)
             aGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGGAAGCAGCTTGCTGCTTcgctgac                                                                                      <  1:3963652/100‑7 (MQ=14)
              gCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTa                                                                                      >  2:3371842/1‑100 (MQ=255)
               cTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAc                                                                                     <  2:1362893/100‑1 (MQ=255)
                ttccGTTTCCAGTAGTTATCCCCCTCCATCGGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCCGCTTCCTGTTan                                                                                     <  2:2029865/97‑2 (MQ=255)
                tACCGTTTCCATTAGTTATCCCCCTCCATCAGGCACCTTCCCAGCCCTTACCCACCCGTCCGCCACTCGTCAGGAAGTAGCTTGCTGCATcgctgacgag                                                                                   <  2:1815667/100‑10 (MQ=2)
                tACCGTTTCCAGTAGTTATCCCCCTCCATTAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAcc                                                                                    <  2:2522313/100‑1 (MQ=255)
                tACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAcc                                                                                    <  1:1559089/100‑1 (MQ=255)
                   cGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGtt                                                                                 <  1:1931753/100‑1 (MQ=255)
                     tttCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                                                                               <  2:494762/100‑1 (MQ=255)
                      ttCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                                                                               >  2:2329001/1‑99 (MQ=255)
                       tCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGAc                                                                             <  1:2240639/100‑1 (MQ=255)
                        ccAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACt                                                                            <  1:1350576/100‑1 (MQ=255)
                         cAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACtt                                                                           >  1:2801915/1‑100 (MQ=255)
                           gtagtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGc                                                                         <  1:486840/100‑1 (MQ=255)
                            tagtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCa                                                                        <  1:3035168/100‑1 (MQ=40)
                             agtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGc                                                                         >  2:3227037/1‑98 (MQ=255)
                             agtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAt                                                                       >  1:982301/1‑100 (MQ=40)
                                tATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtgtg                                                                    >  1:894535/1‑100 (MQ=255)
                                tATCCCCCTCCAGCAGCCAGCTGCCCACACCTTCCTCGCCCGTCCGCCACTCGTCCGCGAAGCAGCGAGCTGCTTCCTGGTACCGTTCGACTTGc                                                                         <  2:3366759/95‑1 (MQ=14)
                                    ccccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAg                                                                <  1:3800824/100‑1 (MQ=40)
                                     cccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAgg                                                               <  1:645376/100‑1 (MQ=255)
                                     cccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAgg                                                               <  1:2905294/100‑1 (MQ=255)
                                     cccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAg                                                                <  2:3133691/99‑1 (MQ=40)
                                       cTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGc                                                              <  2:2417912/99‑1 (MQ=255)
                                         ccAGCAGGCAGCTGCCCAGACATTACTCACCCGGCCGCCACTCGTCAGCGAAGCAGCAAGCTGCGCCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTg                                                           <  2:1539249/100‑1 (MQ=14)
                                          cATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgc                                                          <  1:1351709/100‑1 (MQ=255)
                                          cATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgc                                                          <  1:946079/100‑1 (MQ=255)
                                            tCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgccg                                                        >  2:2048715/1‑100 (MQ=40)
                                               ggCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCCCCa                                                     <  2:2488689/100‑1 (MQ=255)
                                                gCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCa                                                     <  2:1978152/99‑1 (MQ=255)
                                                 cAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGc                                                   <  2:2928760/100‑1 (MQ=40)
                                                  aGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCg                                                  <  1:2360817/100‑1 (MQ=255)
                                                   gCTTCCCAGACAGTACGCCCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTGACCGTTCGACTTGCATGTGTTAGGCCTgccgcc                                                      <  2:502522/95‑1 (MQ=14)
                                                     ttCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTc                                               <  1:2962537/100‑1 (MQ=38)
                                                       cccAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCaa                                             <  1:1042591/100‑1 (MQ=38)
                                                       cccAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCaa                                             >  1:1678050/1‑100 (MQ=38)
                                                       cccAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCaa                                             >  1:2919626/1‑100 (MQ=38)
                                                          aGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCt                                          <  2:3081697/100‑1 (MQ=38)
                                                           gACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTg                                         <  1:1010839/100‑1 (MQ=38)
                                                           gACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTg                                         >  1:2444510/1‑100 (MQ=38)
                                                           gACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTg                                         <  1:1798885/100‑1 (MQ=38)
                                                             cATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAg                                       <  1:1196396/100‑1 (MQ=35)
                                                               ttACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGcc                                     <  1:1786333/100‑1 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                  >  1:2863174/1‑100 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                  <  1:930857/100‑1 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                  <  1:2510065/100‑1 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                  >  2:3581859/1‑100 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                  <  2:4170988/100‑1 (MQ=38)
                                                                    cACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGAt                                >  2:2770458/1‑100 (MQ=38)
                                                                     aCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATc                               >  1:4068893/1‑100 (MQ=35)
                                                                           ccGACACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCAAGCGTTCAATCTGAGCCATGATCAAAct                          >  1:3445891/1‑99 (MQ=255)
                                                                                acggtaacaggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCaa                    >  2:3134219/11‑100 (MQ=11)
                                                                                 cTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAt                   <  2:2533304/100‑1 (MQ=38)
                                                                                  ggtaacaggAAGCAGCTTGCTGCTTCCTGATACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAtt                  >  1:2457514/9‑100 (MQ=11)
                                                                                   cGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGATTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAttt                 >  1:2020122/1‑100 (MQ=255)
                                                                                   cGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAttt                 >  1:3243098/1‑100 (MQ=38)
                                                                                   cGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAttt                 >  2:4117223/1‑100 (MQ=38)
                                                                                    gTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTa                >  2:2243886/1‑100 (MQ=38)
                                                                                      cAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGTCTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaa               >  2:4022101/1‑99 (MQ=25)
                                                                                      cAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaa              >  1:426359/1‑100 (MQ=38)
                                                                                      cAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaa              >  1:1502990/1‑100 (MQ=38)
                                                                                      acaggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaa              >  2:3719441/5‑100 (MQ=255)
                                                                                       caggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaaa             >  1:4112593/4‑100 (MQ=11)
                                                                                        gCGAATCAGTAAGCCGCGTCCGGTTACCGTTCCGCTGGCATGTGTTAGGCCGGCCACCAGCGTTCAAGCTGAGCCATGATCAAACGCTTCAACTTaaa              <  2:1809615/98‑1 (MQ=2)
                                                                                        gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg            >  2:2214663/1‑100 (MQ=38)
                                                                                        gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg            <  1:3547014/100‑1 (MQ=38)
                                                                                        gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg            >  2:2579532/1‑100 (MQ=38)
                                                                                        gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg            >  1:1076858/1‑100 (MQ=38)
                                                                                        aggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg            >  2:4034087/3‑100 (MQ=255)
                                                                                        aggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATGTAAAAg            >  2:615018/3‑100 (MQ=11)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  1:3429381/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  1:3494472/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  1:1098698/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  1:1144018/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  1:3949516/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  1:554443/2‑100 (MQ=255)
                                                                                         cGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt           >  2:3817234/1‑100 (MQ=38)
                                                                                          gAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGtt          >  2:2075843/1‑100 (MQ=255)
                                                                                           aaGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGttt         <  1:375462/100‑1 (MQ=38)
                                                                                               agcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATg     >  2:3845812/1‑100 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc    >  1:3629015/1‑100 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc    >  2:2986710/1‑100 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc    <  2:1288994/100‑1 (MQ=32)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc    <  2:1725504/100‑1 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc    <  1:1880943/100‑1 (MQ=33)
                                                                                                  aaGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTc  <  1:2283990/100‑1 (MQ=32)
                                                                                                  |                                                                                                   
ATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTC  >  NC_000913/2730975‑2731172

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