Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,885,561 (C)5→4 intergenic (‑319/+40) ygcB ← / ← cysH hypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,885,5570C.99.3% 658.2 / ‑3.2 146intergenic (‑315/+44)ygcB/cysHhypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (88/57);  total (89/57)
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.38e-01

Only 100 of 159 total aligned reads displayed.

GCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACG  >  NC_000913/2885466‑2885645
                                                                                           |                                                                                        
 ccTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccc                                                                                        <  1:3941606/93‑1 (MQ=255)
  tttaaGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGTGTTCACGTTTAAAATCAGGAATGCTGCTCTCAAAATAAAAACGcccc                                                                                       >  2:4210007/4‑93 (MQ=255)
  cTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccc                                                                                       >  2:1867252/1‑93 (MQ=255)
  cTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAACCAGGAATGCTGCTCGCAAAATAAAAACGcccc                                                                                       >  2:1867251/1‑93 (MQ=255)
     aaGCATCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcgccatg                                                                                    >  2:155509/1‑90 (MQ=255)
      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                   >  1:4358915/1‑89 (MQ=255)
      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                   <  1:2632197/93‑5 (MQ=255)
        cACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCAtt                                                                                <  1:672185/93‑1 (MQ=255)
         aCCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCAttt                                                                               >  2:2015528/1‑93 (MQ=255)
          cccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGACCCATCAtttt                                                                              >  1:449648/1‑92 (MQ=255)
           ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCt                                                                             <  1:3312465/93‑1 (MQ=255)
           ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCt                                                                             <  2:4306532/93‑1 (MQ=255)
           ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCt                                                                             >  1:319072/1‑93 (MQ=255)
             gACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGa                                                                           >  2:952648/1‑93 (MQ=255)
             gACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGa                                                                           <  2:4044242/93‑1 (MQ=255)
               cGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACa                                                                         >  1:2438259/1‑93 (MQ=255)
               cGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACa                                                                         <  2:3723431/93‑1 (MQ=255)
               cGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACa                                                                         <  2:3723433/93‑1 (MQ=255)
                gCCATTCTATACCTGATAATTCTTCCCGCGGAGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACag                                                                        >  1:1732348/1‑93 (MQ=255)
                 ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                       <  1:4476448/93‑1 (MQ=255)
                 ccATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACaga                                                                       <  2:529800/93‑1 (MQ=255)
                  cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACagag                                                                      >  1:3382632/1‑93 (MQ=255)
                    ttCTATACCTGATAATTCTTCCCGCGGGGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGc                                                                    <  2:2981984/93‑1 (MQ=255)
                    ttCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGc                                                                    >  1:1110835/1‑93 (MQ=255)
                      cTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAAGCGt                                                                  <  1:2942252/93‑1 (MQ=255)
                       tataCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGtt                                                                 <  1:456651/93‑1 (MQ=255)
                          aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGtttca                                                              >  2:2532627/1‑91 (MQ=255)
                          aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTaa                                                              >  2:2532626/1‑93 (MQ=255)
                          aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTaa                                                              <  2:2180563/93‑1 (MQ=255)
                           ccTGATAATTCTTCCCGCGGCGTACACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAt                                                             <  2:2578325/93‑1 (MQ=255)
                             tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAAttt                                                           >  2:3883018/1‑93 (MQ=255)
                              ggtAATTCTTCCCACGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                          >  2:3290853/3‑93 (MQ=255)
                              gATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                          >  1:2483036/1‑93 (MQ=255)
                                 aaTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                       <  1:4323050/93‑1 (MQ=255)
                                   ttcttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCgg                                                     >  1:3203901/1‑93 (MQ=255)
                                   ttcttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCgg                                                     >  1:4410638/1‑93 (MQ=255)
                                    tcttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGACCCATCATTTCTGACAGAGGCGTTTAATTTGTTCGGc                                                    >  2:1832035/1‑93 (MQ=255)
                                      ttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCaa                                                  >  1:4706056/1‑93 (MQ=255)
                                      ttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCaa                                                  >  2:4234537/1‑93 (MQ=255)
                                      ttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCaa                                                  >  2:2377827/1‑93 (MQ=255)
                                        cccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGCCCGTCAata                                                <  2:1185628/93‑1 (MQ=255)
                                         ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                               <  2:676697/93‑1 (MQ=255)
                                         ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                               >  1:2140445/1‑93 (MQ=255)
                                         ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                               >  2:1163796/1‑93 (MQ=255)
                                           gcggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAttt                                             >  2:2015733/1‑93 (MQ=255)
                                             ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                           >  2:525277/1‑93 (MQ=255)
                                             ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                           >  2:1950353/1‑93 (MQ=255)
                                             ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                           <  2:905030/93‑1 (MQ=255)
                                             ggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                           >  1:3640675/1‑93 (MQ=255)
                                              gcgTTCACGTTTAAAATCAGGAATGATGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCTGGCAATATTTAcc                                          >  1:1623457/1‑93 (MQ=255)
                                                gTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                        >  1:106625/1‑93 (MQ=255)
                                                gTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                        <  1:1485927/93‑1 (MQ=255)
                                                 ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                       >  2:1506114/1‑93 (MQ=255)
                                                 ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                       >  2:1587983/1‑93 (MQ=255)
                                                 ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                       >  1:848273/1‑93 (MQ=255)
                                                 ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                       >  2:862390/1‑93 (MQ=255)
                                                  tCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                      >  1:4664730/1‑93 (MQ=255)
                                                  tCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                      >  1:395163/1‑93 (MQ=255)
                                                     cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtg                                   <  1:1006140/93‑1 (MQ=255)
                                                      gTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtgt                                  <  2:4717953/93‑1 (MQ=255)
                                                      gTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtgt                                  >  1:4227324/1‑93 (MQ=255)
                                                       tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATTATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                                 <  2:1072436/93‑1 (MQ=255)
                                                       tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                                 <  2:1072437/93‑1 (MQ=255)
                                                       tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                                 >  1:3547665/1‑93 (MQ=255)
                                                        ttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTaa                                <  2:1251936/93‑1 (MQ=255)
                                                          aaaaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAAcc                              <  2:1376327/93‑1 (MQ=255)
                                                          aaaaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAAcc                              >  2:923529/1‑93 (MQ=255)
                                                          aaaaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTAAATTTGTCCGGCAATATTTACCCTTCGTGTAAcc                              <  2:2976629/93‑1 (MQ=255)
                                                               cAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAt                         >  2:1715257/1‑93 (MQ=255)
                                                               cAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAt                         <  1:336509/93‑1 (MQ=255)
                                                               cAGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAt                         <  2:4679099/93‑1 (MQ=255)
                                                                aGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                        >  1:2595838/1‑93 (MQ=255)
                                                                aGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                        >  2:2712448/1‑93 (MQ=255)
                                                                aGGAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                        >  1:2239809/1‑93 (MQ=255)
                                                                 ggAATGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTc                       >  2:2064113/1‑93 (MQ=255)
                                                                    aTGCTGCTCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCt                    >  1:1507611/1‑93 (MQ=255)
                                                                     tgctgcTCGCAAAATAAAAACGCCCCATCATTTATGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtt                   >  2:1965014/1‑93 (MQ=255)
                                                                      gctgctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCttt                  <  1:2860682/93‑1 (MQ=255)
                                                                      gctgctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCttt                  <  1:2752037/93‑1 (MQ=255)
                                                                       ctgctCGCAAAATAAAAACGCCCCATCATGTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtttt                 <  2:2573019/93‑1 (MQ=255)
                                                                          ctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTaa              >  2:772583/1‑93 (MQ=255)
                                                                          ctCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTaa              >  1:578853/1‑93 (MQ=255)
                                                                           tCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg             >  1:451022/1‑93 (MQ=255)
                                                                           tCGCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg             >  2:3847510/1‑93 (MQ=255)
                                                                             gCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc           <  2:4322883/93‑1 (MQ=255)
                                                                             gCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc           <  1:4390082/93‑1 (MQ=255)
                                                                             gCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc           >  2:2570012/1‑93 (MQ=255)
                                                                             gCAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTACGGGAATATTTACCCTTCGTGTAACCCACATACCCTTTTTAAGcc           <  1:3268305/93‑1 (MQ=255)
                                                                              cAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCa          <  2:3556986/93‑1 (MQ=255)
                                                                              cAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCa          <  2:3556984/93‑1 (MQ=255)
                                                                              cAAAATAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCa          >  2:2141223/1‑93 (MQ=255)
                                                                                 aaTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAg       <  1:645079/93‑1 (MQ=255)
                                                                                  aTAAAAACGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGa      <  2:131420/93‑1 (MQ=255)
                                                                                    aaaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa    >  1:2601071/1‑93 (MQ=255)
                                                                                    aaaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa    <  2:2647332/93‑1 (MQ=255)
                                                                                    aaaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa    <  1:2087207/93‑1 (MQ=255)
                                                                                     aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc   <  1:3509498/93‑1 (MQ=255)
                                                                                     aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc   <  2:439921/93‑1 (MQ=255)
                                                                                     aaaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc   <  1:2277724/93‑1 (MQ=255)
                                                                                      aaaCGCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAcg  <  1:1134339/93‑1 (MQ=255)
                                                                                           |                                                                                        
GCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACG  >  NC_000913/2885466‑2885645

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