Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,606,548 A→T 5.3% H150L (CAT→CTT)  lsrF → putative autoinducer‑2 (AI‑2) aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,606,5480AT 5.3% 244.0 / 3.6 75H150L (CAT→CTT) lsrFputative autoinducer‑2 (AI‑2) aldolase
Reads supporting (aligned to +/- strand):  ref base A (39/32);  new base T (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

ATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGG  >  NC_000913/1606401‑1606692
                                                                                                                                                   |                                                                                                                                                
aTAGCCAGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAACTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTTCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACa                                                                                                                                                  >  1:9908/1‑148 (MQ=255)
        gtggtACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGa                                                                                                                                          <  7:297044/148‑1 (MQ=255)
          ggtACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGAt                                                                                                                                         >  6:407481/1‑147 (MQ=255)
          ggtACTGCGGGCCTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGAt                                                                                                                                         >  2:193634/1‑147 (MQ=255)
              cTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCaaa                                                                                                                                     >  1:138740/1‑147 (MQ=255)
              cTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCaaa                                                                                                                                     >  8:272429/1‑147 (MQ=255)
                 cGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGGGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTGTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAat                                                                                                                                  <  2:32279/147‑1 (MQ=255)
                        cAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATattatt                                                                                                                            >  6:381002/1‑146 (MQ=255)
                          ggTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCAGTACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGGAGCGAAGATGAACAGCAGTCGATCAAAAATAGTATTCAGc                                                                                                                        >  6:338554/1‑148 (MQ=255)
                            tGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACTTCAGTCGATCAAAAATattatt                                                                                                                            <  3:151814/142‑1 (MQ=255)
                            tGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTg                                                                                                                      <  1:193634/148‑1 (MQ=255)
                            tGCGAACTCTAGTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACTTCAGTCGATCAAAAATattatt                                                                                                                            >  4:151814/1‑142 (MQ=255)
                             gCGAACTCAATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGCTCGCGGCGCAGGTATATATCGGCAGCGAAAATGAACATCAGTCGATCAAAAATATTATTCAGcatg                                                                                                                     >  4:360968/1‑145 (MQ=255)
                              cGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGt                                                                                                                    >  8:325520/1‑148 (MQ=255)
                                aaCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGtt                                                                                                                   >  5:286452/1‑147 (MQ=255)
                                aaCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTg                                                                                                                  <  8:143217/148‑1 (MQ=255)
                                aaCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTg                                                                                                                  <  5:407481/148‑1 (MQ=255)
                                 aCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTgg                                                                                                                     >  7:366193/1‑144 (MQ=255)
                                 aCTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTgg                                                                                                                     <  8:366193/144‑1 (MQ=255)
                                    ctATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGc                                                                                                              >  7:156093/1‑148 (MQ=255)
                                        tCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCgg                                                                                                           >  7:132996/1‑147 (MQ=255)
                                        tCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACAGCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCgg                                                                                                           >  6:419464/1‑147 (MQ=255)
                                           ggCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTgg                                                                                                                     >  2:377390/1‑134 (MQ=255)
                                           ggCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTgg                                                                                                                     <  1:377390/134‑1 (MQ=255)
                                             cGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATgc                                                                                                      >  8:129222/1‑146 (MQ=255)
                                                     aGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGa                                                                                              <  1:12254/147‑1 (MQ=255)
                                                     aGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGa                                                                                              <  7:115305/147‑1 (MQ=255)
                                                        aaTGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTgtacgg                                                                                            >  5:359745/1‑141 (MQ=255)
                                                          tGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGcc                                                                                         >  7:318005/1‑147 (MQ=255)
                                                              gCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGaa                                                                                                         >  5:190973/1‑127 (MQ=255)
                                                              gCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGaa                                                                                                         <  6:190973/127‑1 (MQ=255)
                                                                cGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAAtattat                                                                                                                             <  1:244470/105‑1 (MQ=255)
                                                                cGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAAtattat                                                                                                                             >  2:244470/1‑105 (MQ=255)
                                                                cGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATg                                                                                  <  8:194905/148‑1 (MQ=255)
                                                                          tCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACTTCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACt                                                                         >  8:256595/1‑147 (MQ=255)
                                                                          tCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACTTCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACt                                                                         <  7:256595/147‑1 (MQ=255)
                                                                             atggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCg                                                                     >  1:268424/1‑148 (MQ=255)
                                                                               ggatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTg                                                                   >  5:232366/1‑148 (MQ=255)
                                                                                gatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGa                                                                                              <  5:101492/120‑1 (MQ=255)
                                                                                gatgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGa                                                                                              >  6:101492/1‑120 (MQ=255)
                                                                                 atgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTg                                                                            >  2:396640/1‑137 (MQ=255)
                                                                                 atgACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTg                                                                            <  1:396640/137‑1 (MQ=255)
                                                                                    aCGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGCACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAAGGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCa                                                               >  8:323772/1‑147 (MQ=255)
                                                                                    aCGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCa                                                               >  1:12315/1‑147 (MQ=255)
                                                                                        cGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAg                                                                                                   <  1:353184/107‑1 (MQ=255)
                                                                                        cGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAg                                                                                                   >  2:353184/1‑107 (MQ=255)
                                                                                            cgcCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATgg                                                      <  5:381002/148‑1 (MQ=255)
                                                                                                tGAACAGTTGAGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGACGGAATTAAAGTGTGAATGACGACAATTGCAGTGACTGGAGTTGGCAAAGATATGGTgcg                                                  >  8:375832/1‑148 (MQ=255)
                                                                                                  aaCAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCTAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGcgcg                                                <  2:138740/148‑1 (MQ=255)
                                                                                                  aaCAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTc                                                                                                                           >  1:254794/1‑73 (MQ=255)
                                                                                                  aaCAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTc                                                                                                                           <  2:254794/73‑1 (MQ=255)
                                                                                                   aCAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGa                                               >  8:86520/1‑148 (MQ=255)
                                                                                                     aGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCATCTAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATc                                             <  2:9908/148‑1 (MQ=255)
                                                                                                     aGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAata                                                                                                                                 <  3:303245/64‑1 (MQ=255)
                                                                                                     aGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAata                                                                                                                                 >  4:303245/1‑64 (MQ=255)
                                                                                                          cgcgGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTggg                                                                                               >  8:83346/1‑93 (MQ=255)
                                                                                                          cgcgGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTggg                                                                                               <  7:83346/93‑1 (MQ=255)
                                                                                                           gcgGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTCGGCAAAGATATGGTGCGCGATCAGCGtt                                       <  8:132996/148‑1 (MQ=255)
                                                                                                            cgGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTa                                      >  8:184056/1‑148 (MQ=255)
                                                                                                              gTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTAtt                                    <  7:325520/148‑1 (MQ=255)
                                                                                                                      gcAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTctcgctc                            <  7:272429/148‑1 (MQ=255)
                                                                                                                       cAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCTGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTctcgctcg                           >  5:323171/1‑148 (MQ=255)
                                                                                                                         ggTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCtcgctcgc                          >  4:312944/1‑147 (MQ=255)
                                                                                                                           tttATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACt                      >  7:79988/1‑148 (MQ=255)
                                                                                                                               tatCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGaa                  <  8:79988/148‑1 (MQ=255)
                                                                                                                                atCGGCAGCGAATATGAACATCTGTCGATCATAAATATTATTCAGGTGGTTGATGCCGGATTGAATGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGa                   <  6:232366/147‑1 (MQ=255)
                                                                                                                                 tCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCTGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGaa                  <  6:323171/147‑1 (MQ=255)
                                                                                                                                   ggCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGc                                                                                >  7:67109/1‑83 (MQ=255)
                                                                                                                                   ggCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGc                                                                                <  8:67109/83‑1 (MQ=255)
                                                                                                                                      aGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATcgccg            <  2:12315/148‑1 (MQ=255)
                                                                                                                                      aGCGAATATGAACATCAGGCGATCAAAAATATTATTCAGCTGGTTGCTGCCGGAATGAAAGTGGGATTGCCGACCATGGCCTTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGCCTCGAATcgcc             >  6:101364/1‑147 (MQ=255)
                                                                                                                                        cGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCt          >  5:171277/1‑148 (MQ=255)
                                                                                                                                           atatGAACATCAGGCGATCAAAAAAATTATTCAGCTGGGTGATGCCGGAATGAAAGGGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGGTATTTCTCGCTCGCGACTCGAATCGCCGCTGaa       <  2:268424/148‑1 (MQ=255)
                                                                                                                                                aaCATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATcgcc             >  5:260705/1‑137 (MQ=255)
                                                                                                                                                aaCATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATcgcc             <  6:260705/137‑1 (MQ=255)
                                                                                                                                                aaCATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATggg  <  6:286452/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTGGCGTTATCGATGGATGACGCCGTGCGCCTGAACAGTTGCGCGGTGGCGGCGCAGGTTTATATCGGCAGCGAATATGAACATCAGTCGATCAAAAATATTATTCAGCTGGTTGATGCCGGAATGAAAGTGGGAATGCCGACCATGGCCGTGACTGGCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGG  >  NC_000913/1606401‑1606692

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: