Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,953,437 A→G 6.2% 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,4370AG 6.2% 237.4 / 5.2 80intergenic (+24/‑41)ilvE/ilvDbranched‑chain amino acid aminotransferase/dihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (30/45);  new base G (3/2);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

TTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAT  >  NC_000913/3953295‑3953584
                                                                                                                                              |                                                                                                                                                   
ttCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCACACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGCCTGGGTAGATCAAGTTAATCAATAAATACAAAAAATGTGACGGCACGCACCGTcc                                                                                                                                               >  1:1333262/1‑149 (MQ=255)
ttCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGTTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                                                                               >  2:1324320/1‑149 (MQ=255)
ttCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                                                                               <  1:1022029/149‑1 (MQ=255)
 tCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGTCCTCTTCACTGGCGAAAGCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTccc                                                                                                                                              <  2:112565/149‑1 (MQ=255)
 tCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTccc                                                                                                                                              <  2:988568/149‑1 (MQ=255)
 tCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTccc                                                                                                                                              <  1:1237194/149‑1 (MQ=255)
  cAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCa                                                                                                                                             >  2:1346844/1‑149 (MQ=255)
           gAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgaga                                                                                                                                    <  1:708714/149‑1 (MQ=255)
           gAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgaga                                                                                                                                    <  1:1221505/149‑1 (MQ=255)
                     gtgGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                          <  1:440464/149‑1 (MQ=255)
                     gtgGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                          <  2:71699/149‑1 (MQ=255)
                     gtgGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                          <  2:227414/149‑1 (MQ=255)
                     gtgGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                          <  1:567193/149‑1 (MQ=255)
                       ggCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTACTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGa                                                                                                                        <  2:316737/149‑1 (MQ=255)
                          ccGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGta                                                                                                                     <  1:1421091/149‑1 (MQ=255)
                           cGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaa                                                                                                                    <  2:1227496/149‑1 (MQ=255)
                               tACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaataa                                                                                                                >  2:1312695/1‑149 (MQ=255)
                                   aaaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTa                                                                                                            >  2:1506891/1‑149 (MQ=255)
                                    aaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATTGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTAt                                                                                                           >  1:1154539/1‑149 (MQ=255)
                                    aaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTAt                                                                                                           <  1:912544/149‑1 (MQ=255)
                                    aaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTAt                                                                                                           <  1:66773/149‑1 (MQ=255)
                                      cGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGc                                                                                                         <  2:908724/149‑1 (MQ=255)
                                       gCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGcc                                                                                                        <  2:963965/149‑1 (MQ=255)
                                       gCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGcc                                                                                                        <  2:182858/149‑1 (MQ=255)
                                         aTTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTa                                                                                                      >  1:1386382/1‑149 (MQ=255)
                                              gCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                                                                 >  2:1128640/1‑149 (MQ=255)
                                              gCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                                                                 <  1:1098609/149‑1 (MQ=255)
                                               cAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAcc                                                                                                >  2:989006/1‑149 (MQ=255)
                                                 aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCCCGCACCGTCCCATTTACGAGACAGACACTGGGAGTATATAAAGAATGCCTAAGTACCGt                                                                                              >  1:1079458/1‑149 (MQ=255)
                                                    cttcttCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                                                           >  2:1022447/1‑149 (MQ=255)
                                                      tcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGc                                                                                         <  2:918329/149‑1 (MQ=255)
                                                           ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacca                                                                                    <  2:1386382/149‑1 (MQ=255)
                                                           ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacca                                                                                    <  1:1525415/149‑1 (MQ=255)
                                                             ccTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACTGACACTGGGAGTAAATAAAGTATGCCTAAGTACAGTTCCGacaccacc                                                                                  >  2:612170/1‑149 (MQ=255)
                                                             ccTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacc                                                                                  <  1:1458428/149‑1 (MQ=255)
                                                                ctTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAACTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACt                                                                               <  2:190668/149‑1 (MQ=255)
                                                                 ttCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTc                                                                              >  2:725069/1‑149 (MQ=255)
                                                                   cACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCAt                                                                            >  2:985609/1‑149 (MQ=255)
                                                                    aCTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                           <  1:1022447/149‑1 (MQ=255)
                                                                    aCTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAAAAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                           <  1:149648/149‑1 (MQ=255)
                                                                      tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGTCAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGt                                                                         <  2:1087882/149‑1 (MQ=255)
                                                                      tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGt                                                                         <  1:418962/149‑1 (MQ=255)
                                                                         cGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGt                                                                      <  1:1180237/149‑1 (MQ=255)
                                                                         cGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGt                                                                      <  2:123642/149‑1 (MQ=255)
                                                                           aaaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa                                                                    <  1:159951/149‑1 (MQ=255)
                                                                           aaaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa                                                                    <  1:619995/149‑1 (MQ=255)
                                                                           aaaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGCGACGGCACGCACCGTCCCATTTACGAGAAATACACAGGGAGTAAAAAAAGTATGCCAAAGAACCGTTCCGACACCACCACTCATGGACGTaa                                                                    >  1:1588446/1‑149 (MQ=255)
                                                                              ccGAAGATAACTGGGGCTGGTTAGATCATGTTCATCTATCACTACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAatat                                                                 <  2:234887/149‑1 (MQ=255)
                                                                               cGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATg                                                                >  1:115199/1‑149 (MQ=255)
                                                                                 aaGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGc                                                              <  2:216699/149‑1 (MQ=255)
                                                                                   gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCgg                                                            >  1:1390833/1‑149 (MQ=255)
                                                                                      aaaTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTg                                                         <  2:1326393/149‑1 (MQ=255)
                                                                                         tGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAATTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTc                                                      <  1:1364973/149‑1 (MQ=255)
                                                                                            ggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAGGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTg                                                   >  2:1179489/1‑149 (MQ=255)
                                                                                                 gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtg                                              >  2:1366501/1‑149 (MQ=255)
                                                                                                 gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtg                                              >  2:783992/1‑149 (MQ=255)
                                                                                                     gATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgc                                          >  2:929370/1‑149 (MQ=255)
                                                                                                       tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgcgc                                        <  1:612170/149‑1 (MQ=255)
                                                                                                        cAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCCCCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgcgcg                                       >  1:682670/1‑149 (MQ=255)
                                                                                                           gatAATCAATAAATACAAAAAAAGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCa                                    <  1:235259/147‑1 (MQ=255)
                                                                                                           gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTACGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGccc                                    >  2:1510481/1‑148 (MQ=255)
                                                                                                            ttAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCAc                                   >  2:1246431/1‑149 (MQ=255)
                                                                                                              aaTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCg                                 >  1:1073752/1‑149 (MQ=255)
                                                                                                                tCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCTCCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGa                               <  2:446474/149‑1 (MQ=255)
                                                                                                                tCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCCCCCCTCATGGTCGTAATTTGGCGGGTGCTCGTGCGCTGTGGCTCGCCACCGGa                               >  2:1001816/1‑149 (MQ=255)
                                                                                                                 cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCTCCGGaa                              <  1:592166/149‑1 (MQ=255)
                                                                                                                 cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCTCCGGaa                              <  1:597826/149‑1 (MQ=255)
                                                                                                                  aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTTATATGGCGGGTGCTCGTGCGCTNTGGCGCGCCACCngnnn                             <  2:27839/149‑4 (MQ=255)
                                                                                                                  aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAAt                             <  1:257288/149‑1 (MQ=255)
                                                                                                                  aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGACACCAACACTCATTGTCGAAATATGGCGGGTGCTCGTGCTCTGTGGCGCGCCACCGGAAt                             >  2:997460/1‑149 (MQ=255)
                                                                                                                      aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGCCACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGAcc                         >  1:1162239/1‑149 (MQ=255)
                                                                                                                          tcgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACg                     >  2:996563/5‑149 (MQ=255)
                                                                                                                            aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGcc                   >  2:306431/1‑149 (MQ=255)
                                                                                                                            aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGcc                   >  2:1146448/1‑149 (MQ=255)
                                                                                                                             taaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCg                  >  2:712087/2‑149 (MQ=255)
                                                                                                                              aaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCATAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGa                 <  2:1154539/149‑1 (MQ=255)
                                                                                                                               aaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGAt                <  1:1144938/149‑1 (MQ=255)
                                                                                                                               aaTGGGACCGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGAt                <  1:412772/149‑1 (MQ=255)
                                                                                                                                aTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGCCGCCGAtt               >  2:160975/1‑149 (MQ=255)
                                                                                                                                 tGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGAttt              <  1:101144/149‑1 (MQ=255)
                                                                                                                                          acgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGcc     >  2:1179873/1‑149 (MQ=255)
                                                                                                                                             cacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt  <  2:115199/149‑1 (MQ=255)
                                                                                                                                             cacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt  <  1:1324320/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAT  >  NC_000913/3953295‑3953584

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: