Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,820,565 A→G 6.1% I217T (ATC→ACC)  chbC ← N,N'‑diacetylchitobiose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,820,5650AG 6.1% 191.3 / 3.1 65I217T (ATC→ACC) chbCN,N'‑diacetylchitobiose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base G (2/2);  total (31/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 = 3.44e-01

CCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACC  >  NC_000913/1820417‑1820702
                                                                                                                                                    |                                                                                                                                         
ccACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGATCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCCACGATGCCAGTGGGCTTGAGATGGTATCCCTAATGa                                                                                                                                           >  1:117410/1‑149 (MQ=255)
 taCGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGAt                                                                                                                                          >  7:119226/2‑149 (MQ=255)
  aCGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATc                                                                                                                                         <  1:364917/149‑1 (MQ=255)
    ggCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCt                                                                                                                                        >  4:269560/1‑148 (MQ=255)
     gCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCTAGTGGGGTTGAGATGGTATCCATAATGATCTGa                                                                                                                                      <  4:343535/149‑1 (MQ=255)
          aTAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATgg                                                                                                                                   >  5:289101/1‑147 (MQ=255)
          aTAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGaa                                                                                                                                 <  5:52184/149‑1 (MQ=255)
             aTGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGt                                                                                                                               <  8:119226/148‑1 (MQ=255)
               gCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTgg                                                                                                                            >  7:79098/1‑149 (MQ=255)
                     gTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAAGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGccc                                                                                                                       <  3:258426/148‑1 (MQ=255)
                        cAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAg                                                                                                                    <  2:191228/148‑1 (MQ=255)
                                     gccagcgcGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCa                                                                                                      <  1:20087/149‑1 (MQ=255)
                                                 tggaTACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATccc                                                                                          >  4:128469/4‑149 (MQ=255)
                                                  gAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATccc                                                                                          >  4:118411/1‑148 (MQ=255)
                                                       ccGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCAc                                                                                    <  5:37031/149‑1 (MQ=255)
                                                        cGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGGTCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACg                                                                                   >  7:61018/1‑149 (MQ=255)
                                                        cGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCGCATCACg                                                                                   >  7:131599/1‑149 (MQ=255)
                                                        cGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACg                                                                                   >  3:149893/1‑149 (MQ=255)
                                                        agaagaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGGTCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACg                                                                                   >  3:188156/2‑149 (MQ=255)
                                                            gaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGaaa                                                                               <  7:69159/149‑1 (MQ=255)
                                                            gaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCAACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATACTGATCTGATGGAAGTTGGTGCCCCATGTATTCAACGCCCAGGCAATAATCCCCATCACGGaaa                                                                               <  4:45275/149‑1 (MQ=255)
                                                               cctgagCAGTGGAAAAAAGATCACATAGGCCCTGCCCACCACGCTACGCAACGATGCCAGTTGGGTTGAGATGGTATCCATAATGGTCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGaa                                                                            <  4:188156/146‑1 (MQ=255)
                                                               ccAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGGTCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGaa                                                                            <  8:61018/149‑1 (MQ=255)
                                                               ccAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTg                                                                                                                             >  2:372867/1‑100 (MQ=255)
                                                               ccAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTg                                                                                                                             <  1:372867/100‑1 (MQ=255)
                                                                cAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGaa                                                                            >  1:334883/1‑148 (MQ=255)
                                                                      aGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATaaa                                                                     <  7:281904/149‑1 (MQ=255)
                                                                        tGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTc                                                                                                              <  1:193096/106‑1 (MQ=255)
                                                                        tGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTc                                                                                                              >  2:193096/1‑106 (MQ=255)
                                                                        tGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGCTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTAATCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGc                                                                   <  6:289101/149‑1 (MQ=255)
                                                                           aaCAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCgg                                                                <  2:334883/149‑1 (MQ=255)
                                                                            aCAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCgg                                                                >  7:353224/1‑148 (MQ=255)
                                                                             cAAAGATCACATAGGCCCAGCCCACCCCGCTACCCATCGATGCCAGTGGGGTTGAGATGGTATCCATAATGAGCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCggg                                                               >  3:277748/1‑148 (MQ=255)
                                                                             cAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAg                                                                                                                    >  4:45182/1‑95 (MQ=255)
                                                                              aaaGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGaa                                                             >  3:290088/1‑149 (MQ=255)
                                                                              aaaGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGaa                                                             >  6:351201/1‑149 (MQ=255)
                                                                                aGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAAt                                                            >  6:72785/1‑148 (MQ=255)
                                                                                   tCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAAt                                                        <  7:147652/149‑1 (MQ=255)
                                                                                       aTAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGcc                                                     >  4:264986/1‑148 (MQ=255)
                                                                                        tAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCga                                                   <  3:269560/149‑1 (MQ=255)
                                                                                          ggCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCgag                                                  <  8:79098/148‑1 (MQ=255)
                                                                                           gcccagcccaCCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGaa                                                <  8:131599/149‑1 (MQ=255)
                                                                                                    accacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGa                                       <  5:351201/149‑1 (MQ=255)
                                                                                                    accacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGa                                       <  4:290088/149‑1 (MQ=255)
                                                                                                     ccacGCTACCCATCGATTCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGTCAATAATACCCATCACGGAAAGAATAATAAAGCCGGGaa                                                             <  7:5160/126‑1 (MQ=255)
                                                                                                     ccacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGAt                                      <  7:205552/149‑1 (MQ=255)
                                                                                                     ccacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCACATCAAGGAAAGAATAATAAAGCCGGGaa                                                             >  8:5160/1‑126 (MQ=255)
                                                                                                      cacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAAAGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATa                                     >  6:147278/1‑149 (MQ=255)
                                                                                                             cccAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGc                              <  3:264986/149‑1 (MQ=255)
                                                                                                              ccAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAg                               >  5:109304/1‑147 (MQ=255)
                                                                                                               cAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAg                                               >  7:256760/1‑130 (MQ=255)
                                                                                                               cAAAGATGCGAGTGGGGTTGAGATGTTTTCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAg                                               <  8:256760/130‑1 (MQ=255)
                                                                                                                     tGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGc                      >  5:124355/1‑149 (MQ=255)
                                                                                                                      gCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATATTCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCAATACTGAAGCATGTACGCt                     >  1:213263/1‑149 (MQ=255)
                                                                                                                         aGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTc                  <  6:109304/149‑1 (MQ=255)
                                                                                                                         aGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTc                  <  3:118411/149‑1 (MQ=255)
                                                                                                                           tGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTaa                                                         >  1:233455/1‑108 (MQ=255)
                                                                                                                           tGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTaa                                                         <  2:233455/108‑1 (MQ=255)
                                                                                                                           tGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCgg                <  5:72785/149‑1 (MQ=255)
                                                                                                                             cggTTGAGATCGTATACATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCACTAATCACCATCACGGAAAGAATAATAAAGCCGGGAATTAATGAGGACAAGCAACCAGATACTGAAGCAGGTACGCTGTCggc               >  4:211826/2‑147 (MQ=255)
                                                                                                                                  gagaTGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGtt          <  8:353224/148‑1 (MQ=255)
                                                                                                                                     aTGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCGGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAAt      >  3:324740/1‑149 (MQ=255)
                                                                                                                                        gTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACTGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGa    <  5:147278/148‑1 (MQ=255)
                                                                                                                                         tATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGGTTAATGAcc  >  4:93930/1‑149 (MQ=255)
                                                                                                                                              aTAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCa                                                                                     >  1:31272/1‑61 (MQ=255)
                                                                                                                                              aTAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCa                                                                                     <  2:31272/61‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
CCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACC  >  NC_000913/1820417‑1820702

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: ‑