Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,953,437 A→G 15.5% intergenic (+24/‑41) 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,4370AG15.5% 251.3 / 34.3 109intergenic (+24/‑41)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 115 total aligned reads displayed.

CAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTG  >  NC_000913/3953339‑3953535
                                                                                                  |                                                                                                  
cAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGgcacgcac                                                                                                   >  2:2102649/1‑100 (MQ=255)
     aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                              <  1:348509/100‑1 (MQ=255)
     aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                              >  1:1197365/1‑100 (MQ=255)
     aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                              <  1:353305/100‑1 (MQ=255)
     aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                              >  1:168728/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)
          tcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCAttt                                                                                         >  1:1977963/1‑100 (MQ=255)
          tcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCAttt                                                                                         <  2:1211280/100‑1 (MQ=255)
             tcGGCTTCTTCACTGGCGAAACAGAAGAGAAATTGGGTTGGTTAGATCAATTTAATCAATAAATACAAAAAATGTGACGGCACGCACCGTCCCATTTACg                                                                                      <  1:496622/100‑1 (MQ=255)
             tcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAtttttt                                                                                      <  1:1501743/100‑4 (MQ=37)
             tcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACg                                                                                      <  2:2345901/100‑1 (MQ=255)
             tcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACg                                                                                      >  2:1040272/1‑100 (MQ=255)
              cGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttg                                                                                     <  2:2331693/100‑5 (MQ=37)
              cGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttg                                                                                     <  2:1126176/100‑5 (MQ=37)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttgt                                                                                    <  1:1241194/100‑6 (MQ=37)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttgt                                                                                    <  2:113855/100‑6 (MQ=37)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttttttgt                                                                                    <  2:2482761/100‑6 (MQ=37)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                                    <  2:1376974/100‑1 (MQ=255)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                                    <  1:1987214/100‑1 (MQ=255)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                                    <  1:88122/100‑1 (MQ=255)
               ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgag                                                                                    <  1:2221852/100‑1 (MQ=255)
                    ctTCACTGGCGAAACCGAAGGTAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagacaga                                                                               >  1:1922756/1‑100 (MQ=255)
                     ttCACTGGCTAAACCGAAGATAAATGGGGCTGGTTGGATCAAGTTAATCAATAAAGACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagac                                                                                  <  1:956857/96‑1 (MQ=255)
                       cACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGacac                                                                            >  2:1938195/1‑100 (MQ=255)
                          tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                         <  1:412457/100‑1 (MQ=255)
                          tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                         >  2:2485190/1‑100 (MQ=255)
                           ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                        <  2:2505349/100‑1 (MQ=255)
                           ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                        <  2:390778/100‑1 (MQ=255)
                           ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                        <  2:116251/100‑1 (MQ=255)
                           ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTggg                                                                        <  1:1963758/100‑1 (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)
                                      aGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAg                                                             <  1:1532688/100‑1 (MQ=255)
                                      aGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCACATTTACGAGACAGACACTGGGAGTAAATAAAg                                                             <  1:731708/100‑1 (MQ=255)
                                       gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGt                                                            <  1:1838157/100‑1 (MQ=255)
                                       gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGt                                                            <  2:2006917/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                                                      <  1:223550/100‑1 (MQ=255)
                                              ggggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTa                                                     <  2:1977200/100‑1 (MQ=255)
                                                ggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                                   <  1:2593690/100‑1 (MQ=255)
                                                 gCTGGTTAGATCAAGTTAATCAATAAAAAAAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                                   <  1:2222799/99‑1 (MQ=255)
                                                   tGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                <  2:1067848/100‑1 (MQ=255)
                                                      ttAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGt                                             <  2:2383519/100‑1 (MQ=255)
                                                        aGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTc                                           <  2:1139847/100‑1 (MQ=255)
                                                          aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCg                                         >  1:564012/1‑100 (MQ=255)
                                                          aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCg                                         <  2:351156/100‑1 (MQ=255)
                                                           tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGc                                        >  1:1552963/1‑100 (MQ=255)
                                                            cAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                       >  2:233503/1‑100 (MQ=255)
                                                            cAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                       <  1:282945/100‑1 (MQ=255)
                                                             aaGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                       >  1:2078997/1‑99 (MQ=255)
                                                               gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacc                                    <  2:1441711/100‑1 (MQ=255)
                                                                ttAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacca                                   >  1:2502273/1‑100 (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:829527/1‑100 (MQ=255)
                                                                   aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                                >  2:2305092/1‑100 (MQ=255)
                                                                   aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                                >  2:829590/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)
                                                                          aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACATTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                         <  1:120014/100‑1 (MQ=255)
                                                                          aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                         <  1:363592/100‑1 (MQ=255)
                                                                           aaaaaaaaaaTGGGACGGCACGGACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACGCATGGt                        <  1:2139874/96‑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                    <  2:55348/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)
                                                                                     tGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATgg              >  1:2304622/1‑100 (MQ=255)
                                                                                     tGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATgg              >  1:922475/1‑100 (MQ=255)
                                                                                      gggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGc             >  1:1967974/1‑100 (MQ=255)
                                                                                       ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCg            <  2:1965256/100‑1 (MQ=255)
                                                                                           ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTg        >  1:621247/1‑100 (MQ=255)
                                                                                            gcacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGc       <  2:1824508/100‑1 (MQ=255)
                                                                                            gcacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGc       <  1:584888/100‑1 (MQ=255)
                                                                                               cgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCg    >  1:1840021/1‑100 (MQ=255)
                                                                                                 cacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTg  <  1:2332956/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTG  >  NC_000913/3953339‑3953535

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: