Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,836 A→G 5.2% T188A (ACC→GCC)  rlmB → 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,8360AG 5.2% 230.3 / 3.0 78T188A (ACC→GCC) rlmB23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (2/2);  total (40/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 = 7.26e-01

AAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTT  >  NC_000913/4409689‑4409982
                                                                                                                                                   |                                                                                                                                                  
aaaaaaGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGa                                                                                                                                                    >  8:260598/1‑148 (MQ=255)
         gCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGcc                                                                                                                                           <  5:257982/148‑1 (MQ=255)
          ccTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCt                                                                                                                                          <  4:288380/148‑1 (MQ=255)
               ggcgcggcAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGgcgc                                                                                                                                     >  1:355550/1‑148 (MQ=255)
                gcgcggcAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCt                                                                                                                                    <  7:260598/148‑1 (MQ=255)
                 cgcggcAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTg                                                                                                                                   <  7:8923/148‑1 (MQ=255)
                    ggcAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTg                                                                                                                                <  4:264571/148‑1 (MQ=255)
                    ggcAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTg                                                                                                                                <  8:256506/148‑1 (MQ=255)
                        gAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATg                                                                                                                             <  1:371089/147‑1 (MQ=255)
                        gAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATg                                                                                                                             >  2:371089/1‑147 (MQ=255)
                          aaGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGt                                                                                                                          >  5:277966/1‑148 (MQ=255)
                                  cGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGgaag                                                                                                                   <  2:49308/147‑1 (MQ=255)
                                       aTTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGgaaggtgaa                                                                                                              <  5:223584/147‑1 (MQ=255)
                                          cGGGTGATTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGgaaggtgaaggt                                                                                                           <  1:369532/147‑1 (MQ=255)
                                            ggTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAATATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATg                                                                                                        <  8:292237/148‑1 (MQ=255)
                                             gTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGc                                                                                                       >  6:222977/1‑148 (MQ=255)
                                                aCTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTc                                                                                                    >  7:205699/1‑148 (MQ=255)
                                                aCTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCGGGTGAGGGGTGCGGAAGGTGAAGGTATGCGTc                                                                                                    <  4:276608/148‑1 (MQ=255)
                                                 cTAACATTGCGCGCACCATGCGTATGTTGCAGGAAGAGATGTTGTGGATCGTCGGTACGGCAGGCGAGGCGGAGCAGTCGCTCTATCAGAGCGAAATGCCCGGACGCGTGGCGCTGGTGATGGTTGGGGAAGGTGAAGGTATGCGACg                                                                                                   >  8:147920/1‑148 (MQ=255)
                                                   aaCCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGcc                                                                                                 >  7:282276/1‑148 (MQ=255)
                                                   aaCCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGcc                                                                                                 >  5:12339/1‑148 (MQ=255)
                                                    aCCTTGCGCGCACCATGCGTATGTTGAAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGc                                                                                                  <  7:190168/146‑1 (MQ=255)
                                                     ccTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATCTCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTg                                                                                               <  1:179726/148‑1 (MQ=255)
                                                        tGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACt                                                                                            >  8:232016/1‑148 (MQ=255)
                                                         gcgcgcACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACt                                                                                            <  6:277966/147‑1 (MQ=255)
                                                                  aTGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACAtt                                                                                  >  8:71961/1‑148 (MQ=255)
                                                                  aTGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTAGCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGCCTCGGGAACAtt                                                                                  >  5:196619/1‑148 (MQ=255)
                                                                          gTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAg                                                                          >  3:325506/1‑148 (MQ=255)
                                                                           ttGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTg                                                                                                                                >  7:247866/1‑93 (MQ=255)
                                                                           ttGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTg                                                                                                                                >  7:253716/1‑93 (MQ=255)
                                                                           ttGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTg                                                                                                                                <  8:253716/93‑1 (MQ=255)
                                                                           ttGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTg                                                                                                                                <  8:247866/93‑1 (MQ=255)
                                                                              cAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACAGTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAACGTATGCGTCGCCGGACTCGTGAACATTGCGATGAGTTGa                                                                      >  2:222009/1‑148 (MQ=255)
                                                                                ggAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATc                                                                    <  8:272349/148‑1 (MQ=255)
                                                                                ggAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATc                                                                    <  2:20048/148‑1 (MQ=255)
                                                                                 gAAGAGAATATCTGGATCGTCGGTACGGCAGGAGAGGCGGATCATACACACTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGAt                                                                     >  5:328368/1‑146 (MQ=255)
                                                                                  aaGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCTAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAg                                                                  >  3:310507/1‑148 (MQ=255)
                                                                                  aaGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAg                                                                  >  2:102791/1‑148 (MQ=255)
                                                                                   agagAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAg                                                                  >  2:242960/1‑147 (MQ=255)
                                                                                    gagaATATCTGGATCGTCGGGACGGCCGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGGGGGTGATGGGTGGGGAAGGGGTAGGTATGCGTGGCCTGCCTCGGGTACATTGGGATGAGTTGATCAGcg                                                                >  6:110093/1‑147 (MQ=255)
                                                                                         tatCTTGATCGTCGGTACGGCAGGCGAGGCGGATCATACAATCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTTAACATTGCGATGAGTTGATCAGCATCCCg                                                           <  6:328368/148‑1 (MQ=255)
                                                                                            cTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATg                                                        >  1:211174/1‑148 (MQ=255)
                                                                                               gATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGGCCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCt                                                     >  8:58992/1‑148 (MQ=255)
                                                                                               gATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGGCCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCt                                                     >  8:65098/1‑148 (MQ=255)
                                                                                               gATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCt                                                     <  4:269410/148‑1 (MQ=255)
                                                                                                aTCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTg                                                    <  1:45320/148‑1 (MQ=255)
                                                                                                aTCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTg                                                    <  5:282913/148‑1 (MQ=255)
                                                                                                 tcgtcgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGCCGCCTGGCGCTGGTGATGAGTGCGGAAGGTGTAGGTATGCGTCGCCGGAGTCGGGAACATTTCGATGAGTTGATCAGCATCCCGATGGCTgg                                                   >  2:353260/1‑148 (MQ=255)
                                                                                                  cgtcgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCg                                                                               >  3:155070/1‑119 (MQ=255)
                                                                                                  cgtcgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCg                                                                               <  4:155070/119‑1 (MQ=255)
                                                                                                     cgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGc                                               >  6:349277/1‑148 (MQ=255)
                                                                                                     cgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGc                                               >  3:189602/1‑148 (MQ=255)
                                                                                                       gTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGt                                             >  6:127775/1‑148 (MQ=255)
                                                                                                        tACGGCAGGGGAGGCGGATCATACTCTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGCTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGtt                                            <  4:310507/148‑1 (MQ=255)
                                                                                                          cGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTttc                                          <  2:211174/148‑1 (MQ=255)
                                                                                                            gCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTttctt                                        <  8:282276/148‑1 (MQ=255)
                                                                                                            gCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTttctt                                        >  3:155534/1‑148 (MQ=255)
                                                                                                             cAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTttcttc                                       <  1:102791/148‑1 (MQ=255)
                                                                                                                gcgaggcgGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCAGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCt                                    >  6:370939/1‑148 (MQ=255)
                                                                                                                   aggcgGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGaa                                 <  5:222977/148‑1 (MQ=255)
                                                                                                                         aTCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTc                           >  6:84188/1‑148 (MQ=255)
                                                                                                                            aTACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGt                        >  6:144256/1‑148 (MQ=255)
                                                                                                                            aTACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGt                        <  4:189602/148‑1 (MQ=255)
                                                                                                                              acacTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGGGTCGCCTGCCTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTg                      >  1:136630/1‑148 (MQ=255)
                                                                                                                              acacTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTg                      >  1:301728/1‑148 (MQ=255)
                                                                                                                               cacTCTATCAGAGCAAAATGGCCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTg                      <  7:65098/147‑1 (MQ=255)
                                                                                                                               cacTCTATCAGAGCAAAATGGCCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGGCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTg                      <  7:58992/147‑1 (MQ=255)
                                                                                                                               cacTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCa                                                                <  1:51052/105‑1 (MQ=255)
                                                                                                                               cacTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCa                                                                >  2:51052/1‑105 (MQ=255)
                                                                                                                                 ctctATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGa                   <  7:71961/148‑1 (MQ=255)
                                                                                                                                   ctATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGAcc                 >  2:2826/1‑148 (MQ=255)
                                                                                                                                     aTCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGCCCg                >  3:183638/1‑147 (MQ=255)
                                                                                                                                        agagCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGaa             <  1:2826/147‑1 (MQ=255)
                                                                                                                                         gagCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAAtt           <  5:349277/148‑1 (MQ=255)
                                                                                                                                         gagCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAAtt           >  7:152817/1‑148 (MQ=255)
                                                                                                                                          agCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGACCGTTTCGGTTGCGACCGGAAttt          >  6:266355/1‑148 (MQ=255)
                                                                                                                                            cAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGc        <  5:370939/148‑1 (MQ=255)
                                                                                                                                                aTGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTAt    >  3:209702/1‑148 (MQ=255)
                                                                                                                                                aTGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTAt    <  7:119674/148‑1 (MQ=255)
                                                                                                                                                   aCCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTAttt  <  4:209702/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTT  >  NC_000913/4409689‑4409982

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