Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,130,146 A→C 5.0% intergenic (‑195/+109) gcvH‑I ← / ← PP_0990 glycine cleavage system protein H/DoxX family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,130,1460AC 5.0% 749.1 / 22.0 218intergenic (‑195/+109)gcvH‑I/PP_0990glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base A (88/119);  new base C (8/3);  total (96/122)
Fisher's exact test for biased strand distribution p-value = 6.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 218 total aligned reads displayed.

TCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCC  >  NC_002947/1129999‑1130293
                                                                                                                                                   |                                                                                                                                                   
tctTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc                                                                                                                                                   <  1:3192653/150‑1 (MQ=255)
tctTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATc                                                                                                                                                   >  2:4090788/1‑150 (MQ=255)
   tAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGACTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCaaa                                                                                                                                                >  2:2270195/1‑150 (MQ=255)
   tAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCaaa                                                                                                                                                <  2:78912/150‑1 (MQ=255)
        tAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCtt                                                                                                                                           <  1:4243034/150‑1 (MQ=255)
        tAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCtt                                                                                                                                           <  1:908463/150‑1 (MQ=255)
        tAGGAGTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGATGCTCAAACCCtt                                                                                                                                           <  2:1527880/150‑1 (MQ=255)
          ggATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTtgt                                                                                                                                         <  2:851682/150‑1 (MQ=255)
             ttCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCa                                                                                                                                      <  2:2662999/150‑1 (MQ=255)
             ttCCCCTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCa                                                                                                                                      <  2:495174/150‑1 (MQ=255)
              tCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCAc                                                                                                                                     <  2:2129509/150‑1 (MQ=255)
                ccTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCg                                                                                                                                   >  1:1616911/1‑150 (MQ=255)
                 cTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGa                                                                                                                                  >  2:4364533/1‑150 (MQ=255)
                 cTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATGTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGa                                                                                                                                  >  1:3591673/1‑150 (MQ=255)
                 cTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGACGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACACTTGTGCACAGa                                                                                                                                  >  1:4166348/1‑150 (MQ=255)
                     gCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGc                                                                                                                              >  2:255126/1‑150 (MQ=255)
                      cAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCg                                                                                                                             >  1:4095939/1‑150 (MQ=255)
                      cAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGATACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCGTGTGCACCGACGGCg                                                                                                                             >  2:1071843/1‑150 (MQ=255)
                         taatGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCgagg                                                                                                                          <  2:2029111/150‑1 (MQ=255)
                          catGCGGCCAACCTATCACCTGTCTTCTTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTCTTGTGTGGAATGATCTCGTTACAGGGGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCgagga                                                                                                                         <  2:1321789/149‑1 (MQ=255)
                           atGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCgaggag                                                                                                                        <  1:4005417/150‑1 (MQ=255)
                              cGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCt                                                                                                                     <  1:2911562/150‑1 (MQ=255)
                                gCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCtgt                                                                                                                   <  1:3503116/150‑1 (MQ=255)
                                gCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCtgt                                                                                                                   <  2:753181/150‑1 (MQ=255)
                                  cAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTgtgt                                                                                                                 >  2:2233819/1‑150 (MQ=255)
                                   aaCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTa                                                                                                                <  1:536873/150‑1 (MQ=255)
                                   aaCCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTa                                                                                                                <  1:1192167/150‑1 (MQ=255)
                                     ccTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATc                                                                                                              <  1:5459/150‑1 (MQ=255)
                                       tataAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGt                                                                                                            >  1:709071/1‑150 (MQ=255)
                                        ataAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGtt                                                                                                           >  1:2427838/1‑150 (MQ=255)
                                          aaaaTTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGtttt                                                                                                         <  1:372842/150‑1 (MQ=255)
                                           aaaTTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAATCCGTATTTATTGTGTGGACTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTc                                                                                                        >  2:4257422/1‑150 (MQ=255)
                                            aaTTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCg                                                                                                       <  2:3677134/150‑1 (MQ=255)
                                                   tatttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTacgac                                                                                                >  2:90222/1‑150 (MQ=255)
                                                     tttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGt                                                                                              >  2:2963379/1‑150 (MQ=255)
                                                      ttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTa                                                                                             <  2:2103168/150‑1 (MQ=255)
                                                      ttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTa                                                                                             <  2:1509409/150‑1 (MQ=255)
                                                        atCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGc                                                                                           <  1:3603338/150‑1 (MQ=255)
                                                           aaTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaa                                                                                        <  2:1648168/150‑1 (MQ=255)
                                                            aTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaa                                                                                       <  2:4338892/150‑1 (MQ=255)
                                                            aTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaa                                                                                       >  2:3009390/1‑150 (MQ=255)
                                                             tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                      <  1:2763024/150‑1 (MQ=255)
                                                             tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                      <  1:2786580/150‑1 (MQ=255)
                                                             tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                      <  1:2664389/150‑1 (MQ=255)
                                                             tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                      <  2:3728640/150‑1 (MQ=255)
                                                             tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                      <  1:4507173/150‑1 (MQ=255)
                                                              aGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTATGTATCGTTTTCGCTACGTCGTAGCGAAAAc                                                                                     <  1:774425/150‑1 (MQ=255)
                                                               gATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAAcc                                                                                    >  1:3285920/1‑149 (MQ=255)
                                                                     gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAt                                                                              >  1:2531674/1‑150 (MQ=255)
                                                                     gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAt                                                                              >  2:223614/1‑150 (MQ=255)
                                                                     gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAt                                                                              <  1:318300/150‑1 (MQ=255)
                                                                     gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAt                                                                              <  2:3802466/150‑1 (MQ=255)
                                                                       aaaTTTACGGAGAAGCCGTTTTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAATGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtgt                                                                            <  1:1872222/150‑1 (MQ=255)
                                                                          tttACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATgtggt                                                                         <  1:4003597/150‑1 (MQ=255)
                                                                          tttACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATgtggt                                                                         <  1:3757023/150‑1 (MQ=255)
                                                                               ggAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTg                                                                    >  1:2928390/1‑150 (MQ=255)
                                                                                 agaAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGa                                                                  >  2:2695889/1‑150 (MQ=255)
                                                                                  gaAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGAc                                                                 <  2:2901869/150‑1 (MQ=255)
                                                                                   acgCCGTATTTTTTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAATGGGTTTGAGCATCAAACCGTTGGGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACt                                                                <  1:3911984/148‑1 (MQ=255)
                                                                                     gCCGTATTTATTGTGTGGAATGATATCGTGACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTgg                                                              <  2:2436112/150‑1 (MQ=255)
                                                                                      ccGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                             >  1:922953/1‑150 (MQ=255)
                                                                                        gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                           >  1:1335387/1‑150 (MQ=255)
                                                                                        gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                           <  2:41617/150‑1 (MQ=255)
                                                                                        gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                           >  1:1843112/1‑150 (MQ=255)
                                                                                        gTATTTATTGTGTGGAATGATATCGTTACAGGGGCCTCTGAGAGGAAAAGGGTTTGCGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                           <  1:4247185/150‑1 (MQ=255)
                                                                                         tattttttGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAg                                                          <  1:2869576/150‑1 (MQ=255)
                                                                                         tatttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATAAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAg                                                          <  1:2903459/150‑1 (MQ=255)
                                                                                          atttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAgg                                                         >  2:3640269/1‑150 (MQ=255)
                                                                                           tttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTTGACTGGCGAGGt                                                        >  2:2010135/1‑150 (MQ=255)
                                                                                               ttGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCa                                                    >  1:2357356/1‑150 (MQ=255)
                                                                                                tgtgtgTAATGATATCGTTACCGGTGCCTCTGAGTGGAAAAGGGTTTGCGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCaa                                                   <  1:1915973/150‑1 (MQ=255)
                                                                                                  tgtgTAATTATATCGTTACAGGTGCCTCTGAGAGCAAAAGGGTTTGCGCATAAAACCCTTGTGCACCGACGGCGTGGCGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                 <  1:1538071/150‑1 (MQ=255)
                                                                                                    tgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGc                                               <  2:3866353/150‑1 (MQ=255)
                                                                                                     ggAATGATATCGTTACAGGGGCCTCTGAGAGTAAACGGGTTTGAGCATCAACCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACTTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                              <  1:1172164/150‑1 (MQ=255)
                                                                                                        aTGATATCGTTACAGGTGCCTCTGAGTGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGAGGGCGAGGAGGCTGTGTATCGTTTTCGCTACGAGGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTc                                           <  1:1297313/150‑1 (MQ=255)
                                                                                                         ttatatCGTTACAGGTGCCTCTGAGAGGAAAAGGTTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCg                                          <  2:1887552/148‑1 (MQ=255)
                                                                                                         tGATCTCGTTCCCGGTGCCGCTGAGAGGAAACGGGTTGGAGCAGCAAACCCTTGTGCACCGCCGGCGAGGAGGCGGTGTATCGTTTGCGCTCCGCCGTCGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCg                                          <  1:2296926/150‑1 (MQ=255)
                                                                                                           atatCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGt                                        >  2:633973/1‑150 (MQ=255)
                                                                                                                gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtct                                   >  1:175969/1‑150 (MQ=255)
                                                                                                                  tACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCa                                 <  1:3381200/150‑1 (MQ=255)
                                                                                                                      cgtcccTCTGAGAGCAAAAGGCATTCAGCAACAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGATTTCGCAACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGATCATCTCACttt                             <  2:2140321/146‑1 (MQ=255)
                                                                                                                             ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAga                      >  2:1304159/1‑150 (MQ=255)
                                                                                                                             ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAga                      <  1:3549801/150‑1 (MQ=255)
                                                                                                                              tGAGAGTAAATGGTTTTGAGCATCAAACCCTTGTGCACCGACGGCGTGGAGGCTGTTTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgat                     <  1:2044996/150‑1 (MQ=255)
                                                                                                                                    gAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCa               >  1:307927/1‑150 (MQ=255)
                                                                                                                                     acaaGGGTTTGATGCTCAAACCCTTTTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACATTCAGCAGATGATCCAg              >  1:3820348/3‑150 (MQ=255)
                                                                                                                                     aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACCTGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAg              >  2:3427171/1‑150 (MQ=255)
                                                                                                                                     aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAg              >  2:1752290/1‑150 (MQ=255)
                                                                                                                                     aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTGCAGCTGATGATCCAg              >  2:1351598/1‑150 (MQ=255)
                                                                                                                                       acgggTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGACAACGAGCCAGGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTGTCAGCAGCTGATCCAGCg            >  2:3095522/3‑150 (MQ=255)
                                                                                                                                       aaGGGTTTGATGCTCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCCATGGCCGTCGGTCATCTCACTTTCAGCAGATGATCCcgcg            >  1:3401310/1‑150 (MQ=255)
                                                                                                                                       aaGGGTTTGATGCTCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCg            >  1:3166388/1‑150 (MQ=255)
                                                                                                                                       aaGGGTTTGATGCTCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGGAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCg            >  1:1477904/1‑150 (MQ=255)
                                                                                                                                        aGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGa           <  1:3031494/150‑1 (MQ=255)
                                                                                                                                         gggTTTGATGCTCAAACCCTTTTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGTTCCAGCgta          >  1:3468860/1‑148 (MQ=255)
                                                                                                                                          ggTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGcaa         >  1:3472725/1‑150 (MQ=255)
                                                                                                                                           gTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATATCACTTTCAGCAGATGATCCAGCGaaaa        >  2:1495221/1‑150 (MQ=255)
                                                                                                                                            tttGTGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATTCATGTGGTGGCTGGACTGGCGAGGTGGCAATGGCCTTCGGTCATCTCTCTTTCAGCTGATGTTGCAGCGtaatc       >  2:3219135/1‑148 (MQ=255)
                                                                                                                                                aGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACgcgc   >  2:2833022/1‑150 (MQ=255)
                                                                                                                                                 gCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGcc  <  2:3347918/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCTTAGGGTAGGATTCCCTTAGCAATAATGCGGCCAACCTATAAAATTTCTTATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCC  >  NC_002947/1129999‑1130293

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