Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,987,438 T→C 7.1% D231D (GAT→GAC PP_1778 → lipopolysaccharide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,987,4380TC 7.1% 206.8 / 4.3 70D231D (GAT→GACPP_1778lipopolysaccharide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base C (2/3);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTAC  >  NC_002947/1987295‑1987584
                                                                                                                                               |                                                                                                                                                    
ctgcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCTCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGa                                                                                                                                                  <  1:306435/148‑1 (MQ=255)
ctgcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGa                                                                                                                                                  <  8:212764/148‑1 (MQ=255)
ctgcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGa                                                                                                                                                  <  6:150692/148‑1 (MQ=255)
ctgcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGa                                                                                                                                                  <  7:5709/148‑1 (MQ=255)
     ggTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAgg                                                                                                                                            <  3:183745/149‑1 (MQ=255)
     ggTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAgg                                                                                                                                            <  1:306812/149‑1 (MQ=255)
         ttgttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCgg                                                                                                                                         <  7:45946/148‑1 (MQ=255)
           gttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGtt                                                                                                                                      <  4:305788/148‑1 (MQ=255)
            ttCTGGTTTACCCCTATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCGCttta                                                                                                                                    <  6:41278/149‑3 (MQ=255)
            ttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGa                                                                                                                                                  >  5:43597/1‑136 (MQ=255)
            ttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGt                                                                                                                                   >  2:40850/1‑149 (MQ=255)
            ttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCGCttta                                                                                                                                    <  6:220072/149‑3 (MQ=255)
            ttCTGGTTTACACCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGa                                                                                                                                                  <  6:43597/136‑1 (MQ=255)
                 gTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCAt                                                                                                                                <  1:207374/147‑1 (MQ=255)
                     aCCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCg                                                                                                                           <  1:147946/148‑1 (MQ=255)
                      ccccAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCgg                                                                                                                          >  6:363509/1‑148 (MQ=255)
                          aaTCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGAttt                                                                                                                     <  5:363509/149‑1 (MQ=255)
                               tttACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCttt                                                                                                                <  2:304479/149‑1 (MQ=255)
                                     tGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCtg                                                                                                          <  1:40850/149‑1 (MQ=255)
                                               gATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTGGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTg                                                                                                >  7:100114/1‑149 (MQ=255)
                                                  aCCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcc                                                                                             >  8:196117/1‑149 (MQ=255)
                                                    cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCtgt                                                                                                         <  6:170536/135‑1 (MQ=255)
                                                    cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCtgt                                                                                                         >  5:170536/1‑135 (MQ=255)
                                                    cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCAAGACCGGGATTTCCTTTGGACTGTTGTTGCTTGGactg                                                                                           >  1:323823/1‑149 (MQ=255)
                                                     tGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgt                                                                                          >  5:263474/1‑149 (MQ=255)
                                                         cacCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCTTGATTTGAGTCAGGCGGTTGTc                                                                                                                                  <  2:342351/105‑1 (MQ=255)
                                                          acCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcctg                                                                                     >  8:112259/1‑149 (MQ=255)
                                                            cTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGAACGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcctgtt                                                                                   >  4:27464/1‑149 (MQ=255)
                                                                    cACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGt                                                                           <  7:321910/149‑1 (MQ=255)
                                                                       ccTTAGTTTTAACCCAATCTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTCTCATGAACGTGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTCTTTCGTCGGTCtg                                                                        >  4:334603/1‑149 (MQ=255)
                                                                                 aaCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAt                                                              >  1:343974/1‑149 (MQ=255)
                                                                                  aCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtg                                                             <  1:355189/149‑1 (MQ=255)
                                                                                  aCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTGGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAt                                                              <  8:100114/148‑1 (MQ=255)
                                                                                       aTGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGa                                                        >  4:303449/1‑149 (MQ=255)
                                                                                          tATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgtt                                                                                         >  3:217444/1‑113 (MQ=255)
                                                                                          tATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgtt                                                                                         <  4:217444/113‑1 (MQ=255)
                                                                                          tATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGACTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATg                                                      <  6:145241/148‑1 (MQ=255)
                                                                                               cATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTTGATCATGAACGGGATAACATTTTTTATTTTTTTGATTGGCATGTTAATGTTTAGTCTTTCTGTGCCTTAAATGTTGTATGTACt                                                  >  1:298225/1‑149 (MQ=255)
                                                                                               cATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATTACCGGGATTTGCTTTGGTCTCTTGTTGATTGGCCTGTTCCTGTTTCGTCGGTCAGTGCATGAAATGGTGGATGTACTAt                                                >  7:128261/1‑149 (MQ=255)
                                                                                                aTTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATg                                               >  8:219295/1‑149 (MQ=255)
                                                                                                      aCTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtt                                         >  7:89623/1‑149 (MQ=255)
                                                                                                        tGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                       >  2:124882/1‑149 (MQ=255)
                                                                                                          ccTTTCAGGACGTACTCCTGTACGCTTATGCTCCTGATTTTAGTCAGGCGGTTGTCTTGAGGGGGTTTTGCTTTGGTCTTTTTTTGCTTGGCCTGTTCCTTTTTCGTCGGTCTGTGCGTGTAATGGTGGATGTACTATGAGCCtttt                                       <  3:334603/147‑1 (MQ=255)
                                                                                                          ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggt                                                           >  1:245738/1‑127 (MQ=255)
                                                                                                          ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggt                                                           <  2:245738/127‑1 (MQ=255)
                                                                                                          ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                       >  6:149408/1‑147 (MQ=255)
                                                                                                          ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                       <  7:196117/147‑1 (MQ=255)
                                                                                                              tCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAgt                                   <  3:93859/147‑1 (MQ=255)
                                                                                                                acgacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCa                               <  3:315117/149‑1 (MQ=255)
                                                                                                                acgacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGAGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTAGGCCTGTTCATGTTTCGTCGGTCTGTGCGTGAAATGGTGGATGTACTATGAGACTTTTGAGTGTCa                               >  4:107084/1‑149 (MQ=255)
                                                                                                                  gacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGt                             >  8:165777/1‑149 (MQ=255)
                                                                                                                   acgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATTACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGc                                            <  4:309156/133‑1 (MQ=255)
                                                                                                                   acgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGc                                            >  3:309156/1‑133 (MQ=255)
                                                                                                                    cgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGa                           <  6:263474/149‑1 (MQ=255)
                                                                                                                     gTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGAt                          >  8:147895/1‑149 (MQ=255)
                                                                                                                         tCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTgg                      <  4:251182/149‑1 (MQ=255)
                                                                                                                           cTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTggg                     >  8:81792/1‑148 (MQ=255)
                                                                                                                           cTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTggg                     >  7:37416/1‑148 (MQ=255)
                                                                                                                           cTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGc                    <  5:149408/149‑1 (MQ=255)
                                                                                                                                cGCTAAAGCGCCTGTTTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTTCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGc               <  2:298225/149‑1 (MQ=255)
                                                                                                                                   tAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCt             >  5:2642/1‑148 (MQ=255)
                                                                                                                                     aaccgcCTGACTCAAATCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                     >  4:267664/4‑122 (MQ=255)
                                                                                                                                     aaccgcCTGACTCAAATCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                     <  3:267664/119‑1 (MQ=255)
                                                                                                                                     aaGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCt             <  6:2642/146‑1 (MQ=255)
                                                                                                                                     aaGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTAc           <  3:303449/148‑1 (MQ=255)
                                                                                                                                      aGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACa          >  7:81446/1‑148 (MQ=255)
                                                                                                                                       ccgcCTGACTCAAATCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGACTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACaa         >  4:341863/2‑148 (MQ=255)
                                                                                                                                          ccTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACg     >  6:342857/1‑149 (MQ=255)
                                                                                                                                             gATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGCCAAAGCCTACAAAACGTAc  >  1:306098/1‑149 (MQ=255)
                                                                                                                                              aTTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTAc  <  7:112259/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTGCAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTAC  >  NC_002947/1987295‑1987584

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: