Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,287,7710CA42.4% 204.0 / inf 316intergenic (‑171/+58)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (91/90);  new base A (79/54);  total (171/145)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 320 total aligned reads displayed.

CTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTA  >  NC_002947/2287622‑2287920
                                                                                                                                                     |                                                                                                                                                      
cTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                                                                                                                                                        <  2:7337951/150‑1 (MQ=255)
    aTGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt                                                                                                                                                    <  2:2274101/150‑1 (MQ=255)
     tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg                                                                                                                                                   <  2:4573764/150‑1 (MQ=255)
      ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc                                                                                                                                                  <  1:5173148/150‑1 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                                                                                                 <  1:68594/150‑1 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                                                                                                 <  1:3866924/150‑1 (MQ=37)
        cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGt                                                                                                                                               <  2:149312/150‑1 (MQ=18)
         gcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGt                                                                                                                                               <  1:7613707/150‑1 (MQ=255)
         gcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGt                                                                                                                                               <  2:5021014/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGGTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                                                                              <  1:1092455/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                                                                              <  2:4227902/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                                                                                              <  2:513390/150‑1 (MQ=255)
           tctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTg                                                                                                                                             <  2:3781663/150‑1 (MQ=255)
                  ggCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg                                                                                                                                      <  1:2839391/150‑1 (MQ=255)
                  ggCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg                                                                                                                                      <  2:3248932/150‑1 (MQ=255)
                   gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTgg                                                                                                                                     <  2:442680/150‑1 (MQ=255)
                         aGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCgg                                                                                                                               <  1:5518580/150‑1 (MQ=255)
                         aGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCgg                                                                                                                               <  2:6031247/150‑1 (MQ=255)
                              tCATTTTAATGGCTACCGTAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCAACAGTCAAACTGCACAGTCTTGTAAGCCTGTGAATTTGGACTGTGGGAGCGGGCGTg                                                                                                                          <  1:2198738/150‑1 (MQ=37)
                                 ttttAATGGCTACCGGAACCCGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCCCAGGCCACCTGCACAGCCTTGGAAGCCTGGGCAGTTGGCCTGTGGGAGCGGGCGTGccc                                                                                                                       <  2:452335/150‑1 (MQ=255)
                                 ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTg                                                                                                                                             <  1:3201938/128‑1 (MQ=255)
                                           tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                                                                                                 <  2:5826350/114‑1 (MQ=25)
                                           tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                                                                                                 >  1:5826350/1‑114 (MQ=25)
                                             ccGGAAACCGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCACCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCg                                                                                                           <  2:2476519/150‑1 (MQ=255)
                                                     aGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGc                                                                                                  >  2:5868945/1‑150 (MQ=255)
                                                          tttGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAAATGCACAGTCTTGGAAGCCGGTGCAGTTGGACAGTGGGAGCGGGCGTGCCCCCGAAGGGGCCGGGGCAGGCGACaa                                                                                             >  2:160535/1‑149 (MQ=25)
                                                             gcgcCCGGGCTGTCAAAATTGTGGGGCCCGCCACCGGCCCCTCCGCGGGCACCCCCCCCCCCACAGCCCACCTGCACCGTGTGGGAAGCCGGGGCAGTTGGACTGTGGGAGGGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAg                                                                                          <  1:160535/150‑1 (MQ=255)
                                                                          gAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGc                                                                                                                                                  >  1:4996890/1‑82 (MQ=17)
                                                                               aTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCCCAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGc                                                                        <  1:1197940/150‑1 (MQ=255)
                                                                                       ccagccaCCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgt                                                                                                                                       <  2:8110942/80‑1 (MQ=255)
                                                                                       ccagccaCCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgt                                                                                                                                       >  1:8110942/1‑80 (MQ=255)
                                                                                                 gCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCGGTGCAGGTGGACGGGGGGTGCGGGCGGGCCCGCGAAGGGGCCGGAGCAGGCGCCAACACTTACCGGGCAGATAGGGCCCGGGCTCACccaccaca                                                      >  2:2744374/1‑148 (MQ=21)
                                                                                                   ccctTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg                                                                                                                                                   >  1:1496116/3‑56 (MQ=255)
                                                                                                    cctTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc                                                                                                                                                  >  2:3656163/2‑56 (MQ=255)
                                                                                                           gggCACGCCCGCTCCCACAGTCCAACTGCACCGGCTTCCACGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGtt                                            >  2:7502284/2‑150 (MQ=21)
                                                                                                             gCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                          <  2:6145323/150‑1 (MQ=255)
                                                                                                                   ccGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctcc                                    >  2:6294607/1‑150 (MQ=255)
                                                                                                                   ccGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctcc                                    >  2:490854/1‑150 (MQ=255)
                                                                                                                    cGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctcct                                   <  2:5718092/150‑1 (MQ=255)
                                                                                                                        cccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCa                               <  2:2461491/150‑1 (MQ=255)
                                                                                                                         ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGGCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGc                                        >  2:5455563/1‑140 (MQ=255)
                                                                                                                         ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg                              >  2:7993460/1‑150 (MQ=255)
                                                                                                                         ccACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAggg                                                                         >  1:7309078/1‑107 (MQ=18)
                                                                                                                          cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                             >  1:7911992/1‑150 (MQ=255)
                                                                                                                          cacaGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                             >  1:4813507/1‑150 (MQ=21)
                                                                                                                           acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                            >  2:259580/1‑150 (MQ=255)
                                                                                                                           acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCCTTCTGCCTCCTGCCAGGc                            >  1:906630/1‑150 (MQ=255)
                                                                                                                           acaGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                            >  2:2527139/1‑150 (MQ=255)
                                                                                                                             aGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCAc                          >  2:3518908/1‑150 (MQ=255)
                                                                                                                               tCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAacac                                                 >  2:685355/1‑125 (MQ=255)
                                                                                                                               tCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa                        >  1:2326665/1‑150 (MQ=255)
                                                                                                                               tCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa                        >  2:2229463/1‑150 (MQ=255)
                                                                                                                               tCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa                        >  1:7982316/1‑150 (MQ=255)
                                                                                                                               tCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa                        >  1:7673997/1‑150 (MQ=255)
                                                                                                                               tCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCa                        >  1:1919469/1‑150 (MQ=255)
                                                                                                                                 cAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                      >  2:6311125/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTg                                                                                                                          >  1:4096139/1‑50 (MQ=255)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacc                                                             >  1:7218510/1‑110 (MQ=21)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacc                                                             <  2:7218510/110‑1 (MQ=21)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                          >  2:72956/1‑129 (MQ=255)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                     >  1:6679840/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                     >  1:4491545/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                     >  2:2999465/1‑150 (MQ=37)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                      >  2:3578842/1‑150 (MQ=18)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                     >  1:6047823/1‑150 (MQ=25)
                                                                                                                                    cTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGc                   >  1:1788226/1‑150 (MQ=255)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGTGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                  >  1:3531642/1‑150 (MQ=255)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCa                               <  2:5303208/137‑1 (MQ=255)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                  >  1:1836297/1‑150 (MQ=255)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                  >  1:7123494/1‑150 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                            >  1:1297844/1‑140 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                  >  2:7552263/1‑150 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                  >  1:6526331/1‑150 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                  >  2:1385438/1‑150 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGtcatca                                                               <  1:6581041/105‑1 (MQ=18)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGtcatca                                                               >  2:6580871/1‑105 (MQ=18)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg                              >  2:4759812/1‑138 (MQ=25)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg                              <  1:4759812/138‑1 (MQ=25)
                                                                                                                                      gCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                               <  2:7960287/57‑1 (MQ=255)
                                                                                                                                      gCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                 >  1:2516039/1‑150 (MQ=255)
                                                                                                                                      gCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCAt                       <  2:1854822/144‑1 (MQ=255)
                                                                                                                                      gCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCAACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGc                                        >  1:530351/1‑128 (MQ=17)
                                                                                                                                       cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacccac                                                          <  1:57147/108‑1 (MQ=255)
                                                                                                                                       cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                      <  2:2812271/144‑1 (MQ=255)
                                                                                                                                       cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc                >  1:5893394/1‑150 (MQ=255)
                                                                                                                                       cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc                >  2:1213041/1‑150 (MQ=255)
                                                                                                                                       cacaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacccacc                                                         <  1:7441561/109‑1 (MQ=25)
                                                                                                                                       cacaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc                >  1:3987854/1‑150 (MQ=255)
                                                                                                                                        acaGTCTTGGACGCCTTTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa               <  1:2804270/150‑1 (MQ=255)
                                                                                                                                        acaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa               >  2:2803753/1‑150 (MQ=255)
                                                                                                                                        acaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa               <  1:6311125/150‑1 (MQ=255)
                                                                                                                                          aGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAgg             >  2:2035458/1‑150 (MQ=255)
                                                                                                                                              ttGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc                <  1:2071402/143‑1 (MQ=255)
                                                                                                                                               tccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCa        <  1:7552263/147‑1 (MQ=255)
                                                                                                                                                ccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAg                                                                           <  2:4350686/80‑1 (MQ=255)
                                                                                                                                                ccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAaca                                                  <  2:6682984/105‑1 (MQ=255)
                                                                                                                                                ccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAg       >  1:6173604/3‑150 (MQ=255)
                                                                                                                                                ccaaGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAg       >  1:6173191/3‑150 (MQ=255)
                                                                                                                                                     ccTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTa  <  1:2005565/150‑1 (MQ=255)
                                                                                                                                                     ccTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTa  >  2:466549/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGTA  >  NC_002947/2287622‑2287920

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