Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 213,354 T→A 13.0% L6287Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947213,3540TA13.0% 183.6 / 19.7 69L6287Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base A (4/5);  total (34/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 = 3.76e-01

TCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGG  >  NC_002947/213207‑213501
                                                                                                                                                   |                                                                                                                                                   
ccggctcGTCGGTCACCGACGTGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGACTCACGACCTgg                                                                                                                                                   >  1:1197021/7‑150 (MQ=17)
      cGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGa                                                                                                                                             <  2:1329814/150‑1 (MQ=255)
       gTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGAc                                                                                                                                            >  1:1092614/1‑150 (MQ=255)
       gTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGag                                                                                                                                            <  2:909087/150‑2 (MQ=255)
       gTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGag                                                                                                                                            <  1:708369/150‑2 (MQ=255)
          ggTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCGAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                         <  2:1082305/150‑1 (MQ=255)
            tGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCCAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                       <  2:125905/150‑1 (MQ=255)
             gACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                      <  2:207306/150‑1 (MQ=255)
                cgacgaGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                   >  2:701355/1‑150 (MQ=255)
                       ccGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGcc                                                                                                                            <  2:212747/150‑1 (MQ=255)
                          ggTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGag                                                                                                                                            <  1:663959/131‑2 (MQ=255)
                          ggTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGag                                                                                                                                            >  2:663959/1‑130 (MQ=255)
                               cccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                               <  2:18950/139‑1 (MQ=255)
                               cccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                               >  1:18950/1‑139 (MQ=255)
                               cccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGtcat                                                                                                                    >  2:737300/1‑150 (MQ=255)
                                ccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTgg                                                                                                                                                   <  1:731740/118‑1 (MQ=255)
                                ccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTgg                                                                                                                                                   >  2:731740/1‑118 (MQ=255)
                                ccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGtcatc                                                                                                                   <  2:1092614/150‑1 (MQ=255)
                                  aGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGtcatcat                                                                                                                 >  1:86762/1‑150 (MQ=255)
                                   ggCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCTGAGCAACAATGCCCAGGtcatcatca                                                                                                               <  2:416410/150‑1 (MQ=255)
                                      aaCGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAg                                                                                                             >  1:909492/1‑150 (MQ=255)
                                      aaCGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGACTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAg                                                                                                             <  2:1197021/150‑1 (MQ=255)
                                       aCGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAgg                                                                                                            <  1:852333/150‑1 (MQ=255)
                                          aaGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTg                                                                                                         <  2:909492/150‑1 (MQ=255)
                                           aGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTgg                                                                                                        <  1:73016/150‑1 (MQ=255)
                                            ggCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGc                                                                                                       <  1:465994/150‑1 (MQ=255)
                                                 cTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAAc                                                                                                  <  1:701355/150‑1 (MQ=255)
                                                     ttCACGGTCACCATCCCGGCCGACCAGACCTCGGTGGCCAAGAACGTCAAGCCGAGGTTCACCGGGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCACGAATGCCCAGGTCATCATGAAGGCTGGCGAAACCAgc                                                                                              >  2:322584/2‑150 (MQ=17)
                                                     gTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGa                                                                                                                                             <  2:128546/103‑1 (MQ=255)
                                                     gTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGa                                                                                                                                             >  1:128546/1‑103 (MQ=255)
                                                         aGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCGTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGCCTGGCGAAACCAgcgcgc                                                                                          <  2:957971/150‑1 (MQ=255)
                                                          ggTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                         >  2:43124/1‑102 (MQ=255)
                                                          ggTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                         <  1:43124/102‑1 (MQ=255)
                                                                cATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATacac                                                                                   >  1:162372/1‑150 (MQ=255)
                                                                  tCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACAccc                                                                                 >  1:1126148/1‑150 (MQ=255)
                                                                         cGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGGCCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTgc                                                                          <  1:737300/150‑1 (MQ=255)
                                                                         cGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTgc                                                                          <  2:86762/150‑1 (MQ=255)
                                                                                         gCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGa                                                                                                                                             >  2:222275/1‑67 (MQ=255)
                                                                                         gCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGa                                                                                                                                             <  1:222275/67‑1 (MQ=255)
                                                                                           cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGt                                                                                                                                               <  1:311457/63‑1 (MQ=255)
                                                                                           cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGt                                                                                                                                               >  2:311457/1‑63 (MQ=255)
                                                                                            gagaACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTAc                                                       >  1:1183826/1‑150 (MQ=255)
                                                                                             agaACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACa                                                      >  1:913566/1‑150 (MQ=255)
                                                                                              gaACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACaa                                                     <  1:719454/150‑1 (MQ=255)
                                                                                                  gTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAg                                                                                                             <  2:1288350/90‑1 (MQ=255)
                                                                                                  gTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAg                                                                                                             >  1:1288350/1‑90 (MQ=255)
                                                                                                   tcggcttgACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                >  2:1010294/8‑150 (MQ=255)
                                                                                                   tCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                >  1:1317310/1‑150 (MQ=255)
                                                                                                   tCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                >  2:1074925/1‑150 (MQ=255)
                                                                                                   tCAAGCCCACGTTCACCGTGCACATCAACCCCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                <  1:451678/150‑1 (MQ=255)
                                                                                                        ttgACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGc                                           <  1:1010294/148‑1 (MQ=255)
                                                                                                         cGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGcc                                          >  2:388492/1‑150 (MQ=255)
                                                                                                             gTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGAt                                      >  1:7192/1‑150 (MQ=255)
                                                                                                                cACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGc                                                                                                           >  2:1414531/1‑78 (MQ=255)
                                                                                                                cACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGc                                                                                                           <  1:1414531/78‑1 (MQ=255)
                                                                                                                 aCCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                   >  1:911560/1‑117 (MQ=255)
                                                                                                                 aCCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                   <  2:911560/117‑1 (MQ=255)
                                                                                                                   cGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                >  2:683711/1‑134 (MQ=255)
                                                                                                                   cGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                <  1:683711/134‑1 (MQ=255)
                                                                                                                   cGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCt                                              <  2:1272533/136‑1 (MQ=255)
                                                                                                                   cGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCt                                              >  1:1272533/1‑136 (MQ=255)
                                                                                                                        acaTCAACACCGCCTTGGCTCACGACCTAGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCa                           >  1:398927/1‑150 (MQ=255)
                                                                                                                         caTCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCAt                          >  2:947308/1‑150 (MQ=255)
                                                                                                                              acacCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCATCCTGGGCATTAAtt                     >  1:327424/1‑150 (MQ=255)
                                                                                                                                 ccGCCCTGGCCCTCGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCgg                  <  2:7192/150‑1 (MQ=255)
                                                                                                                                    ccTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAg               <  2:1081673/150‑1 (MQ=255)
                                                                                                                                    ccTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAg               >  2:100842/1‑150 (MQ=255)
                                                                                                                                        gggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGa           >  2:1104778/4‑150 (MQ=255)
                                                                                                                                                 ccTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTgg  >  1:423927/1‑150 (MQ=255)
                                                                                                                                                   tGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGAt                                                 <  2:303409/101‑1 (MQ=255)
                                                                                                                                                   tGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGAt                                                 >  1:303409/1‑101 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGG  >  NC_002947/213207‑213501

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: