Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,436,589 C→A 6.1% L161M (CTG→ATG)  gpI → tail formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,436,5890CA 6.1% 209.8 / 2.8 65L161M (CTG→ATG) gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

CACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCGGGGGCTGCACCGGCAGACGATGGAA  >  NC_002947/3436444‑3436732
                                                                                                                                                 |                                                                                                                                               
cACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCtggt                                                                                                                                             <  2:1367901/150‑1 (MQ=255)
 aCGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCtggtg                                                                                                                                            <  2:868099/150‑1 (MQ=255)
 aCGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTCCCCATCATCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCtggtg                                                                                                                                            <  1:798053/150‑1 (MQ=255)
          gCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGa                                                                                                                                   <  1:1176986/150‑1 (MQ=255)
            cAACGATAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGc                                                                                                                                 <  2:440770/150‑1 (MQ=255)
                gagaACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGTTTCGACGACGGCCTGGTGGCTACCGGCGCCAGt                                                                                                                             >  2:78539/1‑150 (MQ=255)
                  gaACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCa                                                                                                                           >  1:185052/1‑150 (MQ=255)
                    aCCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAgg                                                                                                                         <  2:180206/150‑1 (MQ=255)
                     ccggccggGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAggg                                                                                                                        <  2:777694/150‑1 (MQ=255)
                       ggccggGTGAAGGCTCGATCATTTCTCCGCAGCAGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGcc                                                                                                                      >  1:1045843/1‑150 (MQ=255)
                          cggGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGa                                                                                                                                   <  2:1104776/134‑1 (MQ=255)
                          cggGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGa                                                                                                                                   >  1:1104776/1‑134 (MQ=255)
                              tGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAgctgct                                                                                                               <  2:671068/150‑1 (MQ=255)
                                  ggCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCATGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAc                                                                                                           >  2:1180766/1‑150 (MQ=255)
                                   gCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAcc                                                                                                          >  2:186477/1‑150 (MQ=255)
                                   gCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAcc                                                                                                          <  1:859836/150‑1 (MQ=255)
                                   gCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGATGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAcc                                                                                                          >  2:73162/1‑150 (MQ=255)
                                    cTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCAtggc                                                                                                                                             >  1:1055496/1‑113 (MQ=255)
                                    cTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCAtggc                                                                                                                                             <  2:1055496/114‑2 (MQ=255)
                                     tCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGc                                                                                                        >  2:80202/1‑150 (MQ=255)
                                        aTCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAg                                                                                                     >  2:761057/1‑150 (MQ=255)
                                           aTTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACg                                                                                                  >  1:482667/1‑150 (MQ=255)
                                              tctcCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTg                                                                                               <  2:1265861/150‑1 (MQ=255)
                                                          gaggagCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGccc                                                                                   <  1:1337770/150‑1 (MQ=255)
                                                           aggagCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCg                                                                                  >  2:659830/1‑150 (MQ=255)
                                                                  cTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGa                                                                                                                                   >  1:734249/1‑94 (MQ=255)
                                                                  cTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGa                                                                                                                                   <  2:734249/94‑1 (MQ=255)
                                                                   ttGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGtt                                                                          <  2:482667/150‑1 (MQ=255)
                                                                   ctGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGtt                                                                          >  1:1119575/2‑150 (MQ=255)
                                                                       gcgcCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGAcc                                                                      >  1:908648/1‑150 (MQ=255)
                                                                             tGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGt                                                                                                <  1:356373/118‑1 (MQ=255)
                                                                             tGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGt                                                                                                >  2:356373/1‑118 (MQ=255)
                                                                             tGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCAGTTCCGGTTGACCCCTCgg                                                                >  2:259720/1‑150 (MQ=255)
                                                                                 cGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCa                                                            <  1:80940/150‑1 (MQ=255)
                                                                                    tGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGc                                                         >  2:941329/1‑150 (MQ=255)
                                                                                       gacgaAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCt                                                      >  2:1067357/1‑150 (MQ=255)
                                                                                       gacgaAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCGGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGTGCAGccccc                                                      >  2:836924/1‑149 (MQ=255)
                                                                                        acgaAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTg                                                     <  2:1119575/150‑1 (MQ=255)
                                                                                         cgaAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGc                                                    <  2:185052/150‑1 (MQ=255)
                                                                                          gaagaaCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCa                                                   <  1:80202/150‑1 (MQ=255)
                                                                                          gaagaaCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCa                                                   >  2:418667/1‑150 (MQ=255)
                                                                                          gaagaaCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCa                                                   >  2:71141/1‑150 (MQ=255)
                                                                                              aaCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCATGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTc                                               <  1:1180766/150‑1 (MQ=255)
                                                                                               aCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGAc                                                                       <  2:1261254/125‑1 (MQ=255)
                                                                                               aCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGAc                                                                       >  1:1261254/1‑125 (MQ=255)
                                                                                                  aGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCa                                           >  2:310942/1‑150 (MQ=255)
                                                                                                       aGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGcc                                      >  1:314853/1‑150 (MQ=255)
                                                                                                          cACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCAcc                                   <  1:1004318/150‑1 (MQ=255)
                                                                                                            cTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCa                                 >  2:990446/1‑150 (MQ=255)
                                                                                                            cTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCa                                 >  1:610225/1‑150 (MQ=255)
                                                                                                                  cTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGGGGCAGCAGCGCGCGCCCGTTCCGGTTGACCCCTCGAAGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGcgcg                           >  1:1314534/1‑150 (MQ=255)
                                                                                                                   ttCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCgg                          >  2:675037/1‑150 (MQ=255)
                                                                                                                   ttCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCgg                          >  1:385334/1‑150 (MQ=255)
                                                                                                                   ttCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCgg                          <  1:761057/150‑1 (MQ=255)
                                                                                                                      cGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGCCCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCgtggg                       <  2:524696/150‑1 (MQ=255)
                                                                                                                       ggCTCGATGCCCGGTTCGACTACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCGGGGGc                      <  1:73162/150‑1 (MQ=255)
                                                                                                                            ggTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACg                                                                                                  <  1:190181/69‑1 (MQ=255)
                                                                                                                            ggTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACg                                                                                                  >  2:190181/1‑69 (MQ=255)
                                                                                                                                 ccGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCGGGGGCTGCAccatcg            <  1:675037/150‑5 (MQ=255)
                                                                                                                                    cggtctcgaccGCCTGCTTTCTACCTGATCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCGGGGGCTGCACCGGCAGAc         <  2:660620/144‑1 (MQ=255)
                                                                                                                                         acgacgGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAAc                                         <  1:269263/113‑1 (MQ=255)
                                                                                                                                         acgacgGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAAc                                         >  2:269263/1‑113 (MQ=255)
                                                                                                                                           gacgGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCGGGGGCTGCACCGGCAGACGagatcg  >  1:660620/1‑145 (MQ=255)
                                                                                                                                             cgGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGcc                                            >  2:566450/1‑106 (MQ=255)
                                                                                                                                             cgGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGcc                                            <  1:566450/106‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTTGACCCCTCGGCGCAGGCCCTGCAGTTCGCCAACGCCACCCAGGCGCGGGGGCTGCACCGGCAGACGATGGAA  >  NC_002947/3436444‑3436732

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