Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,237,790 G→C 5.2% D68H (GAC→CAC)  PP_4614 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,237,7900GC 5.2% 244.4 / 4.7 78D68H (GAC→CAC) PP_4614hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/32);  new base C (2/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

AACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGA  >  NC_002947/5237641‑5237938
                                                                                                                                                     |                                                                                                                                                     
aaCGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCcgacg                                                                                                                                                       >  1:1128348/1‑150 (MQ=255)
      cGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCg                                                                                                                                                 <  1:352124/150‑1 (MQ=255)
       ggTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTACGACGACGCCgg                                                                                                                                                >  2:356287/1‑150 (MQ=255)
          cAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACCCCCAAAACACACGAATCAGTGCTGCATTCGCCCCTGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGt                                                                                                                                             <  1:1887164/150‑1 (MQ=255)
           aaCCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCACCCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTg                                                                                                                                            <  1:450846/150‑1 (MQ=255)
           aaCCTCTGCATCACTCACCCCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTg                                                                                                                                            <  2:1272649/150‑1 (MQ=255)
               tctGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCg                                                                                                                                        >  2:68075/1‑150 (MQ=255)
               tctGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCg                                                                                                                                        >  1:656847/1‑150 (MQ=255)
               tctGCACCGTTCACCGCAGAAGCAAGCCCATGAACAACCCATACGAACGAATCACTTCTGCATCCGCCCCCGAATACCGGGTTAACCTGATCAATGCCAACCTCGACAGCATCGTCATGCTGACCATGGCCGACGACGCCGTCGAGCTCg                                                                                                                                        >  2:816086/1‑150 (MQ=255)
                   cATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCaa                                                                                                                                    <  2:1076653/150‑1 (MQ=255)
                   cATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCaa                                                                                                                                    <  2:1076645/150‑1 (MQ=255)
                    aTCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAg                                                                                                                                   >  2:1577219/1‑150 (MQ=255)
                     tcactcacCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTTGTCGCCAAgc                                                                                                                                  >  2:32006/1‑150 (MQ=255)
                        ctcacCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGcc                                                                                                                               >  1:32112/1‑150 (MQ=255)
                         tcacCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCt                                                                                                                              <  1:414600/150‑1 (MQ=255)
                                    gCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCgg                                                                                                                   <  1:1931054/150‑1 (MQ=255)
                                      aaGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGcc                                                                                                                 >  2:1091834/1‑150 (MQ=255)
                                      aaGACCATGAACAACCCATTCGAACGAAGCAGTGCTGCATTCGCCCCCGAATACCGGATGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCCCCAACCGCCTCCTCACCCCGGcc                                                                                                                 >  2:290000/1‑150 (MQ=255)
                                       aGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCACGAATCCCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCGGGccc                                                                                                                >  1:638305/1‑150 (MQ=255)
                                         aCCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCaa                                                                                                              >  1:1632515/1‑150 (MQ=255)
                                                 caacCCATTCGAACGAGTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAAGCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGa                                                                                                      >  1:1593687/1‑150 (MQ=255)
                                                     ccATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCcgacga                                                                                                                                                      >  1:924447/1‑98 (MQ=255)
                                                     ccATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCcgacga                                                                                                                                                      <  2:924447/98‑1 (MQ=255)
                                                     ccATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCtca                                                                                                                            >  2:1844177/1‑124 (MQ=255)
                                                     ccATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCtca                                                                                                                            <  1:1844177/124‑1 (MQ=255)
                                                         tCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTgc                                                                                              <  2:455020/150‑1 (MQ=255)
                                                            aacgaaTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGCGCATTGCCAACCTCGACGGCAGCATCAGGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGcc                                                                                           >  1:1932415/1‑150 (MQ=255)
                                                            aacgaaTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGcc                                                                                           >  2:942205/1‑150 (MQ=255)
                                                            aacgaaTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACCACGCCGGCGtcgtcggaca                                                                                                                                    <  2:304240/109‑5 (MQ=255)
                                                            aacgaaTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACCACGCCGGCGtcgtcggaca                                                                                                                                    >  1:304240/1‑105 (MQ=255)
                                                             acgaaTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCt                                                                                          <  1:1267741/150‑1 (MQ=255)
                                                               gaaTCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGc                                                                                        >  2:1603502/1‑150 (MQ=255)
                                                                     gTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACggg                                                                                  <  1:942205/150‑1 (MQ=255)
                                                                     gTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACggg                                                                                  >  1:206038/1‑150 (MQ=255)
                                                                        ctgcATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGa                                                                               <  1:22254/150‑1 (MQ=255)
                                                                         tgcACTCGCCCCCGTATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGAt                                                                              <  2:1593687/150‑1 (MQ=255)
                                                                             ttCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGAc                                                                          >  2:194883/1‑150 (MQ=255)
                                                                             ttCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGAc                                                                          >  1:1000110/1‑150 (MQ=255)
                                                                             ttCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGAc                                                                          <  2:32112/150‑1 (MQ=255)
                                                                              tCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACa                                                                         <  1:1905792/150‑1 (MQ=255)
                                                                              tCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACa                                                                         <  2:1000110/150‑1 (MQ=255)
                                                                              tCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACa                                                                         >  2:1905792/1‑150 (MQ=255)
                                                                              tCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTACGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACa                                                                         <  1:356287/150‑1 (MQ=255)
                                                                                gCCCCCGAATACCGGGTGAGCCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGt                                                                                 >  2:443807/1‑140 (MQ=255)
                                                                                gCCCCCGAATACCGGGTGAGCCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGCGt                                                                                 <  1:443807/140‑1 (MQ=255)
                                                                                 cccccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCGCGGCCCAACGCAATGACCCCGTGAGCCTGCAAAGGGAGGTCggaagag                                                                      >  2:415765/1‑144 (MQ=255)
                                                                                  ccccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCtcga                                                                                                                           >  1:941258/1‑94 (MQ=255)
                                                                                  ccccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCtcga                                                                                                                           <  2:941258/96‑3 (MQ=255)
                                                                                  ccccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCAt                                                                     >  1:602714/1‑150 (MQ=255)
                                                                                  ccccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCAt                                                                     >  2:1388485/1‑150 (MQ=255)
                                                                                   cccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGcc                                                                                                                 <  2:1154469/105‑1 (MQ=255)
                                                                                   cccGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGcc                                                                                                                 >  1:1154469/1‑105 (MQ=255)
                                                                                      gAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGAGGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGc                                                                 >  1:830648/1‑150 (MQ=255)
                                                                                      gAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGc                                                                 >  2:284109/1‑150 (MQ=255)
                                                                                          aCCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTgg                                                             <  1:194883/150‑1 (MQ=255)
                                                                                              ggTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCt                                                         >  2:1370779/1‑150 (MQ=255)
                                                                                                    ccTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACTACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCCATGACCCGGCGCGCCTTCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGAt                                                   >  1:699556/1‑150 (MQ=255)
                                                                                                    ccTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGAt                                                   <  2:656847/150‑1 (MQ=255)
                                                                                                    ccTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGAt                                                   >  2:202903/1‑150 (MQ=255)
                                                                                                    ccTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGAt                                                   >  2:25331/1‑150 (MQ=255)
                                                                                                     cTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATc                                                  <  2:1128348/150‑1 (MQ=255)
                                                                                                       ctcaaTTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGa                                               >  2:971797/5‑150 (MQ=255)
                                                                                                           aTTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTc                                           >  1:603308/1‑150 (MQ=255)
                                                                                                            ttGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGcc                                                                                                                               <  2:1653281/66‑1 (MQ=255)
                                                                                                            ttGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGcc                                                                                                                               >  1:1653281/1‑66 (MQ=255)
                                                                                                              gCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGc                                        <  1:1091834/150‑1 (MQ=255)
                                                                                                              gCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGc                                        <  2:1632515/150‑1 (MQ=255)
                                                                                                               ccAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGcc                                       >  2:806688/1‑150 (MQ=255)
                                                                                                                cAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACCACGCCGGCgtcgtc                                                                                                                                         <  2:1013452/52‑1 (MQ=38)
                                                                                                                cAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACCACGCCGGCgtcgtc                                                                                                                                         >  1:1013452/1‑52 (MQ=37)
                                                                                                                  aCCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGa                                    >  1:1350373/1‑150 (MQ=255)
                                                                                                                    cgccaccgcagccTCATTCTGTCCCTGGCCGACGACGCCGTGGGGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATCACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAAc                                  <  2:830648/143‑1 (MQ=255)
                                                                                                                    cTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAAc                                  <  2:602714/150‑1 (MQ=255)
                                                                                                                       gACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAAcccc                               <  2:206038/150‑1 (MQ=255)
                                                                                                                             agcaTCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCaa                         >  2:432957/1‑150 (MQ=255)
                                                                                                                                    tGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTgg                  >  1:1092870/1‑150 (MQ=255)
                                                                                                                                       tGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCcg               >  1:377810/1‑150 (MQ=255)
                                                                                                                                       tGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCcg               >  1:1815218/1‑150 (MQ=255)
                                                                                                                                          cccTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCgcgc            <  1:1370779/150‑1 (MQ=255)
                                                                                                                                             tGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGa         >  1:976492/1‑150 (MQ=255)
                                                                                                                                               tCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGAcc       <  2:976492/150‑1 (MQ=255)
                                                                                                                                                  gacgacGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGt    >  1:402825/1‑150 (MQ=255)
                                                                                                                                                  gacgacGCCGGCGTGGTCGACAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGt    <  1:432957/150‑1 (MQ=255)
                                                                                                                                                   acgacgCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATTTACAGCGTCAGCATGGGGCCCCACAGCTGCACCGTGCAGTTCCCCCTGCAAGTGCTGGCCGCGCTGACCCGTg   >  1:435442/1‑150 (MQ=255)
                                                                                                                                                    cgacgCCGGCGTGGTCGCCAAGCGCCTCATCCCCCCGGCCCAACGCAATGACCCGGTGCGCCGGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGa  <  1:816086/150‑1 (MQ=255)
                                                                                                                                                    cgacgCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGAt                                                   <  2:1489075/102‑1 (MQ=255)
                                                                                                                                                    cgacgCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGAt                                                   >  1:1489075/1‑102 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGTCCGGTCAACCTCTGCATCACTCACCGCAGAAGCAAGACCATGAACAACCCATTCGAACGAATCAGTGCTGCATTCGCCCCCGAATACCGGGTGAACCTGAGCATTGCCAACCTCGACGGCAGCATCATGCTGACCCTGTCCGACGACGCCGGCGTGGTCGCCAAGCGCCTCATCACCCCGGCCCAACGCAATGACCCGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGA  >  NC_002947/5237641‑5237938

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: