Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 219,777 C→A 6.0% P8428Q (CCA→CAA)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947219,7770CA 6.0% 241.9 / 3.3 77P8428Q (CCA→CAA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (38/33);  new base A (3/2);  total (42/35)
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.86e-01

TACTGGAACTACGTCTACTCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGTCGTCAGCGAGG  >  NC_002947/219629‑219926
                                                                                                                                                    |                                                                                                                                                     
tACTGGAACTACGTCTACTCGGCCGCAGGTACCAGCGACGCAACCGCATGTTTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCca                                                                                                                                                      <  1:1599420/150‑1 (MQ=255)
    ggAACTACGTCTACTCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGa                                                                                                                                                  <  1:1131036/150‑1 (MQ=255)
    ggAACTACGTCTACTCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGa                                                                                                                                                  <  2:1613787/150‑1 (MQ=255)
     gAACTACGTCTACTCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGACGGCACCAATGGCTTGAGCAAAGAAGCCAATACGGGCCGGAACAGCGATGTCACCCTGCACGATCTGAAACCTTCCGACAGTGTCGGCAAGCCACAGAc                                                                                                                                                 >  1:1104086/1‑150 (MQ=255)
            gTCTACTCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATc                                                                                                                                          <  1:382135/150‑1 (MQ=255)
                  tCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATc                                                                                                                                          <  2:972615/144‑1 (MQ=255)
                  tCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATc                                                                                                                                          >  1:972615/1‑144 (MQ=255)
                  tCGGCCGCAGGTTCTACTGAAAGTATCACAGTATTGGGCGCCACCACTTGCTTGACCAAAGTCCAGTCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCGTACGACAGTGCCGGCAAGCCACAGACCACTATCGAACCg                                                                                                                                    >  2:696318/1‑150 (MQ=255)
                    ggCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTc                                                                                                                                  <  2:1148458/150‑1 (MQ=255)
                      ccGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCg                                                                                                                                >  1:122615/1‑150 (MQ=255)
                                tACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGc                                                                                                                         >  1:1321124/1‑147 (MQ=255)
                                tACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGc                                                                                                                         <  2:1321124/147‑1 (MQ=255)
                                  cTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAgg                                                                                                                    <  2:122615/150‑1 (MQ=255)
                                   tGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAatat                                                                                                                                           >  1:1482235/1‑126 (MQ=255)
                                   tGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAatat                                                                                                                                           <  2:1482235/126‑1 (MQ=255)
                                   tGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGc                                                                                                                   <  1:1341077/150‑1 (MQ=255)
                                   tGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTAAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGc                                                                                                                   >  1:1386982/1‑150 (MQ=255)
                                       aGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGAtt                                                                                                               <  2:213413/150‑1 (MQ=255)
                                         tttCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATtct                                                                                                             >  1:389871/1‑150 (MQ=255)
                                                  aTTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTc                                                                                                    <  2:1104086/150‑1 (MQ=255)
                                                      ggcggcACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCTGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCaa                                                                                                                      <  1:1519593/128‑1 (MQ=255)
                                                      ggcggcACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCaa                                                                                                                      >  2:1519593/1‑128 (MQ=255)
                                                      ggcggcACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGACCGATTTGGCCAAGGCGATTCTCGGTCACTCCGaa                                                                                                <  2:415754/150‑1 (MQ=255)
                                                         ggcACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCg                                                                                                           >  2:1411259/1‑136 (MQ=255)
                                                         ggcACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCg                                                                                                           <  1:1411259/136‑1 (MQ=255)
                                                            aCCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAg                                                                                                                     >  2:292577/1‑123 (MQ=255)
                                                            aCCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAg                                                                                                                     <  1:292577/123‑1 (MQ=255)
                                                              cAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGt                                                                                        <  1:1069231/150‑1 (MQ=255)
                                                                     ttGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAgg                                                                                                                    <  1:796716/115‑1 (MQ=255)
                                                                     ttGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAgg                                                                                                                    >  2:796716/1‑115 (MQ=255)
                                                                         gCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGAGACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCg                                                                             >  2:1622505/1‑150 (MQ=255)
                                                                               tCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCAACAGACCAATATCGACTCGTCCGAGTTGGCCAAGGCGATTCTCGCTAACGAGGAGGTAAGGGTGCCAGGTGCTGAGACCa                                                                       >  2:1440354/1‑150 (MQ=255)
                                                                                 cAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCAtt                                                                     >  2:1077858/1‑150 (MQ=255)
                                                                                  aGGCTATCGGCCTGAACAGCGATGTCCCGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGGTTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTg                                                                    >  1:1013854/1‑150 (MQ=255)
                                                                                  aGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCAACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGt                                                                                        <  1:1336852/130‑1 (MQ=255)
                                                                                  aGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCAACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGt                                                                                        <  1:1337056/130‑1 (MQ=255)
                                                                                  aGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCAACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGt                                                                                        >  2:1337056/1‑130 (MQ=255)
                                                                                  aGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCAACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGt                                                                                        >  2:1336852/1‑130 (MQ=255)
                                                                                   ggCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGAttt                                                                                                                            >  1:114427/1‑93 (MQ=255)
                                                                                   ggCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGAttt                                                                                                                            <  2:114427/93‑1 (MQ=255)
                                                                                   ggCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGa                                                                   >  1:1417165/1‑150 (MQ=255)
                                                                                     cTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATa                                                                 >  1:91684/1‑150 (MQ=255)
                                                                                      tATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGt                                                                                                         <  1:686095/109‑1 (MQ=255)
                                                                                      tATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGt                                                                                                         >  2:686095/1‑109 (MQ=255)
                                                                                       aTCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTTGCCAAGGCGATTCTCGGTCACTCCGAAGCCACGGTGCCAGGTGCCGACACCATGCATAgc                                                               >  2:451493/1‑149 (MQ=255)
                                                                                       aTCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGt                                                               <  2:91684/150‑1 (MQ=255)
                                                                                          ggCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCAt                                                                      >  1:1551246/1‑140 (MQ=255)
                                                                                          ggCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCAt                                                                      <  2:1551246/140‑1 (MQ=255)
                                                                                          ggCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGt                                                            >  1:1275302/1‑150 (MQ=255)
                                                                                           gCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTa                                                           >  2:68981/1‑150 (MQ=255)
                                                                                              tGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACg                                                        <  2:1417165/150‑1 (MQ=255)
                                                                                                      gATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGAc                                                >  2:942153/1‑150 (MQ=255)
                                                                                                       aTGTCACGCCGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGAca                                               <  1:696318/150‑1 (MQ=255)
                                                                                                          tCACGCTGAACGATCTAAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGAcatca                                            <  2:1386982/150‑1 (MQ=255)
                                                                                                            aCGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGAca                                               >  2:801512/1‑145 (MQ=255)
                                                                                                            aCGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGAca                                               <  1:801512/145‑1 (MQ=255)
                                                                                                                ggAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGCCCGATTTGGCCAAGGCGATTCTCGGTCCCTCCGAAGCATCGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCg                                      <  1:451493/149‑1 (MQ=255)
                                                                                                                  aaCTATCTTAAACCTTACGACAGTTCCGGCAAGCCACACACCAATATCGCCCCGTCCGATTTGGCCAATGCGATTCTCGGTCACTCCGAATCAACCTTCCCAGGTGACAACACCATTCATAGTGGTAACGGCAACGACATCATCTTCGGt                                    <  1:1440354/150‑1 (MQ=255)
                                                                                                                  aaCGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGt                                    >  2:1544637/1‑150 (MQ=255)
                                                                                                                  aaCGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGt                                    >  1:1310193/1‑150 (MQ=255)
                                                                                                                    cGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACg                                                        >  2:280478/1‑128 (MQ=255)
                                                                                                                    cGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACg                                                        <  1:280478/128‑1 (MQ=255)
                                                                                                                     gATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGAc                                 <  2:1275302/150‑1 (MQ=255)
                                                                                                                     gATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGAc                                 >  1:1451738/1‑150 (MQ=255)
                                                                                                                      aTCTGAAACCTTACGACAGTGCCTGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTACCAGATGCCGACACCATTCAGAGTGGTAACGGCAACGACATCATCTTCGGTGAcc                                >  2:1334630/1‑150 (MQ=255)
                                                                                                                            aaCCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTACGGGGACCTGATCa                          >  1:895507/1‑150 (MQ=255)
                                                                                                                             aCCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGAc                                                                                                                                       <  2:1419144/40‑1 (MQ=255)
                                                                                                                             aCCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGAc                                                                                                                                       >  1:1419144/1‑40 (MQ=255)
                                                                                                                                ttACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTc                                                                                                    <  1:1049567/72‑1 (MQ=255)
                                                                                                                                ttACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTc                                                                                                    >  2:1049567/1‑72 (MQ=255)
                                                                                                                                      cAGTGCAGGCAAGCCACAGACAAATATCGACCCGTGCGATTTTGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACgg                >  2:17695/1‑150 (MQ=255)
                                                                                                                                       aGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGACGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGGGACCTGATCACCCTCAACGgt               >  1:1538782/1‑150 (MQ=255)
                                                                                                                                       aGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGgt               >  1:913529/1‑150 (MQ=255)
                                                                                                                                         tGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATtct                                                                                                             <  1:852031/54‑1 (MQ=255)
                                                                                                                                         tGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATtct                                                                                                             >  2:852031/1‑54 (MQ=255)
                                                                                                                                                aaGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTTGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTgtc            <  2:1410514/144‑1 (MQ=255)
                                                                                                                                                aaGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTTGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTgtc            >  1:1410514/1‑144 (MQ=255)
                                                                                                                                                aaGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGCTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGTCGTCAGc      >  2:470283/1‑150 (MQ=255)
                                                                                                                                                   ccACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGTCGTCAGCGAg   >  1:748935/1‑150 (MQ=255)
                                                                                                                                                    cacaGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGaca                                                                          >  1:969982/1‑78 (MQ=255)
                                                                                                                                                    cacaGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGaca                                                                          <  2:969982/78‑1 (MQ=255)
                                                                                                                                                    cacaGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGTCGTCAGCGAgg  >  1:449718/1‑150 (MQ=255)
                                                                                                                                                    cacaGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGTCGTCAGCGAgg  >  2:979821/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACTGGAACTACGTCTACTCGGCCGCAGGTTCTACTGAAAGTTTCGCAGTATTGGGCGGCACCAATGGCTTGAGCAAAGTCCAGGCTATCGGCCTGAACAGCGATGTCACGCTGAACGATCTGAAACCTTACGACAGTGCCGGCAAGCCACAGACCAATATCGACCCGTCCGATTTGGCCAAGGCGATTCTCGGTCACTCCGAAGCAACGGTGCCAGGTGCCGACACCATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGTCGTCAGCGAGG  >  NC_002947/219629‑219926

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