Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 709,861 C→G 8.7% A201A (GCC→GCG ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947709,8610CG 8.7% 310.4 / 16.0 103A201A (GCC→GCGribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (53/41);  new base G (3/6);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 107 total aligned reads displayed.

GCCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGG  >  NC_002947/709715‑710010
                                                                                                                                                  |                                                                                                                                                      
gccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGc                                                                                                                                                     >  2:1214358/1‑150 (MQ=255)
gccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGc                                                                                                                                                     >  1:2215675/1‑150 (MQ=255)
gccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGc                                                                                                                                                     >  1:1367190/1‑150 (MQ=255)
 ccAACACGGTCATCCAGGATGGCCTGCGGGTCAGAAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACCCCCGCTGGGCCGCTCGTAGCGCATTACATGCCGGGTGCTGCACGGCCAGCTTCTGGccgtgc                                                                                                                                                    <  2:277081/150‑5 (MQ=255)
      acGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTc                                                                                                                                               >  1:1696083/1‑150 (MQ=255)
       cGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              <  1:829101/150‑3 (MQ=255)
       cGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              <  1:369099/150‑3 (MQ=255)
        ggTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGACGCTGGCGGGCCAGCttct                                                                                                                                             <  2:93548/150‑4 (MQ=255)
        ggTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCgg                                                                                                                                             >  1:1297097/1‑150 (MQ=255)
         gTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttctg                                                                                                                                            <  1:1799826/150‑5 (MQ=255)
          tCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGca                                                                                                                                           >  2:923036/1‑150 (MQ=255)
          tCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGca                                                                                                                                           <  2:35764/150‑1 (MQ=255)
          tCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTACACGGCCAGAAGCTGGCCCGCCAGCTCGGca                                                                                                                                           <  2:1637832/150‑1 (MQ=255)
             tCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCacac                                                                                                                                        <  1:923036/150‑1 (MQ=255)
              ccAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCCTTACGGCCACCAAGCTACACGTCCAGAAGCTGTCCTGCCAGCTCGGCACAcc                                                                                                                                       >  2:1232171/1‑150 (MQ=255)
              ccAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACAcc                                                                                                                                       >  1:1810336/1‑150 (MQ=255)
              ccAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCCCGGCACAcc                                                                                                                                       >  2:1741423/1‑150 (MQ=255)
               cAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCa                                                                                                                                      >  1:1961713/1‑150 (MQ=255)
               cAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCAAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCTGCACACCa                                                                                                                                      >  1:1868575/1‑150 (MQ=255)
                  gATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAAcc                                                                                                                                   <  1:2347248/150‑1 (MQ=255)
                       ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCaa                                                                                                                              >  2:649674/1‑150 (MQ=255)
                              gTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAg                                                                                                                       <  2:861684/150‑1 (MQ=255)
                              gTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAg                                                                                                                       >  2:2119997/1‑150 (MQ=255)
                                         cGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCCCCAACATCCAGCTCAAGcgccg                                                                                                            >  2:1317343/1‑150 (MQ=255)
                                               gcgcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGt                                                                                                      <  2:1810336/150‑1 (MQ=255)
                                                 gcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              <  1:1410770/108‑3 (MQ=255)
                                                 gcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              >  2:1410735/1‑106 (MQ=255)
                                                 gcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              >  2:1410770/1‑106 (MQ=255)
                                                 gcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              <  1:1410735/108‑3 (MQ=255)
                                                 gcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGc                                                                                                    >  2:1788127/1‑150 (MQ=255)
                                                 gcAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCAGTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGc                                                                                                    >  2:213523/1‑150 (MQ=255)
                                                   aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCCGTTGAAGCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCg                                                                                                  >  1:1933357/1‑150 (MQ=255)
                                                        cTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTc                                                                                             >  1:289854/1‑150 (MQ=255)
                                                        cTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTc                                                                                             <  2:2215675/150‑1 (MQ=255)
                                                            tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  2:2100835/150‑1 (MQ=255)
                                                            tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  2:1367190/150‑1 (MQ=255)
                                                            tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  2:1297097/150‑1 (MQ=255)
                                                            tCCGAGGGCGACTTAGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  1:1317343/150‑1 (MQ=255)
                                                                aGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGtat                                                                                     <  2:289854/150‑1 (MQ=255)
                                                                 gggCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGtata                                                                                    >  1:1580765/1‑150 (MQ=255)
                                                                   gCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         >  1:466582/1‑143 (MQ=255)
                                                                   gCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  2:466582/143‑1 (MQ=255)
                                                                       cTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc                                                                                                                <  1:1146382/116‑1 (MQ=255)
                                                                       cTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc                                                                                                                >  2:1146382/1‑116 (MQ=255)
                                                                        ttCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGtata                                                                                    <  2:1424045/143‑1 (MQ=255)
                                                                        ttCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGtata                                                                                    >  1:1424045/1‑143 (MQ=255)
                                                                        aaggAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGcc                                                                             >  1:2269529/4‑150 (MQ=255)
                                                                          cGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAg                                                                           <  1:1605557/150‑1 (MQ=255)
                                                                               tGGCCGAACACCTGCTTTGACGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAAAATCCACCTCAAGCGCCGCCGTGTGCCGCTGTCCGGGGCATACCTGGCCAGCATCg                                                                      >  2:1955285/1‑150 (MQ=255)
                                                                                 gCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGcgccgccgc                                                                                                        >  1:1003825/1‑114 (MQ=255)
                                                                                 gCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGcgccgccgc                                                                                                        <  2:1003825/114‑1 (MQ=255)
                                                                                        acCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                             >  2:2156320/1‑150 (MQ=255)
                                                                                        acCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                             >  1:689619/1‑150 (MQ=255)
                                                                                         ccTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACgg                                                            >  2:643242/1‑150 (MQ=255)
                                                                                         ccTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACgg                                                            >  2:2033614/1‑150 (MQ=255)
                                                                                           tgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATCCCTGGCCCGCATCGTGATCGACGGCa                                                          >  2:1145350/1‑150 (MQ=255)
                                                                                           tgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCa                                                          >  2:1590494/1‑150 (MQ=255)
                                                                                             ctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCaaa                                                        <  1:1214358/150‑1 (MQ=255)
                                                                                             ctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCaaa                                                        <  1:981142/150‑1 (MQ=255)
                                                                                                  ccgccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTc                                                                                             <  1:2118498/108‑1 (MQ=255)
                                                                                                  ccgccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTc                                                                                             >  2:2118498/1‑108 (MQ=255)
                                                                                                  ccgccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                   >  1:1866347/1‑150 (MQ=255)
                                                                                                   cgccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCAGCCAGCTCGTCTCACTAACCGCCAACAACCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTgg                                                  <  1:543056/150‑1 (MQ=255)
                                                                                                     cccGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                             <  1:2034200/137‑1 (MQ=255)
                                                                                                     cccGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                             >  2:2034200/1‑137 (MQ=255)
                                                                                                     cccGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGcc                                                >  1:2035792/1‑150 (MQ=255)
                                                                                                       cGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                             <  2:689619/150‑1 (MQ=255)
                                                                                                       cGTACCGCATTACAGGCCGCGTGCTGAACGGCCATAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                             <  2:1537372/150‑1 (MQ=255)
                                                                                                            cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAAcc                                                                                                                                   <  1:785975/60‑1 (MQ=255)
                                                                                                            cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAAcc                                                                                                                                   >  2:785975/1‑60 (MQ=255)
                                                                                                            cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGcgccgc                                                                                                           >  1:32979/1‑84 (MQ=255)
                                                                                                            cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTc                                        >  2:2369111/1‑150 (MQ=255)
                                                                                                             gCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCg                                       <  2:1696083/150‑1 (MQ=255)
                                                                                                               aTTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGc                                     >  1:1382653/1‑150 (MQ=255)
                                                                                                                    aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                       <  2:1266674/112‑1 (MQ=255)
                                                                                                                    aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                       >  1:1266674/1‑112 (MQ=255)
                                                                                                                    aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGa                                                                   <  2:735560/116‑1 (MQ=255)
                                                                                                                    aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGa                                                                   >  1:735560/1‑116 (MQ=255)
                                                                                                                    aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAatat                                >  2:410660/1‑150 (MQ=255)
                                                                                                                    aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCAGCATCGAATCTTGGCCGGGTGTCGGCAatat                                >  1:1810743/1‑150 (MQ=255)
                                                                                                                        cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                       >  1:2309674/1‑108 (MQ=255)
                                                                                                                        cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                       <  2:2309674/108‑1 (MQ=255)
                                                                                                                        cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                            <  2:2178444/150‑1 (MQ=255)
                                                                                                                        cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                            <  1:1105148/150‑1 (MQ=255)
                                                                                                                              ctgcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgcgc                      >  1:2281044/1‑150 (MQ=255)
                                                                                                                               tgcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGcc                     >  2:897443/1‑150 (MQ=255)
                                                                                                                                  aCGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCa                  >  1:226773/1‑150 (MQ=255)
                                                                                                                                   cGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCAc                 <  2:1868575/150‑1 (MQ=255)
                                                                                                                                       cAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGAGgtgttt                                       <  2:684799/125‑3 (MQ=255)
                                                                                                                                       cAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGAGgtgttt                                       >  1:684799/1‑123 (MQ=255)
                                                                                                                                       cAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGtt             <  1:1114411/150‑1 (MQ=255)
                                                                                                                                         cgaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGc           >  1:1147496/3‑150 (MQ=255)
                                                                                                                                          aaGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGcc          >  2:2079220/1‑150 (MQ=255)
                                                                                                                                            gCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAACGCTAGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCgg        <  2:2269529/150‑1 (MQ=255)
                                                                                                                                               ggCCCGCCAGCTCGGCACACCATCCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAACGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGa     <  2:124578/150‑1 (MQ=255)
                                                                                                                                               ggCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGa     <  2:226773/150‑1 (MQ=255)
                                                                                                                                               ggCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGa     <  2:908538/150‑1 (MQ=255)
                                                                                                                                                 cccGCCAGCTCGGCACACCAACCGACAAAATCCAGCTCAAAAGCCGAAGCGTTCCGGGGTCTGGGGTATACCAGGACAGCATCGTGATCGACGTCAAAGCTTGGCCGGGTGTCGGCAATATCCGCGTGCGCCCCACCGTTGCCGGTGATg   <  1:1741423/150‑1 (MQ=255)
                                                                                                                                                  ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg  >  2:1319161/1‑150 (MQ=255)
                                                                                                                                                  ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGACGTGCGCCCCACCGTTGCCGGTGATgg  <  1:1349421/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GCCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGG  >  NC_002947/709715‑710010

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: