Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,101,800 A→C 12.5% A127A (GCA→GCC) ‡ mlaD → phospholipid ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,101,8000AC12.5% 181.6 / 16.7 71A127A (GCA→GCC) ‡mlaDphospholipid ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base C (6/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GTCGGGTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTG  >  NC_002947/1101658‑1101944
                                                                                                                                              |                                                                                                                                                
gTCGGGTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatc                                                                                                                                            <  4:213946/149‑1 (MQ=255)
     gTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGa                                                                                                                                       <  5:121932/149‑1 (MQ=255)
     gTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGa                                                                                                                                       >  2:84219/1‑149 (MQ=255)
      tGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatcc                                                                                                                                           <  5:154117/144‑1 (MQ=255)
      tGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatcc                                                                                                                                           >  6:355290/1‑144 (MQ=255)
      tGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatcc                                                                                                                                           <  5:355290/144‑1 (MQ=255)
      tGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatcc                                                                                                                                           >  6:154117/1‑144 (MQ=255)
              cAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGt                                                                                                                               <  4:240288/148‑1 (MQ=255)
                gCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatccat                                                                                                                                         <  8:14989/136‑1 (MQ=255)
                gCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAAccatccat                                                                                                                                         >  7:14989/1‑136 (MQ=255)
                 cTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTTGCAccatcctt                                                                                                                                         <  8:211175/135‑3 (MQ=255)
                 cTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTTGCAccatcctt                                                                                                                                         >  7:211175/1‑133 (MQ=255)
                    gACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACt                                                                                                                        >  5:170795/1‑149 (MQ=255)
                     aCAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTg                                                                                                                       <  8:64580/149‑1 (MQ=255)
                      cAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTgg                                                                                                                      >  5:235109/1‑149 (MQ=255)
                         gtcggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGc                                                                                                                   >  6:25164/1‑149 (MQ=255)
                         gtcggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGc                                                                                                                   >  6:367132/1‑149 (MQ=255)
                         gtcggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGc                                                                                                                   >  7:245846/1‑149 (MQ=255)
                         gtcggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGc                                                                                                                   >  8:274797/1‑149 (MQ=255)
                         gtcggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTGGCTTGGTGAGAAGTACATAGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCCTACATGACACACAGACGGCCCTGGTGc                                                                                                                   >  4:239019/1‑149 (MQ=255)
                                     ccTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGa                                                                                                       >  2:355835/1‑149 (MQ=255)
                                      cTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAt                                                                                                      <  4:168626/149‑1 (MQ=255)
                                         ccgaccgaCTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTg                                                                                                    <  2:285055/148‑1 (MQ=255)
                                              cgacTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGt                                                                                              >  8:192749/1‑149 (MQ=255)
                                                    cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTg                                                                                                                 <  7:90406/124‑1 (MQ=255)
                                                    cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTg                                                                                                                 >  8:90406/1‑124 (MQ=255)
                                                    cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGc                                                                                        >  2:152450/1‑149 (MQ=255)
                                                    cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGc                                                                                        <  6:170795/149‑1 (MQ=255)
                                                    cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGc                                                                                        >  2:288035/1‑149 (MQ=255)
                                                     aCTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCt                                                                                       <  7:274797/149‑1 (MQ=255)
                                                      cTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTc                                                                                      <  8:245846/149‑1 (MQ=255)
                                                       tGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTc                                                                                      >  4:282745/1‑148 (MQ=255)
                                                         ccTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCg                                                                                                                             <  6:297568/107‑1 (MQ=255)
                                                         ccTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCg                                                                                                                             >  5:297568/1‑107 (MQ=255)
                                                             gATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAActc                                                                                 <  4:323121/147‑1 (MQ=255)
                                                                 cTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAtt                                                                                                     >  1:138550/1‑123 (MQ=255)
                                                                 cTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAtt                                                                                                     <  2:138550/123‑1 (MQ=255)
                                                                 cTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTgg                                                                           >  1:320780/1‑149 (MQ=255)
                                                                    aCCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCaa                                                                        <  5:25164/149‑1 (MQ=255)
                                                                     ccgccgGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCaa                                                                        <  2:320780/148‑1 (MQ=255)
                                                                     ccgccgGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAg                                                                       <  1:152450/149‑1 (MQ=255)
                                                                      cgccgGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAg                                                                       >  8:304974/1‑148 (MQ=255)
                                                                                 ttGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAAc                                                                   <  7:297625/141‑1 (MQ=255)
                                                                                 ttGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAAc                                                                   >  8:297625/1‑141 (MQ=255)
                                                                                   ggTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaag                                                          <  7:192749/148‑1 (MQ=255)
                                                                                   ggTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                           <  6:235109/147‑1 (MQ=255)
                                                                                   ggTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                           <  5:367132/147‑1 (MQ=255)
                                                                                          aGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTCGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCg                                                  >  1:150567/1‑149 (MQ=255)
                                                                                           gTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCgg                                                 <  1:355835/149‑1 (MQ=255)
                                                                                                  ggCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAAtt                                           <  1:288035/148‑1 (MQ=255)
                                                                                                    cATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAg                                        >  8:271932/1‑149 (MQ=255)
                                                                                                     aTCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAg                                        >  3:223417/1‑148 (MQ=255)
                                                                                                                   tGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGAttt                         <  7:271932/149‑1 (MQ=255)
                                                                                                                       gACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGAt                     <  3:282745/149‑1 (MQ=255)
                                                                                                                       gACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGAt                     >  7:206331/1‑149 (MQ=255)
                                                                                                                        aCGAAGTGCTCAAGGATGGTGCAACCATCCATGCCACACAGACGGCACTGGTGCTCGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGGTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCAACATATCGaac                    >  5:355868/1‑147 (MQ=255)
                                                                                                                        aCGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGaa                                                                    >  1:310048/1‑101 (MQ=255)
                                                                                                                        aCGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGaa                                                                    <  2:310048/101‑1 (MQ=255)
                                                                                                                          gAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCt                  >  4:207800/1‑149 (MQ=255)
                                                                                                                              gtgtCATGGATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGAttt                         >  2:191203/4‑138 (MQ=255)
                                                                                                                              gtgtCATGGATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGAttt                         <  1:191203/135‑1 (MQ=255)
                                                                                                                              gtgtCATGGATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGa              >  7:87444/4‑149 (MQ=255)
                                                                                                                                 tCATGGATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                           <  4:118691/101‑1 (MQ=255)
                                                                                                                                 tCATGGATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                           >  3:118691/1‑101 (MQ=255)
                                                                                                                                 tCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATc                    <  6:250150/140‑1 (MQ=255)
                                                                                                                                 tCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATc                    >  5:250150/1‑140 (MQ=255)
                                                                                                                                   aaGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTg          <  4:223417/148‑1 (MQ=255)
                                                                                                                                    tggATGGTTGCACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc        >  1:47331/2‑149 (MQ=255)
                                                                                                                                    tggATGGTTGCACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc        >  5:153388/2‑149 (MQ=255)
                                                                                                                                          gTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCctgctg  >  6:266535/1‑149 (MQ=255)
                                                                                                                                            gCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGc                                                >  2:135799/1‑101 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GTCGGGTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTG  >  NC_002947/1101658‑1101944

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