Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,359 T→A 18.0% intergenic (+53/‑310) yjbM → / → dusA uncharacterized protein/tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,3590TA18.0% 313.1 / 57.0 132intergenic (+53/‑310)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (53/55);  new base A (14/10);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 136 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458
                                                                                                   |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                                     >  2:2248963/1‑100 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                                <  2:909916/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                             <  2:1490562/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                             <  2:198086/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                             <  2:2304117/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                          <  1:1497497/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                          <  2:584364/100‑1 (MQ=38)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttt                                                                                        >  1:696853/1‑100 (MQ=255)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCttt                                                                                        <  1:61594/100‑1 (MQ=39)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCtttt                                                                                       <  2:199818/100‑1 (MQ=39)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttt                                                                                       <  1:1901813/100‑1 (MQ=38)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                          <  1:2580190/97‑1 (MQ=39)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                      <  1:1390100/100‑1 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGATAAATAATATTTAAAGTGCTTCGATAGTATGGAAGGCATGATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                      <  1:527377/99‑1 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                     >  2:1779551/1‑100 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                     >  1:2142060/1‑100 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  1:1329297/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCttttata                                                                                    <  2:913913/100‑4 (MQ=38)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATATGAAGGTCTTCGGTAGTATGGAATGCATTATATAAACGGACGCAGCATTTATTGGGTTCtttct                                                                                      >  1:1862378/1‑96 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                   <  2:2619744/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                   <  2:1300511/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                   >  1:1834302/1‑100 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatat                                                                                   <  1:1492341/100‑5 (MQ=37)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                                 <  2:430930/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                                >  2:1199992/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataatg                                                                               <  2:2851225/100‑9 (MQ=25)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGTAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGCCCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                              <  1:1614382/100‑1 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                              <  1:1689932/100‑1 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                              <  2:1317839/100‑1 (MQ=255)
                        tCGGGGGTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATCGGGTCCttttatataatgcc                                                                             <  1:1776379/100‑11 (MQ=21)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  2:702935/1‑100 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                        >  1:1188282/1‑100 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTGTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     >  1:1668309/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  2:333469/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    <  2:15106/100‑1 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                   >  2:57128/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  <  2:193003/100‑1 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  <  2:1378421/100‑1 (MQ=255)
                                   tCGGAGAAATAATATGTAAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAgn                                                                  <  1:2080238/100‑2 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  <  2:130701/99‑1 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                >  2:1556169/1‑100 (MQ=255)
                                       agaAATAATATGTGACGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                              <  2:2872754/100‑1 (MQ=255)
                                        gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                             >  2:2879885/1‑100 (MQ=255)
                                        gaAATAATATGTGAAGTTCTTAGATAGTATGGAGGCCATTATATAAATGGATACAATATTAACTGGGTTCTTTTTTATCTATATATACTact                                                                     >  1:1004448/1‑90 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                            <  1:2899202/100‑1 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                            <  2:1938037/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  1:2885011/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  2:1165099/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         >  1:2666202/1‑100 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       <  2:1834310/100‑1 (MQ=255)
                                                      aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATAGTACCaga                                               <  1:2026612/100‑1 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATATAAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                           <  1:2789181/100‑1 (MQ=255)
                                                           ttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTa                                          >  2:599192/1‑100 (MQ=255)
                                                            tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATTACCagag                                              >  1:1741732/1‑95 (MQ=255)
                                                               aTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgt                                      >  1:1362023/1‑100 (MQ=255)
                                                                tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  2:1115473/100‑1 (MQ=255)
                                                                tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  2:850402/100‑1 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                   >  2:1174851/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                   >  2:675575/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                 >  2:938078/1‑100 (MQ=255)
                                                                      ggAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                               <  2:2437363/100‑1 (MQ=255)
                                                                       gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                              <  1:507941/100‑1 (MQ=255)
                                                                         aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                            <  1:2783278/100‑1 (MQ=255)
                                                                          ggCATTATATAAAAGGACACAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                            <  1:2762141/99‑1 (MQ=255)
                                                                           gatagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          >  2:113517/10‑100 (MQ=21)
                                                                           gCATTATATACACGGCCCCAATATTTATTGGGTTCTTTTTTCTATATCAATAATATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                           <  2:133445/99‑1 (MQ=255)
                                                                           gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  1:1058935/100‑1 (MQ=255)
                                                                            atagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         >  2:347050/9‑100 (MQ=25)
                                                                            aatgatataAAATGAGCCAATATTTATTGGGGTCTTTTTTCTCTGTCGATACTTTTAGCAGGGAGATATAGTACCAGAGTTTAATGTGTGATTTtttatt                         >  1:243273/5‑100 (MQ=255)
                                                                             tagagaaaaaaGAACCCAATAAATATTGGGTCATTTTTTCTGTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        >  1:1586796/8‑100 (MQ=25)
                                                                             aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGCGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        <  1:1554163/100‑1 (MQ=255)
                                                                              agagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                       >  2:1894478/7‑100 (MQ=25)
                                                                               tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGATATATATCACCAGAGTTTAATGTGTGATTTtttattta                       <  1:773439/99‑1 (MQ=255)
                                                                               tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      <  1:2337108/100‑1 (MQ=255)
                                                                               tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      <  1:2336992/100‑1 (MQ=255)
                                                                               gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:88239/6‑100 (MQ=25)
                                                                                atataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                     <  2:1362028/100‑1 (MQ=255)
                                                                                agaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATAACCAGAGTTTCATGTGTGATTTTTTATTTAtc                     >  2:2470354/5‑100 (MQ=25)
                                                                                 tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCCGGGAGCTATATCCCCAGAATTTAATGTGTGATTTTTTATTTAtc                     >  1:1335654/1‑99 (MQ=255)
                                                                                 gaaaaaataaCCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACCATTAGCAGGGAGATATATAACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                    >  1:1011799/10‑100 (MQ=37)
                                                                                 gaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                    >  1:2197555/4‑100 (MQ=37)
                                                                                   taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  1:1184648/1‑100 (MQ=255)
                                                                                   aaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  2:877793/2‑100 (MQ=39)
                                                                                   aaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  2:2281055/2‑100 (MQ=38)
                                                                                    aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                 >  2:370463/1‑100 (MQ=39)
                                                                                     aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa                >  1:2618465/1‑100 (MQ=255)
                                                                                      aaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa               >  2:1489940/1‑100 (MQ=255)
                                                                                        ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAcc             <  1:688974/100‑1 (MQ=255)
                                                                                          aaccAATATTTATGGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           >  2:1974245/3‑100 (MQ=255)
                                                                                          aCCCAATATTTATTGGGTTCTTGTTTCTCTATCAATACTATTAGCAGTGAGATATATCACCAGAGTTTAATGTGTGATTTTGTATTTATCGTCGAACCTg           <  1:939716/100‑1 (MQ=255)
                                                                                           cccAATATTTATTGGGTTCTTTTTTTTTTTTTTTTCCTATTAGCAGGGCGATATATCACCTTAGTTTAATGTGTGATTTTTTATTTAtcctccccc              >  1:67549/1‑92 (MQ=255)
                                                                                           cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  2:572571/1‑100 (MQ=255)
                                                                                            ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGACCCTGGa         >  2:2644416/1‑100 (MQ=255)
                                                                                            ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa         <  1:2160940/100‑1 (MQ=255)
                                                                                             cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt        >  2:1681251/1‑100 (MQ=255)
                                                                                             cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTCTGATTTTTTATTTTTCTTCTAACCt            >  1:1750606/1‑96 (MQ=255)
                                                                                              aaTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt       <  2:1414258/100‑1 (MQ=255)
                                                                                                  tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   <  2:2738789/100‑1 (MQ=255)
                                                                                                   ttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTa  >  1:1205709/1‑100 (MQ=255)
                                                                                                   ttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTa  <  2:2891992/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458

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