Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,297,196 A→C 6.4% intergenic (‑71/+23) yjcO ← / ← fdhF Sel1 repeat‑containing protein YjcO/formate dehydrogenase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,297,1960AC 6.4% 583.7 / 26.2 189intergenic (‑71/+23)yjcO/fdhFSel1 repeat‑containing protein YjcO/formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base A (86/91);  new base C (7/5);  total (93/96)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 198 total aligned reads displayed.

ATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGT  >  NC_000913/4297048‑4297344
                                                                                                                                                    |                                                                                                                                                    
   cTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggag                                                                                                                                                  <  1:503620/150‑1 (MQ=255)
     ggCTGCCATGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                                <  1:205156/150‑1 (MQ=255)
     ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                                >  2:909517/1‑150 (MQ=255)
     ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                                <  2:1199287/150‑1 (MQ=255)
      gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCt                                                                                                                                               >  2:1034807/1‑150 (MQ=255)
       cTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTg                                                                                                                                              >  2:73763/1‑150 (MQ=255)
        tGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGt                                                                                                                                             <  1:291697/150‑1 (MQ=255)
         gCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTa                                                                                                                                            <  1:666866/150‑1 (MQ=255)
            aGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGttttt                                                                                                                                         <  2:31744/150‑5 (MQ=255)
             ggCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGaaa                                                                                                                                        >  2:1422361/1‑150 (MQ=255)
             ggCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGaaa                                                                                                                                        <  2:1234623/150‑1 (MQ=255)
             ggCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGaaa                                                                                                                                        >  2:1010036/1‑150 (MQ=255)
             ggCTCGGAGTCGTTGGCGTGGACAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGa                                                                                                                                          >  1:1005160/1‑148 (MQ=255)
             ggCTCGGAGTCGTTGGCGTGGACAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGa                                                                                                                                          <  2:1005160/148‑1 (MQ=255)
                 cGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGAc                                                                                                                                    <  2:1051309/150‑1 (MQ=255)
                  ggAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACg                                                                                                                                   <  1:443739/150‑1 (MQ=255)
                    aGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                                                                                 <  1:1012439/150‑1 (MQ=255)
                     gTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTa                                                                                                                                >  1:1112197/1‑150 (MQ=255)
                       cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                              >  1:225995/1‑150 (MQ=255)
                       cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                              >  2:86437/1‑150 (MQ=255)
                       cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                              <  2:775334/150‑1 (MQ=255)
                        gTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTa                                                                                                                             >  1:1251061/1‑150 (MQ=255)
                         ttGGCGTGTGCAAAGAATGTCAAAAACAACATCAATGCTATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATTGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTAc                                                                                                                            <  2:538667/150‑1 (MQ=255)
                          tGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACg                                                                                                                           <  1:1234467/150‑1 (MQ=255)
                          tGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACg                                                                                                                           <  2:70122/150‑1 (MQ=255)
                               tGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGt                                                                                                                      >  2:1075538/1‑150 (MQ=255)
                               tGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGt                                                                                                                      <  1:909517/150‑1 (MQ=255)
                                 ggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGc                                                                                                                    >  1:1383270/1‑150 (MQ=255)
                                  gCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGcc                                                                                                                   >  2:282125/1‑150 (MQ=255)
                                   cAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCAACTGccc                                                                                                                  >  1:429249/1‑149 (MQ=255)
                                    aaaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGc                                                                                                                 >  2:1252619/1‑150 (MQ=255)
                                    aaaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGc                                                                                                                 <  1:961300/150‑1 (MQ=255)
                                    aaaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGc                                                                                                                 >  2:1247760/1‑150 (MQ=255)
                                     aaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCt                                                                                                                >  1:161275/1‑150 (MQ=255)
                                          tGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTcgcg                                                                                                           <  2:782116/150‑1 (MQ=255)
                                              aaaaaCAACATCAATGCGATAATTTTTTTCATTTTGATAAAACGATTGGGTAGGGAGTAACAAGTATAAAGAGTTAGAATGGATAAAAAAACAGACTCCGAAAGGAGGCTTTAGAAAGGACGGAATTACGACAGTGCCGCTACGCGCAgg                                                                                                       >  1:1385499/1‑150 (MQ=255)
                                               aaaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGc                                                                                                      >  1:61331/1‑150 (MQ=255)
                                                     acaTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGttt                                                                                                >  2:753748/1‑150 (MQ=255)
                                                         cAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCaa                                                                                            <  2:1300257/150‑1 (MQ=255)
                                                           aTGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACt                                                                                          <  2:1209430/150‑1 (MQ=255)
                                                           aTGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACt                                                                                          <  2:225995/150‑1 (MQ=255)
                                                            tGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACtt                                                                                         <  2:139645/150‑1 (MQ=255)
                                                              cGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGAGGGTATTACGCCAGTGCTGCTTCGCGCAGGCGAGTTTTCAACttgt                                                                                       >  2:522235/1‑150 (MQ=255)
                                                              cGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgt                                                                                       <  2:161275/150‑1 (MQ=255)
                                                              cGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgt                                                                                       <  1:749215/150‑1 (MQ=255)
                                                               gATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgtt                                                                                      >  1:82475/1‑150 (MQ=255)
                                                                    tttttttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCTCAGGCGAGTTTTCAACTTGTTGTACt                                                                                 <  1:336868/150‑1 (MQ=255)
                                                                     ttttttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTc                                                                                >  2:465095/1‑150 (MQ=255)
                                                                      tttttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCg                                                                               <  1:777071/150‑1 (MQ=255)
                                                                       ttttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGt                                                                              >  1:899239/1‑150 (MQ=255)
                                                                        tttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACg                                                                                                                                   >  2:892104/1‑96 (MQ=255)
                                                                        tttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTc                                                                             >  1:802910/1‑150 (MQ=255)
                                                                         ttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCa                                                                            <  2:1060925/150‑1 (MQ=255)
                                                                             ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACg                                                                                                                                   <  2:838575/91‑1 (MQ=255)
                                                                             ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACg                                                                                                                                   >  1:838575/1‑91 (MQ=255)
                                                                             ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCa                                                                        <  2:1383270/150‑1 (MQ=255)
                                                                               ttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCg                                                                               <  1:1101891/141‑1 (MQ=255)
                                                                               ttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCg                                                                               >  2:1101891/1‑141 (MQ=255)
                                                                                 gATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTcg                                                                                                             >  2:297244/1‑109 (MQ=255)
                                                                                   taataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACt                                                                  >  2:1337039/1‑150 (MQ=255)
                                                                                   taataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACt                                                                  >  2:490872/1‑150 (MQ=255)
                                                                                     ataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGc                                                                >  1:99870/1‑150 (MQ=255)
                                                                                        cGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGt                                                                                                                                             >  2:968656/1‑70 (MQ=255)
                                                                                        cGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGt                                                                                                                                             <  1:968656/70‑1 (MQ=255)
                                                                                           ttGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcggcg                                                          >  1:1363494/1‑150 (MQ=255)
                                                                                             gggTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCggcggc                                                        >  2:1313176/1‑150 (MQ=255)
                                                                                                tAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGgcgcg                                                     <  1:1163849/150‑1 (MQ=255)
                                                                                                        aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcgg                                             <  1:1059485/150‑1 (MQ=255)
                                                                                                         aCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggc                                            >  2:1186720/1‑150 (MQ=255)
                                                                                                           cAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcga                                          >  1:724148/1‑150 (MQ=255)
                                                                                                               ataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcgg                                      >  2:256845/1‑150 (MQ=255)
                                                                                                                  aaGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTc                                   <  2:966403/150‑1 (MQ=255)
                                                                                                                  aaGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTc                                   >  2:1347657/1‑150 (MQ=255)
                                                                                                                  aaGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTc                                   >  2:1352611/1‑150 (MQ=255)
                                                                                                                    gagTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGa                                 <  1:282125/150‑1 (MQ=255)
                                                                                                                     agTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAc                                <  2:424864/150‑1 (MQ=255)
                                                                                                                     agTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAc                                <  1:1349086/150‑1 (MQ=255)
                                                                                                                       ttCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAcgc                              <  1:314172/150‑1 (MQ=255)
                                                                                                                         cGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGa                            >  1:378119/1‑150 (MQ=255)
                                                                                                                            aTGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACg                         <  1:1304823/150‑1 (MQ=255)
                                                                                                                            aTGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACg                         <  2:1363494/150‑1 (MQ=255)
                                                                                                                             tGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACgg                        >  1:1078007/1‑150 (MQ=255)
                                                                                                                             tGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACgg                        >  1:1284587/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcgc                     <  1:236262/150‑1 (MQ=255)
                                                                                                                                  ctttctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAg                   <  1:535902/144‑1 (MQ=255)
                                                                                                                                  aaaaaaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAg                   <  2:1112197/150‑1 (MQ=255)
                                                                                                                                      ctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAtt               >  1:521067/3‑150 (MQ=255)
                                                                                                                                      ctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAtt               >  2:677033/3‑150 (MQ=255)
                                                                                                                                       taCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAttt              >  1:1418203/2‑150 (MQ=255)
                                                                                                                                       taCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAttt              >  2:577222/2‑150 (MQ=255)
                                                                                                                                        aCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTg             <  2:1124657/150‑1 (MQ=255)
                                                                                                                                         cAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGt            <  1:905110/150‑1 (MQ=255)
                                                                                                                                         cAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGt            <  2:1024457/150‑1 (MQ=255)
                                                                                                                                          aGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTa           >  1:558495/1‑150 (MQ=255)
                                                                                                                                           gCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTAc          >  1:1256517/1‑150 (MQ=255)
                                                                                                                                             cTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTc        >  1:422560/1‑150 (MQ=255)
                                                                                                                                              tCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTcc       >  1:858688/1‑150 (MQ=255)
                                                                                                                                                 gAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGc    >  2:689293/1‑150 (MQ=255)
                                                                                                                                                  aaaGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCg   <  1:1034807/150‑1 (MQ=255)
                                                                                                                                                   aaGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGt  >  1:701640/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGT  >  NC_000913/4297048‑4297344

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