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.5% 252.5 / ‑4.5 67intergenic (‑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 (1/0);  new base C (40/26);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAAT  >  NC_000913/3922509‑3922801
                                                                                                                                                   |                                                                                                                                                 
gTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCaaa                                                                                                                                                <  1:613511/151‑1 (MQ=255)
    aGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAAc                                                                                                                                            <  1:275619/151‑1 (MQ=255)
     gCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGTTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAAcg                                                                                                                                           <  1:663488/151‑1 (MQ=255)
       aTCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAAcgcg                                                                                                                                         <  2:374361/151‑1 (MQ=255)
          aacaaGCGGTCACATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTa                                                                                                                                       >  1:60820/4‑150 (MQ=255)
            aaaGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAac                                                                                                                                    <  1:512095/151‑1 (MQ=255)
             aaGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAaca                                                                                                                                   <  2:324194/151‑1 (MQ=255)
                cGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTg                                                                                                                                >  2:574185/1‑151 (MQ=255)
                cGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTg                                                                                                                                <  1:615233/151‑1 (MQ=255)
                   tCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTg                                                                                                                             >  1:640474/1‑151 (MQ=255)
                      aaTTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGAt                                                                                                                           >  2:338780/1‑150 (MQ=255)
                       aTTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATa                                                                                                                          >  2:173633/1‑150 (MQ=255)
                         tataCGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATaaa                                                                                                                        <  1:549983/150‑1 (MQ=255)
                            aCGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGt                                                                                                                    <  2:158594/151‑1 (MQ=255)
                            aCGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGt                                                                                                                    <  1:45138/150‑1 (MQ=255)
                             cGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTg                                                                                                                   >  2:206577/1‑151 (MQ=255)
                                tGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCtt                                                                                                                >  1:302389/1‑151 (MQ=255)
                                    cccccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGa                                                                                                            <  2:443967/151‑1 (MQ=255)
                                     ccccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTc                                                                                                               >  2:96040/1‑147 (MQ=255)
                                       ccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATc                                                                                                          >  1:751855/1‑150 (MQ=255)
                                            aTTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACatat                                                                                                    >  1:33491/1‑151 (MQ=255)
                                             tttCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATAtt                                                                                                   >  1:389105/1‑151 (MQ=255)
                                                   aCAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCAt                                                                                             <  2:285497/151‑1 (MQ=255)
                                                   aCAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCAt                                                                                             >  2:233947/1‑151 (MQ=255)
                                                    cAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATg                                                                                            >  1:70400/1‑151 (MQ=255)
                                                    cAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATg                                                                                            >  2:745375/1‑151 (MQ=255)
                                                    aaaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATg                                                                                            >  2:79149/2‑151 (MQ=255)
                                                      aTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGtt                                                                                          >  1:361020/1‑151 (MQ=255)
                                                       tAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTATCCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAGCCACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGtt                                                                                          <  1:377933/150‑1 (MQ=255)
                                                         aGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTcgc                                                                                       <  1:558867/151‑1 (MQ=255)
                                                         aGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTcgc                                                                                       >  2:616409/1‑151 (MQ=255)
                                                           tAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCgcgc                                                                                     >  1:670724/1‑151 (MQ=255)
                                                           tAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCgcgc                                                                                     >  2:468550/1‑151 (MQ=255)
                                                                aaaaGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGc                                                                                >  2:575591/1‑151 (MQ=255)
                                                                 aaaGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGc                                                                                >  2:6196/1‑150 (MQ=255)
                                                                 aaaGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCa                                                                               <  2:670782/151‑1 (MQ=255)
                                                                      tGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCAtattt                                                                          >  1:569408/1‑151 (MQ=255)
                                                                           aaaTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTgcgc                                                                                               <  2:469772/125‑1 (MQ=255)
                                                                           aaaTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTgcgc                                                                                               >  1:469721/1‑125 (MQ=255)
                                                                           aaaTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATatttattt                                                                      >  2:299297/1‑150 (MQ=255)
                                                                           aaaTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATAtttattta                                                                     <  2:509657/151‑1 (MQ=255)
                                                                            aaTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTAc                                                                    <  1:216913/151‑1 (MQ=255)
                                                                              tGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACtt                                                                  >  1:183663/1‑151 (MQ=255)
                                                                               gTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTg                                                                 >  2:453664/1‑151 (MQ=255)
                                                                                tttAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTgg                                                                >  1:645896/1‑151 (MQ=255)
                                                                                   aaaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCa                                                              >  2:618985/1‑150 (MQ=255)
                                                                                    aaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCaaa                                                            <  2:139996/151‑1 (MQ=255)
                                                                                     atatTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCTTGCTTCAGATCACATATTGCGCATGTTCTCGCACAGCATATTTATTTACTTGGCAAat                                                           <  1:738892/151‑1 (MQ=255)
                                                                                     atatTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAat                                                           >  1:649511/1‑151 (MQ=255)
                                                                                       atTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatg                                                          >  1:394815/1‑150 (MQ=255)
                                                                                       atTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatg                                                          >  2:744941/1‑150 (MQ=255)
                                                                                         ttttCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatgatg                                                       >  1:434084/1‑151 (MQ=255)
                                                                                         ttttCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAatgat                                                        <  2:394864/150‑1 (MQ=255)
                                                                                            tCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCt                                                    >  2:292093/1‑151 (MQ=255)
                                                                                                   cATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCa                                               >  2:429106/1‑149 (MQ=255)
                                                                                                               aCTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATAttt                                 <  2:492962/151‑1 (MQ=255)
                                                                                                                 ttttCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTc                                >  2:648809/1‑150 (MQ=255)
                                                                                                                  tttCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCa                               <  2:302422/150‑1 (MQ=255)
                                                                                                                   ttCTGCTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGt                             >  2:245835/1‑151 (MQ=255)
                                                                                                                        cTAACCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTaa                         >  2:766744/1‑150 (MQ=255)
                                                                                                                           aCCCTGTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAAc                     >  2:284980/1‑151 (MQ=255)
                                                                                                                                gTTCGATCATGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATa                <  2:40457/151‑1 (MQ=255)
                                                                                                                                     aTCCTGAAAAACTGCAAATAATGCGTAAACACTGGTGATAAAGAGTGCTTCAGATCACATATTACGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCATTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGtt           >  1:721941/1‑151 (MQ=255)
                                                                                                                                        aTGAAAAACTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGtttt         >  2:421848/1‑150 (MQ=255)
                                                                                                                                           aaaaaCTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTtaata     <  1:245802/151‑1 (MQ=255)
                                                                                                                                            aaaaCTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTtaataa    <  2:225494/151‑1 (MQ=255)
                                                                                                                                              aaCTGCAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAAt  >  2:557012/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAAT  >  NC_000913/3922509‑3922801

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