Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,953,433 A→G 15.5% intergenic (+20/‑45) ilvE → / → ilvD branched‑chain amino acid aminotransferase/dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,953,4330AG15.5% 251.5 / 33.6 109intergenic (+20/‑45)ilvE/ilvDbranched‑chain amino acid aminotransferase/dihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (36/56);  new base G (9/8);  total (45/64)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGC  >  NC_000913/3953334‑3953530
                                                                                                   |                                                                                                 
gCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGCCGgca                                                                                                   <  1:2248538/100‑1 (MQ=255)
gCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgca                                                                                                   >  2:930011/1‑100 (MQ=255)
gCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgca                                                                                                   >  1:390776/1‑100 (MQ=255)
 cATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgcac                                                                                                  >  1:1228204/1‑100 (MQ=255)
 cATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgcac                                                                                                  >  1:365589/1‑100 (MQ=255)
     cAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgcatgct                                                                                               >  2:1424862/1‑98 (MQ=255)
     cAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgcacgcac                                                                                              >  2:2102649/1‑100 (MQ=255)
          aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                         >  1:1197365/1‑100 (MQ=255)
          aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                         <  1:353305/100‑1 (MQ=255)
           gCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTccc                                                                                        >  1:965064/1‑100 (MQ=255)
            ccTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGTGCTGGTTAGATCCAGCTACTCAATAAATACAAAAAATGGGCCGGCACGCACCGTCCCa                                                                                       >  1:265775/1‑100 (MQ=255)
            ccTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCa                                                                                       <  2:2534808/100‑1 (MQ=255)
            ccTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCa                                                                                       >  1:1577502/1‑100 (MQ=255)
             cttcttCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCAt                                                                                      <  1:1938050/100‑1 (MQ=255)
              ttcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAtt                                                                                     <  1:1018016/100‑1 (MQ=38)
              tccttcGGCCTCTTCACTGGCGAAACCGAAGATAAAGGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGCCGGTGCGTGCCGTCCCAtt                                                                                     <  1:637644/98‑1 (MQ=255)
               tcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCAttt                                                                                    <  2:1211280/100‑1 (MQ=255)
               tcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCAttt                                                                                    >  1:1977963/1‑100 (MQ=255)
                  tcGGCTTCTTCACTGGCGAAACAGAAGAGAAATTGGGTTGGTTAGATCAATTTAATCAATAAATACAAAAAATGTGACGGCACGCACCGTCCCATTTACg                                                                                 <  1:496622/100‑1 (MQ=255)
                  tcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAtttttt                                                                                 <  1:1501743/100‑4 (MQ=37)
                  tcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACg                                                                                 >  2:1040272/1‑100 (MQ=255)
                   cGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttg                                                                                <  2:1126176/100‑5 (MQ=37)
                   cGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttg                                                                                <  2:2331693/100‑5 (MQ=37)
                    ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttgt                                                                               <  2:2482761/100‑6 (MQ=37)
                    ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                               <  1:88122/100‑1 (MQ=255)
                    ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                               <  1:1987214/100‑1 (MQ=255)
                    ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                               <  1:2221852/100‑1 (MQ=255)
                    ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                               <  2:1981717/100‑1 (MQ=255)
                         ctTCACTGGCGAAACCGAAGGTAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagacaga                                                                          >  1:1922756/1‑100 (MQ=255)
                            cACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGacac                                                                       >  2:1938195/1‑100 (MQ=255)
                               tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                    >  2:2485190/1‑100 (MQ=255)
                                ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                   <  2:116251/100‑1 (MQ=255)
                                ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                   <  1:2414330/100‑1 (MQ=255)
                                ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                   <  2:390778/100‑1 (MQ=255)
                                ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                   >  1:2566333/1‑100 (MQ=255)
                                 gCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGa                                                                  >  1:556159/1‑100 (MQ=255)
                                   gAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGt                                                                >  2:1842995/1‑100 (MQ=255)
                                    aaaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAAAAAAAAAAAAAAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGta                                                               <  1:1482748/100‑1 (MQ=255)
                                        cGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaata                                                           <  1:1538391/100‑1 (MQ=255)
                                        cGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaata                                                           >  1:2226506/1‑100 (MQ=255)
                                         gAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaataa                                                          <  2:92106/100‑1 (MQ=255)
                                          aaGATAAAAGGGCCGGTTCAGATCAAGTTAAAAAAAAAAAAAAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaata                                                           <  1:644592/98‑1 (MQ=255)
                                           aGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAg                                                        <  1:2489803/100‑1 (MQ=255)
                                           aGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCACATTTACGAGACAGACACTGGGAGTAAATAAAg                                                        <  1:731708/100‑1 (MQ=255)
                                            gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGt                                                       <  2:2006917/100‑1 (MQ=255)
                                            gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGt                                                       <  1:1838157/100‑1 (MQ=255)
                                               aaaTGGGGCTGGTTAGATCAAGTTAATCAATAAATATAAAAAATGGGACGGCACGCACAGTCACATTTACGAGACAGACACTGGGAGTAAATAAAGTATg                                                    <  1:315055/100‑1 (MQ=255)
                                                aaTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGAAAATAAAGTATGc                                                   >  1:2497337/1‑100 (MQ=255)
                                                 aTGGGTCTGTTTTGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATg                                                    <  2:7120/98‑1 (MQ=255)
                                                 aTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGcc                                                  <  2:1577511/100‑1 (MQ=255)
                                                  tGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCt                                                 >  2:363594/1‑100 (MQ=255)
                                                  tGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCt                                                 <  1:223550/100‑1 (MQ=255)
                                                     ggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                              <  1:2593690/100‑1 (MQ=255)
                                                      gCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCAAGCACCGTCCCATTTACGAGACAGACACTGGGATTAAATAAAGTATGCCTAAGt                                             <  1:1587047/100‑1 (MQ=255)
                                                      gCTGGTTAGATCAAGTTAATCAATAAAAAAAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                              <  1:2222799/99‑1 (MQ=255)
                                                        tGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                           <  2:1067848/100‑1 (MQ=255)
                                                        tGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                           >  1:1869112/1‑100 (MQ=255)
                                                         ggTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAcc                                          >  1:1331633/1‑100 (MQ=255)
                                                           ttAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGt                                        <  2:2383519/100‑1 (MQ=255)
                                                             aGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTc                                      <  2:1139847/100‑1 (MQ=255)
                                                               aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCg                                    <  2:351156/100‑1 (MQ=255)
                                                               aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCg                                    >  1:564012/1‑100 (MQ=255)
                                                                tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGc                                   >  1:1552963/1‑100 (MQ=255)
                                                                 cAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                  <  1:282945/100‑1 (MQ=255)
                                                                 cAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                  >  2:233503/1‑100 (MQ=255)
                                                                  aaGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                  >  1:2078997/1‑99 (MQ=255)
                                                                    gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacc                               <  2:1441711/100‑1 (MQ=255)
                                                                      tAATCAATATATACAAAAAATGGGACGGCACGCACCGTCACATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccan                             <  1:1153537/100‑2 (MQ=255)
                                                                      tAATCAATAAATACAAAAAATGGGACGGCACGCACGGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccac                             <  1:465225/100‑1 (MQ=255)
                                                                       aaTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacc                            >  2:1217919/1‑100 (MQ=255)
                                                                        aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                           >  2:2305092/1‑100 (MQ=255)
                                                                        aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                           >  2:829590/1‑100 (MQ=255)
                                                                        aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                           >  2:829527/1‑100 (MQ=255)
                                                                         tCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCcaccaccac                          <  1:2500183/100‑1 (MQ=255)
                                                                          cgataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACt                         >  1:1619218/3‑100 (MQ=255)
                                                                          cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCcaccaccaca                         >  2:2500202/1‑99 (MQ=255)
                                                                            ataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCa                       <  2:1517832/100‑1 (MQ=255)
                                                                            ataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAACAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCa                       <  1:1759398/100‑1 (MQ=255)
                                                                             taaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCAt                      >  2:1205616/1‑100 (MQ=255)
                                                                               aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACATTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                    <  1:120014/100‑1 (MQ=255)
                                                                               aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                    <  1:363592/100‑1 (MQ=255)
                                                                                  ctcgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCg                 >  1:1086915/6‑100 (MQ=37)
                                                                                   tcgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGt                >  1:172876/5‑100 (MQ=37)
                                                                                   cAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGt                <  2:1444091/100‑1 (MQ=255)
                                                                                    cgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTa               >  2:1023879/4‑100 (MQ=37)
                                                                                    cgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATCCCTAAGTACCGTTCCGCCACAACCACTCATGGTc                  >  1:1410706/4‑97 (MQ=37)
                                                                                    aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTa               <  1:960467/100‑1 (MQ=255)
                                                                                    aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTa               >  1:2439285/1‑100 (MQ=255)
                                                                                     gtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa              >  1:1251933/3‑100 (MQ=38)
                                                                                     gtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa              >  1:2557118/3‑100 (MQ=38)
                                                                                     aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa              >  2:1418562/1‑100 (MQ=255)
                                                                                      taaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAat             >  2:2517024/2‑100 (MQ=38)
                                                                                          tGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAatat           >  1:36629/1‑98 (MQ=39)
                                                                                          tGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATgg         >  1:922475/1‑100 (MQ=255)
                                                                                           gggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGc        >  1:1967974/1‑100 (MQ=255)
                                                                                            ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCg       <  2:1965256/100‑1 (MQ=255)
                                                                                               cGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGt    <  2:2443545/100‑1 (MQ=255)
                                                                                                 gcacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGc  <  2:1824508/100‑1 (MQ=255)
                                                                                                 gcacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGc  <  1:584888/100‑1 (MQ=255)
                                                                                                 gcacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGc  <  1:2102635/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGC  >  NC_000913/3953334‑3953530

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: