Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,822,103 A→C intergenic (‑92/+36) alaS ← / ← recX alanyl‑tRNA synthetase/regulatory protein for RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,822,1030AC98.4% 236.3 / ‑3.9 64intergenic (‑92/+36)alaS/recXalanyl‑tRNA synthetase/regulatory protein for RecA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (32/31);  minor base T (1/0);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

TATGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGTTCG  >  NC_000913/2821960‑2822250
                                                                                                                                               |                                                                                                                                                   
tATGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCt                                                                                                                                              >  1:750807/1‑151 (MQ=255)
 aTGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCtt                                                                                                                                             >  2:917638/1‑151 (MQ=255)
 aTGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCt                                                                                                                                              <  2:711260/150‑1 (MQ=255)
   ggAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTg                                                                                                                                            <  2:709563/150‑1 (MQ=255)
          gTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTggg                                                                                                                                      >  2:830182/1‑149 (MQ=255)
           tCGAGACACGCCTGCCGGATCTCAGCGGGGCTCTTGCTCCTAATTATCCTGAAATCAAGCTAACGAAATTTCGCCACCAGCGCCAGCGTGTCTTAAACGCCGGGCTGGTAACTGCAAAGTGGGAATAAGCTACGTTTTCTTTACTggg                                                                                                                                      >  2:734343/1‑148 (MQ=255)
           tCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATATGTTTTCTTGACTGGGa                                                                                                                                     >  1:793135/1‑149 (MQ=255)
             gagaAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGt                                                                                                                                  >  1:806356/1‑150 (MQ=255)
             gagaAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTa                                                                                                                                 <  1:153792/151‑1 (MQ=255)
                  aCGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAAt                                                                                                                             >  2:917442/1‑150 (MQ=255)
                   cGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATa                                                                                                                            <  2:911565/150‑1 (MQ=255)
                    gCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATAcc                                                                                                                          >  1:39051/1‑151 (MQ=255)
                        gACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTAt                                                                                                                      <  2:116481/151‑1 (MQ=255)
                           ggATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGc                                                                                                                    >  2:449088/1‑150 (MQ=255)
                           ggATCACAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCg                                                                                                                   <  1:501545/151‑1 (MQ=255)
                             aTCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                                                                                 >  1:149334/1‑151 (MQ=255)
                             aTCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                                                                                 >  1:798955/1‑151 (MQ=255)
                              tctcAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTc                                                                                                                >  1:839577/1‑151 (MQ=255)
                               ctcAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTc                                                                                                                >  1:47080/1‑150 (MQ=255)
                                tcAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTc                                                                                                                >  1:250860/1‑149 (MQ=255)
                                 cAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCa                                                                                                               <  1:146339/149‑1 (MQ=255)
                                 cAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCa                                                                                                               <  1:62618/149‑1 (MQ=255)
                                 cAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGAGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGt                                                                                                             >  2:738838/1‑151 (MQ=255)
                                     ggTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                                                                                         <  1:13155/151‑1 (MQ=255)
                                     ggTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                                                                                         <  1:465956/151‑1 (MQ=255)
                                             tGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAttt                                                                                                  <  2:343765/150‑1 (MQ=255)
                                              gCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTc                                                                                                 <  1:361714/150‑1 (MQ=255)
                                                tCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGcc                                                                                              >  1:344469/1‑151 (MQ=255)
                                                    aaTTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAt                                                                                          >  2:335100/1‑151 (MQ=255)
                                                     aTTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtc                                                                                         <  1:141123/151‑1 (MQ=255)
                                                      ttATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtct                                                                                        >  2:697819/1‑151 (MQ=255)
                                                       tATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtctc                                                                                       >  1:669092/1‑151 (MQ=255)
                                                             gAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTgg                                                                                   <  2:301067/149‑1 (MQ=255)
                                                             gAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                                                                 >  1:522524/1‑151 (MQ=255)
                                                              aaaTCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                                                                                <  1:672086/151‑1 (MQ=255)
                                                                  cAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                                                            <  1:375230/151‑1 (MQ=255)
                                                                   aaGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTc                                                                           <  2:846549/151‑1 (MQ=255)
                                                                     gCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTcc                                                                          >  1:871771/1‑150 (MQ=255)
                                                                            aaaTATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                                                                   <  2:141795/150‑1 (MQ=255)
                                                                               tatCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc                                                               >  1:300502/1‑151 (MQ=255)
                                                                                      accaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGAt                                                          >  1:364163/1‑149 (MQ=255)
                                                                                       ccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag                                                        <  2:371368/150‑1 (MQ=255)
                                                                                       ccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag                                                        <  1:179461/150‑1 (MQ=255)
                                                                                        caGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag                                                        >  2:888286/1‑149 (MQ=255)
                                                                                             ccAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAg                                                   <  2:491468/149‑1 (MQ=255)
                                                                                             ccAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGa                                                  >  1:50995/1‑150 (MQ=255)
                                                                                             ccAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGa                                                  <  1:667727/150‑1 (MQ=255)
                                                                                                gCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaaa                                               <  2:41187/150‑1 (MQ=255)
                                                                                                gCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAAc                                              >  2:871611/1‑151 (MQ=255)
                                                                                                    gtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTg                                          <  1:262286/151‑1 (MQ=255)
                                                                                                       ttAGCCGCAGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTACCATCGGATAGCCACGATAGAGCAGaaa                                                <  2:334988/142‑1 (MQ=255)
                                                                                                       ttAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGAt                                       >  2:922447/1‑151 (MQ=255)
                                                                                                        tAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGAt                                       <  1:908251/150‑1 (MQ=255)
                                                                                                         aaCCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCt                                     <  1:495861/151‑1 (MQ=255)
                                                                                                          aCCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCtt                                    >  1:435924/1‑151 (MQ=255)
                                                                                                             gccgGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAAc                                 >  1:738480/1‑151 (MQ=255)
                                                                                                               cgGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACtt                               <  1:917658/151‑1 (MQ=255)
                                                                                                                 ggCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACtttt                             <  1:335108/151‑1 (MQ=255)
                                                                                                                         aCTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGaaa                      >  1:844431/1‑150 (MQ=255)
                                                                                                                           tGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGaaa                      >  2:90496/1‑148 (MQ=255)
                                                                                                                            gAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAg                    <  1:773863/149‑1 (MQ=255)
                                                                                                                              aaaGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGAc                  >  2:606560/1‑149 (MQ=255)
                                                                                                                                aGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGtt              >  1:41286/1‑151 (MQ=255)
                                                                                                                                   gggAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGc           <  1:680615/151‑1 (MQ=255)
                                                                                                                                        taagatacgTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAgg      >  1:867947/1‑151 (MQ=255)
                                                                                                                                            atacgTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGTTCg  <  2:806338/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TATGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGTTCG  >  NC_000913/2821960‑2822250

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