Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,915,375 G→C 5.0% intergenic (‑175/+178) glmU ← / ← atpC fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/F1 sector of membrane‑bound ATP synthase, epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,915,3750GC 5.0% 349.2 / 6.8 99intergenic (‑175/+178)glmU/atpCfused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/F1 sector of membrane‑bound ATP synthase, epsilon subunit
Reads supporting (aligned to +/- strand):  ref base G (37/57);  new base C (3/2);  total (40/59)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTT  >  NC_000913/3915227‑3915520
                                                                                                                                                    |                                                                                                                                                 
gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAg                                                                                                                                                   >  4:622840/1‑149 (MQ=255)
   tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTGATCGTAAATTGTACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGAACATCGTCACGCCTTTTTTCTGCAGAAAATTAAGGCAGTAc                                                                                                                                                >  2:755078/1‑149 (MQ=255)
   tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTAc                                                                                                                                                >  2:436697/1‑149 (MQ=255)
   tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTAc                                                                                                                                                >  1:684338/1‑149 (MQ=255)
      aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCtt                                                                                                                                             >  2:606644/1‑149 (MQ=255)
        cccGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGTACTGCCTTaa                                                                                                                                           <  2:452646/148‑1 (MQ=255)
         ccGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCGGTCCCTTAAt                                                                                                                                          <  2:416572/149‑1 (MQ=255)
         ccGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAAt                                                                                                                                          >  4:654784/1‑149 (MQ=255)
          cGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAAtt                                                                                                                                         >  4:128309/1‑149 (MQ=255)
            cTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTa                                                                                                                                       >  2:275310/1‑149 (MQ=255)
             ttCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTaa                                                                                                                                      <  2:267424/149‑1 (MQ=255)
               cACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGc                                                                                                                                    <  3:278832/149‑1 (MQ=255)
                   atGAAAAAAATGTAGTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCa                                                                                                                                >  2:445167/2‑149 (MQ=255)
                      aaaaaaCTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGt                                                                                                                             <  3:153713/149‑1 (MQ=255)
                         aaaCTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGa                                                                                                                          <  3:233752/149‑1 (MQ=255)
                           aCTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATg                                                                                                                        <  4:360827/149‑1 (MQ=255)
                              aCATTTTTTTCATCGTGAAATGGACCGAGGATAAATTGTTCAATAACGCATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGttt                                                                                                                     >  1:499917/1‑149 (MQ=255)
                                aTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGttttt                                                                                                                   <  2:603024/149‑1 (MQ=255)
                                  ttttttCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTAc                                                                                                                                                >  3:526206/1‑118 (MQ=255)
                                  ttttttCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTAc                                                                                                                                                <  4:526206/118‑1 (MQ=255)
                                   tttttCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGct                                                                                                                <  4:33525/149‑1 (MQ=255)
                                   tttttCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGct                                                                                                                <  1:593116/149‑1 (MQ=255)
                                       tCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTa                                                                                                            <  2:350504/149‑1 (MQ=255)
                                          tCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAttt                                                                                                         <  4:654615/149‑1 (MQ=255)
                                          tCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAttt                                                                                                         >  1:298162/1‑149 (MQ=255)
                                          tCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAttt                                                                                                         >  1:66342/1‑149 (MQ=255)
                                           cGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                        <  3:410798/149‑1 (MQ=255)
                                            gTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAtt                                                                                                          <  2:489282/146‑1 (MQ=255)
                                            gTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAtt                                                                                                          >  1:489282/1‑146 (MQ=255)
                                            gTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCg                                                                                                       >  4:445276/1‑149 (MQ=255)
                                              gAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGAt                                                                                                     <  1:233626/149‑1 (MQ=255)
                                                 aTGGACAGAGGATAAATTGTTCAATAACGGATTATCCTCGCTGACGAAAGCGCCATTTTCGACCTTCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAATGCGAACaa                                                                                                  <  4:416459/149‑1 (MQ=255)
                                                 aTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGGAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCaa                                                                                                  <  3:313038/149‑1 (MQ=255)
                                                 aTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                        >  3:529462/1‑143 (MQ=255)
                                                 aTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                        <  4:529462/143‑1 (MQ=255)
                                                  tGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAAt                                                                                                 <  3:409411/149‑1 (MQ=255)
                                                  tGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAAt                                                                                                 <  4:678073/149‑1 (MQ=255)
                                                      cAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAAt                                                                                             >  1:528128/1‑149 (MQ=255)
                                                        gagGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATaa                                                                                           >  3:656297/1‑149 (MQ=255)
                                                           gATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACa                                                                                        <  1:305783/149‑1 (MQ=255)
                                                           gATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACa                                                                                        >  1:767936/1‑149 (MQ=255)
                                                           gATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTACTTCGATCAATGAATAAACa                                                                                        <  3:445092/149‑1 (MQ=255)
                                                            aTAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAg                                                                                       >  2:395560/1‑149 (MQ=255)
                                                            aTAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAg                                                                                       >  3:693968/1‑149 (MQ=255)
                                                                 ttGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGaaaa                                                                                  >  2:592034/1‑149 (MQ=255)
                                                                   gTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACa                                                                                <  4:156584/149‑1 (MQ=255)
                                                                    ttCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTaa                                                                                                                                      >  1:267452/1‑94 (MQ=255)
                                                                    ttCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTaa                                                                                                                                      <  2:267452/94‑1 (MQ=255)
                                                                    ttCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAg                                                                               >  4:276524/1‑149 (MQ=255)
                                                                     tCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgg                                                                              <  1:362743/149‑1 (MQ=255)
                                                                     tCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgg                                                                              <  3:462358/149‑1 (MQ=255)
                                                                      cAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggg                                                                             <  3:204858/149‑1 (MQ=255)
                                                                      cAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggg                                                                             <  1:606644/149‑1 (MQ=255)
                                                                        aTAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggggg                                                                           <  4:505265/149‑1 (MQ=255)
                                                                            cGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGttt                                                                       >  3:588319/1‑149 (MQ=255)
                                                                              gATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTta                                                                     >  3:205091/1‑149 (MQ=255)
                                                                               aTTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTtat                                                                    <  1:48083/149‑1 (MQ=255)
                                                                                   tCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                                <  1:445167/149‑1 (MQ=255)
                                                                                   tCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                                <  3:515867/149‑1 (MQ=255)
                                                                                     cccGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCACAAAATTAAGGCTGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTCTGTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATcagc                                                              <  4:466270/149‑1 (MQ=255)
                                                                                       cGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATcagcag                                                            <  1:275310/149‑1 (MQ=255)
                                                                                         cTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGa                                                          <  2:690635/149‑1 (MQ=255)
                                                                                          tGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGAt                                                         <  2:300347/149‑1 (MQ=255)
                                                                                          tGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGAt                                                         <  3:53265/149‑1 (MQ=255)
                                                                                              gAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                                     >  1:362440/1‑149 (MQ=255)
                                                                                              gAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                                     <  3:208168/149‑1 (MQ=255)
                                                                                              gAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                                     <  2:403634/149‑1 (MQ=255)
                                                                                                aaGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAACCAGGGGGGTTTTATATCAGCAGGATCTAtgt                                                   <  3:662784/149‑1 (MQ=255)
                                                                                                aaGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAtgt                                                   <  4:588319/149‑1 (MQ=255)
                                                                                                    gcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCAATGTGAAc                                               <  2:415819/149‑1 (MQ=255)
                                                                                                       aTTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCt                                            >  3:533130/1‑149 (MQ=255)
                                                                                                          ttCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAtt                                         <  3:175267/149‑1 (MQ=255)
                                                                                                             gACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAg                                      <  2:546910/149‑1 (MQ=255)
                                                                                                              aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGt                                                                                                                                  <  4:582903/56‑1 (MQ=255)
                                                                                                              aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGt                                                                                                                                  >  3:582903/1‑56 (MQ=255)
                                                                                                              aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                     >  2:485841/1‑149 (MQ=255)
                                                                                                                  tCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                 <  1:395560/149‑1 (MQ=255)
                                                                                                                  tCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                 <  3:577529/149‑1 (MQ=255)
                                                                                                                    gTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGt                               <  4:596792/149‑1 (MQ=255)
                                                                                                                         gCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                    <  1:38288/139‑1 (MQ=255)
                                                                                                                         gCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                    >  2:38288/1‑139 (MQ=255)
                                                                                                                                   tGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATcagcag                                                            >  4:451572/1‑105 (MQ=255)
                                                                                                                                   tGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATcagcag                                                            <  3:451572/105‑1 (MQ=255)
                                                                                                                                     cAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa              <  4:272706/149‑1 (MQ=255)
                                                                                                                                     cAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa              <  2:298162/149‑1 (MQ=255)
                                                                                                                                     cAGAAAATTAAAGCAGTAACTTAATTTAAGCGTCAGGTGTATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTAAATCAGCAGGCTCTATGTGAACGCTAATCAGGACGGGTCACACGCGCaaaaaaaa              <  2:312409/149‑1 (MQ=255)
                                                                                                                                      cttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag             >  3:649055/4‑149 (MQ=255)
                                                                                                                                      aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag             >  3:135743/1‑149 (MQ=255)
                                                                                                                                      aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACAGGTCACACGCGCAAAAAAAag             >  2:274477/1‑149 (MQ=255)
                                                                                                                                       ttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                                     <  4:356593/106‑1 (MQ=255)
                                                                                                                                       ttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                                     >  3:356593/3‑108 (MQ=255)
                                                                                                                                       gAAAATTAAGGCAGTACCTTAATTTAAGCGTGAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc            >  4:497296/1‑149 (MQ=255)
                                                                                                                                       gAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc            >  1:451942/1‑149 (MQ=255)
                                                                                                                                       gAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc            <  1:592034/149‑1 (MQ=255)
                                                                                                                                         aaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca          >  3:648956/1‑149 (MQ=255)
                                                                                                                                           aTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc        <  4:693968/149‑1 (MQ=255)
                                                                                                                                            ttAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc       >  3:145100/1‑149 (MQ=255)
                                                                                                                                             tAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCTATAAATGAATAAACAGGAAAACAGGGTGGTTTTTTTTCTGCAGGATCTATGTGAACGATATTCAGGACGGGTTACACGCGCAAAAAAAAGCCATCCt      <  2:207331/149‑1 (MQ=255)
                                                                                                                                                ggCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGtt   <  2:451942/149‑1 (MQ=255)
                                                                                                                                                 gCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGttt  >  1:527063/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTT  >  NC_000913/3915227‑3915520

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