Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 895,210 G→C 6.2% intergenic (+47/+17) PP_0777 → / ← PP_0778 glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947895,2100GC 6.2% 927.4 / 37.8 274intergenic (+47/+17)PP_0777/PP_0778glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (132/125);  new base C (7/10);  total (139/135)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 292 total aligned reads displayed.

TTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTC  >  NC_002947/895062‑895359
                                                                                                                                                     |                                                                                                                                                     
tttCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCcgcgc                                                                                                                                                      >  2:1450781/1‑150 (MQ=255)
tttCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCcgcgc                                                                                                                                                      >  2:2419560/1‑150 (MQ=255)
tttCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCcgcgc                                                                                                                                                      >  2:3098887/1‑150 (MQ=255)
tttCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCcgcgc                                                                                                                                                      >  2:4285100/1‑150 (MQ=255)
      gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgggcgg                                                                                                                                                <  1:3834474/150‑1 (MQ=255)
       tggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCGCCTGTGAGATCCAGCGCCGCCCGcgcggc                                                                                                                                               <  1:1599120/150‑1 (MQ=255)
          gCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                            <  2:1054548/150‑1 (MQ=255)
           cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCTTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                           >  1:4710609/1‑150 (MQ=255)
           cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                           >  2:3810988/1‑150 (MQ=255)
            aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                          >  1:1030220/1‑150 (MQ=255)
            aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACtcg                                                                                                                                          <  2:697000/150‑3 (MQ=255)
             aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTGCAGCGCCTGTGAGATCCAGCGCCGCCCGCGCGGCACtcg                                                                                                                                          <  1:5080266/150‑3 (MQ=37)
               gATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtc                                                                                                                                       <  1:5252039/150‑1 (MQ=255)
                  ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                    >  2:1846806/1‑150 (MQ=255)
                  ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                    <  1:135901/150‑1 (MQ=255)
                  ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                    <  2:1032377/150‑1 (MQ=255)
                    ccGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTcaca                                                                                                                                  <  2:1684835/150‑1 (MQ=255)
                     cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAg                                                                                                                                 >  2:3276127/1‑150 (MQ=255)
                       ggTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAggt                                                                                                                               <  2:1184058/150‑2 (MQ=38)
                            tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCTAGTGCCGCGCGGGCGGCGCAGGATCTCACAGGCGCTGc                                                                                                                          >  1:2863357/1‑150 (MQ=255)
                            tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                          <  1:2988348/150‑1 (MQ=255)
                            tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                          >  1:4469214/1‑150 (MQ=255)
                             ggCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCa                                                                                                                         >  2:2023907/1‑150 (MQ=255)
                               cgcgTTTCTCGCCGCGGTCGGCGGCGGATGAGCCAGCGGTGGTGCAGGCCATCGAGATAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                       <  2:3621911/150‑1 (MQ=255)
                               cgcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                       <  2:3357737/150‑1 (MQ=255)
                               cgcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                       <  1:597151/150‑1 (MQ=255)
                                  gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACAt                                                                                                                    >  2:509558/1‑150 (MQ=255)
                                   tttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACAtt                                                                                                                   <  1:2082948/150‑1 (MQ=255)
                                    ttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                                  <  2:1557625/150‑1 (MQ=255)
                                     tctcGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgc                                                                                                                 <  2:3756785/150‑1 (MQ=255)
                                      ctcGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                <  1:2697828/150‑1 (MQ=255)
                                       tcGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcgc                                                                                                               >  1:3411130/1‑150 (MQ=255)
                                        cgccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCa                                                                                                              <  1:1210851/150‑1 (MQ=255)
                                        cgccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCa                                                                                                              <  1:456166/150‑1 (MQ=255)
                                          ccgcGCACGGCGACGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCACACATTGCGCAGa                                                                                                            >  1:4981799/1‑150 (MQ=255)
                                           cgcgCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGaa                                                                                                           >  2:558502/1‑150 (MQ=255)
                                                aCGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGGACGccc                                                                                                      >  2:548535/1‑150 (MQ=255)
                                                         gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGc                                                                                             <  1:6950/150‑1 (MQ=255)
                                                           tGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCAc                                                                                           <  1:3250942/150‑1 (MQ=255)
                                                             aCCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCa                                                                                         <  1:4830090/150‑1 (MQ=255)
                                                              cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                        <  1:1861099/150‑1 (MQ=255)
                                                              cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                        <  1:3670760/150‑1 (MQ=255)
                                                              cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                        <  2:1912114/150‑1 (MQ=255)
                                                               ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAgg                                                                                       >  1:1097145/1‑150 (MQ=255)
                                                               ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAgg                                                                                       >  2:1755746/1‑150 (MQ=255)
                                                               ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAgg                                                                                       >  1:4051190/1‑150 (MQ=255)
                                                                 aGCGGTGGTGCAGGCCATCGAGAAATCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGAACCTGTGAGATAGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGtt                                                                                     >  1:216697/1‑150 (MQ=255)
                                                                     gtggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCg                                                                                 <  2:4432459/150‑1 (MQ=255)
                                                                     gtggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCg                                                                                 <  1:4444197/150‑1 (MQ=255)
                                                                        gtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCa                                                                              >  1:5353231/1‑150 (MQ=255)
                                                                             ggCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTCGCGCCAGGTTGCCGGCGTgcgc                                                                         >  2:4602678/1‑150 (MQ=255)
                                                                               ccATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGGTGCCGGCATGCGCTg                                                                       <  2:1597199/150‑1 (MQ=255)
                                                                                cATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGCTCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGt                                                                      >  1:3929648/1‑150 (MQ=255)
                                                                                 aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                     >  1:3192230/1‑150 (MQ=255)
                                                                                 aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                     <  2:2276385/150‑1 (MQ=255)
                                                                                    gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                                  >  2:883410/1‑150 (MQ=255)
                                                                                    gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                                  <  1:3708798/150‑1 (MQ=255)
                                                                                    gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                                  >  1:191361/1‑150 (MQ=255)
                                                                                    gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                                  >  1:4626894/1‑150 (MQ=255)
                                                                                      gaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGAt                                                                >  2:1796300/1‑150 (MQ=255)
                                                                                       aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATa                                                               >  2:413120/1‑150 (MQ=255)
                                                                                            ccTGGGCTGAAGAGTGCCTGTCGTGATGTTTTGGCACCTGTGAGATCGAGTGCGGCGCGGGCGGCGCTGGATCGCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCTTGCGCTGGTCTGATACCCCt                                                          >  1:2382797/1‑150 (MQ=255)
                                                                                            ccTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCt                                                          >  1:4275188/1‑150 (MQ=255)
                                                                                              tGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCg                                                        >  1:2511782/1‑150 (MQ=255)
                                                                                              tGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCg                                                        >  1:3064221/1‑150 (MQ=255)
                                                                                                gCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGaa                                                      <  2:1019867/150‑1 (MQ=255)
                                                                                                 ccTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                                     <  2:80639/150‑1 (MQ=255)
                                                                                                  cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACt                                                    >  1:2393180/1‑150 (MQ=255)
                                                                                                    gcagagTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGt                                                  <  2:2844199/148‑1 (MQ=255)
                                                                                                     aaGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGtt                                                 <  2:4765144/150‑1 (MQ=255)
                                                                                                            ccTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAAt                                          >  1:1433240/1‑150 (MQ=255)
                                                                                                            ccTGACGTGATGTTTTGGCACCTGAGAGATCGAGAGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAAt                                          >  2:5577436/1‑150 (MQ=255)
                                                                                                              tGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCg                                        >  2:3774027/1‑150 (MQ=255)
                                                                                                              tGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCg                                        >  1:4336976/1‑150 (MQ=255)
                                                                                                                  gTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCgg                                    >  2:319496/1‑150 (MQ=255)
                                                                                                                   tGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGc                                   >  2:3487195/1‑150 (MQ=255)
                                                                                                                             ggCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                          <  2:4359039/150‑1 (MQ=255)
                                                                                                                             agcgccTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGGTGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                          >  1:1137268/5‑150 (MQ=37)
                                                                                                                               cACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcg                        <  2:4913967/150‑1 (MQ=255)
                                                                                                                                  cTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcggat                     <  2:4261824/150‑3 (MQ=255)
                                                                                                                                   tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                    >  1:5623928/1‑150 (MQ=255)
                                                                                                                                   tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                    <  2:4523573/150‑1 (MQ=255)
                                                                                                                                   tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                    >  1:3514393/1‑150 (MQ=255)
                                                                                                                                   tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACAGTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                    <  1:4878338/150‑1 (MQ=255)
                                                                                                                                    gtgAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGGTGCCAAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGc                  >  1:4211873/1‑150 (MQ=25)
                                                                                                                                     tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCGGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGc                  <  2:3740230/150‑1 (MQ=255)
                                                                                                                                     tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCACCGGCCAGCCCTGGCGGCTGGc                  >  1:727016/1‑150 (MQ=255)
                                                                                                                                        gaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATg               >  2:2593761/1‑150 (MQ=255)
                                                                                                                                        gaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATg               >  2:5416807/1‑150 (MQ=255)
                                                                                                                                        gaTCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGGTGCCAAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGc              >  1:3987720/1‑150 (MQ=37)
                                                                                                                                           ccagcgccgccCGCGCGGCACTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTg            >  2:3995653/6‑150 (MQ=255)
                                                                                                                                           ccagcgccgccCGCGCGGCACTCGATCTCACAGGTGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTg            >  2:4709219/6‑150 (MQ=39)
                                                                                                                                             aGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTc          <  1:5177731/150‑1 (MQ=255)
                                                                                                                                              gTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcg         <  1:5462966/150‑1 (MQ=255)
                                                                                                                                               tGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgc        <  2:3859163/150‑1 (MQ=255)
                                                                                                                                                 ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCgcgc      >  2:1326844/1‑150 (MQ=255)
                                                                                                                                                 ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCgcgc      >  2:2749716/1‑150 (MQ=255)
                                                                                                                                                    cgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCAtt   >  2:5345809/1‑150 (MQ=255)
                                                                                                                                                     gcgggcggCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTc  <  2:474452/150‑1 (MQ=255)
                                                                                                                                                     gcgggcggCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTc  <  1:2294033/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTC  >  NC_002947/895062‑895359

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