Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,074 A→T 17.3% noncoding (84/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,0740AT17.3% 413.4 / 73.8 168noncoding (84/1542 nt)rrsG16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base A (63/75);  new base T (17/12);  total (80/88)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

Only 100 of 174 total aligned reads displayed.

ATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCA  >  NC_000913/2730975‑2731173
                                                                                                   |                                                                                                   
attattCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTCCTCCCCCGTCCGCCACTCGTCAGCGAAGCAGCaa                                                                                                     <  1:1122877/100‑1 (MQ=255)
  tATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTATTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgc                                                                                                   >  1:1717965/1‑100 (MQ=255)
   ctGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACCCGTCCGCGAAGCAGCAAgct                                                                                                  >  2:1595017/2‑100 (MQ=255)
     gCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgc                                                                                                <  1:2094421/100‑1 (MQ=38)
      cGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgct                                                                                               <  1:3940750/100‑1 (MQ=40)
        gTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTc                                                                                             >  1:1452248/1‑100 (MQ=255)
          aTTAGCTACCGTTTCCAGTAGTTATCGCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgc                                                                                           <  1:3497547/100‑3 (MQ=14)
          aTTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgc                                                                                           <  1:6252/100‑3 (MQ=255)
          aTTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgc                                                                                           <  1:203273/100‑3 (MQ=255)
           ttAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgct                                                                                          <  2:892695/100‑4 (MQ=14)
           ttAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgct                                                                                          <  1:331355/100‑4 (MQ=14)
           ttAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTg                                                                                          <  1:3137874/100‑1 (MQ=255)
            tAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCTTGCTGCTTcgctg                                                                                         <  2:2794259/100‑5 (MQ=14)
            tAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGt                                                                                         <  2:3817506/100‑1 (MQ=255)
              gCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTa                                                                                       >  1:4113125/1‑100 (MQ=255)
              gCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTa                                                                                       >  2:3371842/1‑100 (MQ=255)
               cTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAc                                                                                      <  2:1362893/100‑1 (MQ=255)
               cTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAc                                                                                      >  2:2094340/1‑100 (MQ=255)
                tACCGTTTCCATTAGTTATCCCCCTCCATCAGGCACCTTCCCAGCCCTTACCCACCCGTCCGCCACTCGTCAGGAAGTAGCTTGCTGCATcgctgacgag                                                                                    <  2:1815667/100‑10 (MQ=2)
                tACCGTTTCCAGTAGTTATCCCCCTCCATTAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAcc                                                                                     <  2:2522313/100‑1 (MQ=255)
                tACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAcc                                                                                     <  1:1559089/100‑1 (MQ=255)
                  ccGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGt                                                                                   >  1:3798723/1‑100 (MQ=255)
                   cGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGtt                                                                                  >  1:245749/1‑100 (MQ=255)
                     tttCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                                                                                <  2:494762/100‑1 (MQ=255)
                      ttCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                                                                                >  2:2329001/1‑99 (MQ=255)
                       tCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGAc                                                                              >  1:565735/1‑100 (MQ=255)
                       tCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGAc                                                                              >  2:1933362/1‑100 (MQ=255)
                       tCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGAc                                                                              >  1:1414297/1‑100 (MQ=255)
                        ccAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCCGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACt                                                                             >  1:2951336/1‑100 (MQ=255)
                        ccAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACt                                                                             <  1:1350576/100‑1 (MQ=255)
                        ccAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGCCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACt                                                                             >  1:927120/1‑100 (MQ=255)
                         cAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACtt                                                                            >  1:2801915/1‑100 (MQ=255)
                            tagtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCa                                                                         <  1:1907042/100‑1 (MQ=255)
                            tagtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCa                                                                         <  1:3035168/100‑1 (MQ=40)
                            tagtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCa                                                                         <  2:2851174/100‑1 (MQ=40)
                             agtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGc                                                                          >  2:3227037/1‑98 (MQ=255)
                             agtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAt                                                                        >  2:2988776/1‑100 (MQ=255)
                              gtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtg                                                                       <  2:2797397/100‑1 (MQ=255)
                              gtTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtg                                                                       >  1:1956047/1‑100 (MQ=255)
                              gctATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtg                                                                       <  2:2546774/98‑1 (MQ=255)
                                  tCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGtt                                                                   <  1:4103799/100‑1 (MQ=255)
                                   cccccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTa                                                                  >  1:1794174/1‑100 (MQ=255)
                                   cccccTCCATCAGGCAGCTTCCCAGACATTACGCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGtt                                                                   <  2:2047083/99‑1 (MQ=255)
                                    tcccTCCATCAGGCAGCTTCCCAGACATTACTTACCCGTCCGCCACTCCTAAGCGAAGCAGCAAGCAGCTGCCTGGTACCGTTCGACTTGCATTTGtt                                                                   >  2:3777986/2‑98 (MQ=14)
                                     cccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAgg                                                                <  1:2905294/100‑1 (MQ=255)
                                     cccTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAgg                                                                <  1:645376/100‑1 (MQ=255)
                                          cATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgc                                                           <  1:946079/100‑1 (MQ=255)
                                           aTCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgcc                                                          <  1:2244368/100‑1 (MQ=255)
                                                 cAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGc                                                    <  2:2928760/100‑1 (MQ=40)
                                                  aGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCg                                                   <  1:2360817/100‑1 (MQ=255)
                                                  aGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCg                                                   <  2:2188240/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                                              >  2:1222153/1‑100 (MQ=35)
                                                           gACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTg                                          <  1:1010839/100‑1 (MQ=38)
                                                           gACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTg                                          >  1:2444510/1‑100 (MQ=38)
                                                             cATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAg                                        <  1:1595021/100‑1 (MQ=38)
                                                             cATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAg                                        <  1:1196396/100‑1 (MQ=35)
                                                              aTTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGc                                       <  2:1761228/100‑1 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCagg                                   <  1:2459144/100‑3 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                   <  2:4170988/100‑1 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                   >  2:3581859/1‑100 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                   >  1:2863174/1‑100 (MQ=38)
                                                                  cTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATg                                   >  1:3265071/1‑100 (MQ=38)
                                                                   tCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGa                                  <  2:2406876/100‑1 (MQ=38)
                                                                    cACCCGTCCGCCATTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGAt                                 >  2:1317420/1‑100 (MQ=38)
                                                                    cACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCAGGAt                                 <  1:3982508/100‑1 (MQ=38)
                                                                      cccGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCa                               <  1:937601/100‑1 (MQ=35)
                                                                            cGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAActct                         <  2:1991309/100‑1 (MQ=38)
                                                                                acggtaacaggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCaa                     >  2:3134219/11‑100 (MQ=11)
                                                                                 cTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAt                    <  2:2533304/100‑1 (MQ=38)
                                                                                   cGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAttt                  >  2:4117223/1‑100 (MQ=38)
                                                                                    gTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTa                 >  2:2243886/1‑100 (MQ=38)
                                                                                    gTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTa                 >  1:2099098/1‑100 (MQ=38)
                                                                                      cAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGTCTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaa                >  2:4022101/1‑99 (MQ=25)
                                                                                      cAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaa               >  1:426359/1‑100 (MQ=38)
                                                                                      acaggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaa               >  2:3719441/5‑100 (MQ=255)
                                                                                        gCGAAGCAGCACGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg             <  2:3117458/100‑1 (MQ=38)
                                                                                        gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg             <  1:3547014/100‑1 (MQ=38)
                                                                                        gCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAg             >  2:3260256/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:1144018/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt            >  2:2711700/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt            >  1:3429381/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt            >  1:1098698/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt            >  1:3949516/2‑100 (MQ=255)
                                                                                         ggAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt            >  1:532664/2‑100 (MQ=255)
                                                                                         cgcagcagcAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCCCACTCCTCCATTTaaaa              >  2:2424507/4‑98 (MQ=11)
                                                                                         cGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGt            >  2:1140964/1‑100 (MQ=38)
                                                                                          gAAGCAGCTTGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGtt           >  2:2075843/1‑100 (MQ=255)
                                                                                          gAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGtt           >  1:3481668/1‑100 (MQ=38)
                                                                                               agcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATg      >  2:3845812/1‑100 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc     >  2:2986710/1‑100 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc     >  1:3629015/1‑100 (MQ=33)
                                                                                                gcaAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc     <  1:1880943/100‑1 (MQ=33)
                                                                                                   aGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTTTTAAATTGAAGAGTTTGAtcatgg  <  1:2776523/100‑6 (MQ=11)
                                                                                                   aGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa  >  1:133198/1‑100 (MQ=32)
                                                                                                   |                                                                                                   
ATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCA  >  NC_000913/2730975‑2731173

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