Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,922,656 T→C intergenic (‑216/+401) atpI ← / ← rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,922,6560TC98.7% 293.6 / NA 77intergenic (‑216/+401)atpI/rsmGATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (49/27);  minor base A (1/0);  total (50/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATG  >  NC_000913/3922511‑3922802
                                                                                                                                                 |                                                                                                                                                  
cAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATa                                                                                                                                               <  2:368412/151‑1 (MQ=255)
   gCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAAcg                                                                                                                                            >  2:709013/1‑151 (MQ=255)
     aTCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAAcgcg                                                                                                                                          <  2:382487/151‑1 (MQ=255)
      tCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGt                                                                                                                                         >  1:436075/1‑151 (MQ=255)
       cAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGt                                                                                                                                         <  2:574687/150‑1 (MQ=255)
        aaaaaGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTa                                                                                                                                        >  2:603219/1‑150 (MQ=255)
               ggTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTgg                                                                                                                                >  2:234373/1‑151 (MQ=255)
                  cAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTg                                                                                                                              >  1:695034/1‑150 (MQ=255)
                  cAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGa                                                                                                                             <  1:481858/151‑1 (MQ=255)
                    aaTTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGAt                                                                                                                            >  2:339957/1‑150 (MQ=255)
                    aaTTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGAt                                                                                                                            >  1:244163/1‑150 (MQ=255)
                     aTTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATa                                                                                                                           >  2:9732/1‑150 (MQ=255)
                       tataCGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATaaa                                                                                                                         <  1:301556/150‑1 (MQ=255)
                           cGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTg                                                                                                                    >  1:390766/1‑151 (MQ=255)
                             gTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCt                                                                                                                  <  2:399616/151‑1 (MQ=255)
                               gcgcCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTc                                                                                                                >  1:498716/1‑151 (MQ=255)
                                cgcCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTc                                                                                                                >  1:190051/1‑150 (MQ=255)
                                 gcCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTc                                                                                                                >  2:219749/1‑149 (MQ=255)
                                  cccccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGa                                                                                                             <  2:229467/151‑1 (MQ=255)
                                      cGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATcac                                                                                                         >  1:242919/1‑151 (MQ=255)
                                         gATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACat                                                                                                       >  2:192607/1‑150 (MQ=255)
                                         gATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACat                                                                                                       >  2:698319/1‑150 (MQ=255)
                                          aTTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCGGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACatat                                                                                                     >  2:492886/1‑151 (MQ=255)
                                          aTTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACatat                                                                                                     >  1:421199/1‑151 (MQ=255)
                                          aTTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACatat                                                                                                     >  1:93848/1‑151 (MQ=255)
                                              cAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTgcg                                                                                                 <  1:153187/151‑1 (MQ=255)
                                              cAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTgcg                                                                                                 >  2:47316/1‑151 (MQ=255)
                                               aaaCAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTgcgc                                                                                                >  2:629636/1‑151 (MQ=255)
                                                   aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGt                                                                                            >  2:785249/1‑151 (MQ=255)
                                                   aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGt                                                                                            >  1:109593/1‑151 (MQ=255)
                                                          aGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCgcgc                                                                                      >  2:407194/1‑150 (MQ=255)
                                                            ccAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGtt                                                                                           >  1:75675/1‑143 (MQ=255)
                                                            ccAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGtt                                                                                           <  2:75698/143‑1 (MQ=255)
                                                                aaGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCat                                                                               <  1:586086/151‑1 (MQ=255)
                                                                     gAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCAtattta                                                                          <  2:34347/151‑1 (MQ=255)
                                                                       aTAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCAtatttatt                                                                        >  2:102904/1‑151 (MQ=255)
                                                                       aTAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCAtatttatt                                                                        >  2:745785/1‑151 (MQ=255)
                                                                           aTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGcaca                                                                                   >  1:795126/1‑136 (MQ=255)
                                                                           aTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGcaca                                                                                   <  2:795182/136‑1 (MQ=255)
                                                                                 aaaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCa                                                               <  2:182044/150‑1 (MQ=255)
                                                                                 aaaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCa                                                               >  1:72625/1‑150 (MQ=255)
                                                                                  aaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCaaa                                                             <  2:444677/151‑1 (MQ=255)
                                                                                   atatTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCaaa                                                             >  1:776354/1‑150 (MQ=255)
                                                                                    tatTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATACCCCGTAAACACTGGTGATAAAGCGGGCTCCAGATCAAATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCaatt                                                            >  2:397240/1‑148 (MQ=255)
                                                                                    tatTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatg                                                           >  2:820892/1‑151 (MQ=255)
                                                                                    tatTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatg                                                           >  2:432133/1‑151 (MQ=255)
                                                                                            cAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCttt                                                   <  1:402594/151‑1 (MQ=255)
                                                                                            cAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGAAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCttt                                                   >  1:701244/1‑151 (MQ=255)
                                                                                             aGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTg                                                  >  2:782882/1‑151 (MQ=255)
                                                                                               tGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGc                                                 >  1:412297/1‑150 (MQ=255)
                                                                                                gCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCaa                                               <  1:857925/151‑1 (MQ=255)
                                                                                                 cATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAg                                              <  1:72708/151‑1 (MQ=255)
                                                                                                  aTACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGt                                             >  1:239420/1‑151 (MQ=255)
                                                                                                   tACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGtt                                            >  2:679808/1‑151 (MQ=255)
                                                                                                    aCAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGtt                                            <  1:692261/150‑1 (MQ=255)
                                                                                                      aaTTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatga                                                          <  2:300247/134‑1 (MQ=255)
                                                                                                      aaTTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatga                                                          >  1:300222/1‑134 (MQ=255)
                                                                                                      aaTTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTAt                                         >  2:659650/1‑151 (MQ=255)
                                                                                                      aaTTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTAt                                         <  1:698276/151‑1 (MQ=255)
                                                                                                            gACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATAtt                                   >  1:586083/1‑151 (MQ=255)
                                                                                                             aCTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATAttt                                  >  1:678372/1‑151 (MQ=255)
                                                                                                             aCTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATAttt                                  >  2:613782/1‑151 (MQ=255)
                                                                                                                   cTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCt                            >  1:438927/1‑151 (MQ=255)
                                                                                                                     gCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCt                            <  2:49133/149‑1 (MQ=255)
                                                                                                                     gCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTa                           <  1:140983/150‑1 (MQ=255)
                                                                                                                     gCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTa                           >  2:452558/1‑150 (MQ=255)
                                                                                                                        aaCCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaa                        <  1:50444/150‑1 (MQ=255)
                                                                                                                        aaCCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaa                        >  2:127678/1‑150 (MQ=255)
                                                                                                                        aaCCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaaaa                        >  1:68649/1‑150 (MQ=255)
                                                                                                                              gTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGAt                  >  2:720175/1‑150 (MQ=255)
                                                                                                                                tCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACt               >  1:723773/1‑151 (MQ=255)
                                                                                                                                     cATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGtttt          <  1:27744/151‑1 (MQ=255)
                                                                                                                                     cATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGtttt          >  2:326138/1‑151 (MQ=255)
                                                                                                                                        gAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTtaat       <  2:667153/151‑1 (MQ=255)
                                                                                                                                           aaaCTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATaaa    <  2:671165/151‑1 (MQ=255)
                                                                                                                                            aaCTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTCGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAAt   <  1:117413/151‑1 (MQ=255)
                                                                                                                                              cTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATg  >  2:811218/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATG  >  NC_000913/3922511‑3922802

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