Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,987,440 T→C 7.0% L232S (TTG→TCG) ‡ 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,4400TC 7.0% 248.2 / 6.9 86L232S (TTG→TCG) ‡PP_1778lipopolysaccharide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (43/37);  new base C (3/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCTGCAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTA  >  NC_002947/1987294‑1987583
                                                                                                                                                   |                                                                                                                                               
gctgcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGGCAATGTGATACCTGATCACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTCCGCTAAAGCGCCTGATTTGa                                                                                                                                               <  5:90442/149‑1 (MQ=255)
  tgcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGt                                                                                                                                             <  8:21879/149‑1 (MQ=255)
  tgcAGGTGCTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCATAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTCCTCCTGTACGCTAAAGCGCCTGATTTGAGt                                                                                                                                             >  7:5483/1‑149 (MQ=255)
   gcAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACGCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTACTGCGCCTGATTTGAGTc                                                                                                                                            >  4:236196/1‑149 (MQ=255)
         gttgttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCgg                                                                                                                                      >  6:72412/1‑149 (MQ=255)
          ttgtgCTGGATTACCCCCATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCgg                                                                                                                                      <  4:275440/148‑1 (MQ=255)
            gttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCGCtt                                                                                                                                    <  4:40699/148‑1 (MQ=255)
             ttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTg                                                                                                                                   >  7:254115/1‑148 (MQ=255)
               cTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTTATTTGAGTCAGGCGGTTGTCa                                                                                                                                >  6:374625/1‑149 (MQ=255)
               cTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTc                                                                                                                                 <  1:62644/148‑1 (MQ=255)
                        cccAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTGACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGa                                                                                                                       >  1:124292/1‑149 (MQ=255)
                        cccAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGa                                                                                                                       >  7:147825/1‑149 (MQ=255)
                           aaTCGTTTACATGGTCAATGTGATACCTGAACCCCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGAttt                                                                                                                    >  3:71585/1‑149 (MQ=255)
                                tttACATTGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCttt                                                                                                               <  5:96574/149‑1 (MQ=255)
                                    cATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCTCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTc                                                                                                           >  5:36130/1‑149 (MQ=255)
                                    cATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTc                                                                                                           >  2:104867/1‑149 (MQ=255)
                                     aTGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTcc                                                                                                          >  1:346685/1‑148 (MQ=255)
                                     aTGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCt                                                                                                          >  5:376979/1‑149 (MQ=255)
                                     aTGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCg                                                                                                                                       >  4:270778/1‑120 (MQ=255)
                                     aTGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCg                                                                                                                                       <  3:270778/120‑1 (MQ=255)
                                      tGGTCAATGTGATACCTGAACACCTTCGGGCCACTCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCt                                                                                                          >  1:51046/1‑149 (MQ=255)
                                        gTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCtgt                                                                                                        <  8:254115/148‑1 (MQ=255)
                                         tCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCATGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCtgttg                                                                                                      >  3:100723/1‑149 (MQ=255)
                                              gtgATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGAttt                                                                                                                                                 <  4:253513/101‑1 (MQ=255)
                                              gtgATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGAttt                                                                                                                                                 >  3:253513/1‑101 (MQ=255)
                                              gtgATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTAGACATGACCTGGATTTCCTTTGGTCTGTTGTTGCt                                                                                                 >  4:242488/1‑149 (MQ=255)
                                                 aTACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCCGTTGTTGCTTg                                                                                               <  2:346685/148‑1 (MQ=255)
                                                  tACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGc                                                                                             <  7:258741/149‑1 (MQ=255)
                                                   aCCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCttt                                                                                                               >  4:328654/1‑130 (MQ=255)
                                                    ccTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGa                                                                                                                       >  1:356737/1‑121 (MQ=255)
                                                    ccTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGa                                                                                                                       <  2:356737/121‑1 (MQ=255)
                                                    ccTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGa                                                                                                                       >  3:114812/1‑121 (MQ=255)
                                                    ccTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGa                                                                                                                       <  4:114812/121‑1 (MQ=255)
                                                     cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctg                                                                                          >  4:279066/1‑149 (MQ=255)
                                                        aaCACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttc                                                                                       <  5:72412/149‑1 (MQ=255)
                                                        aaCACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttc                                                                                       <  7:183678/149‑1 (MQ=255)
                                                        aaCACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgt                                                                                         <  6:376979/147‑1 (MQ=255)
                                                          cacCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcct                                                                                     >  2:85971/1‑149 (MQ=255)
                                                           acCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcctg                                                                                    >  4:14384/1‑149 (MQ=255)
                                                             cTTCGGGCCACCCATAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcctgtt                                                                                  >  7:108579/1‑149 (MQ=255)
                                                              ttCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGttt                                                                                 >  6:303472/1‑149 (MQ=255)
                                                                cGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTtc                                                                                <  3:293654/148‑1 (MQ=255)
                                                                    ccACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCgg                                                                           >  2:350496/1‑149 (MQ=255)
                                                                          cTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcct                                                                                           <  7:113243/128‑1 (MQ=255)
                                                                          cTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcct                                                                                           >  8:113243/1‑128 (MQ=255)
                                                                          cTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCtgt                                                                      <  8:147825/149‑1 (MQ=255)
                                                                          cTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCtgt                                                                      >  4:271553/1‑149 (MQ=255)
                                                                               ttttAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCt                                                                  >  3:188457/1‑148 (MQ=255)
                                                                                tttAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTg                                                                 <  6:36130/148‑1 (MQ=255)
                                                                                   aaCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGccgctagg                                                             >  1:324898/1‑143 (MQ=255)
                                                                                    aCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtg                                                            <  1:104867/149‑1 (MQ=255)
                                                                                      cccATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCtt                                                                                                <  8:292509/111‑1 (MQ=255)
                                                                                      cccATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCtt                                                                                                >  7:292509/1‑111 (MQ=255)
                                                                                       ccATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCATGCGGTTGTCATGACCGGGATTTCGTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggtg                                                         <  8:5483/149‑1 (MQ=255)
                                                                                       ccATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggtg                                                         <  3:76890/149‑1 (MQ=255)
                                                                                       ccATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGACATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCGTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggtg                                                         >  2:104501/1‑149 (MQ=255)
                                                                                           gTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATg                                                     >  2:6790/1‑149 (MQ=255)
                                                                                              tCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTa                                                   <  2:271919/148‑1 (MQ=255)
                                                                                               cccATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTAc                                                  >  8:42752/1‑148 (MQ=255)
                                                                                                 cATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTTCTAt                                               >  7:196474/1‑149 (MQ=255)
                                                                                                            ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                      <  8:108579/147‑1 (MQ=255)
                                                                                                            ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                      >  6:355317/1‑147 (MQ=255)
                                                                                                            ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                    <  1:85971/149‑1 (MQ=255)
                                                                                                            ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                    >  2:343564/1‑149 (MQ=255)
                                                                                                             cTATCACGACGTACTCCTGTACGCTAATGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTTGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                    <  3:242488/148‑1 (MQ=255)
                                                                                                                   cgacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAg                             >  4:269911/1‑149 (MQ=255)
                                                                                                                    gacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGt                            >  2:260770/1‑149 (MQ=255)
                                                                                                                          cTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCTTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGAtgt                       <  2:51046/148‑1 (MQ=255)
                                                                                                                          cTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGAtgtg                      <  3:271553/149‑1 (MQ=255)
                                                                                                                          cTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGAtgtg                      <  1:22597/149‑1 (MQ=255)
                                                                                                                           tCCTGTACGCTAAACCGCCTGTTTTGGGTCAGGCGGGTGTGCTGACCGGGATTTGCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAAGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTgg                     <  8:196474/149‑1 (MQ=255)
                                                                                                                                  cGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttc                                                                                       <  6:254831/76‑1 (MQ=255)
                                                                                                                                  cGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttc                                                                                       >  5:254831/1‑76 (MQ=255)
                                                                                                                                  cGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGc              <  3:269911/149‑1 (MQ=255)
                                                                                                                                  cGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGc              >  5:186054/1‑149 (MQ=255)
                                                                                                                                     tAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTa           >  4:353022/1‑149 (MQ=255)
                                                                                                                                     tAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATATCCTTTGGTCTGTTGTTGCTTGGCCTGTTGCTGTTTCGTCGGTCTGTGCCTGAAAAGGTGCATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCt            >  2:54601/1‑148 (MQ=255)
                                                                                                                                       aaccgcCTGACTCAAATCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggt                                                          >  5:368199/4‑100 (MQ=255)
                                                                                                                                       aaccgcCTGACTCAAATCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggt                                                          <  6:368199/97‑1 (MQ=255)
                                                                                                                                       aaGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTa           <  4:188457/147‑1 (MQ=255)
                                                                                                                                       aaGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTAc          <  5:374625/148‑1 (MQ=255)
                                                                                                                                       aaGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATCGTGGATGTACTATGAGCCTTTTGAGTGGCAGTGATGTGGGCAAAGCCTa           <  4:71585/147‑1 (MQ=255)
                                                                                                                                         ccgcCTGACTCAAATCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACaa        >  7:83501/2‑148 (MQ=255)
                                                                                                                                           gcCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAAc     <  1:150920/149‑1 (MQ=255)
                                                                                                                                           gcCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAAc     <  3:14384/149‑1 (MQ=255)
                                                                                                                                             cTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACg    >  1:345540/1‑148 (MQ=255)
                                                                                                                                             cTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACg    >  5:297600/1‑148 (MQ=255)
                                                                                                                                                aTTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTa  <  1:343564/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCTGCAGGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTGTCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTA  >  NC_002947/1987294‑1987583

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: