Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,335,655 T→A 5.3% intergenic (‑65/+39) eptA ← / ← adiC lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,335,6550TA 5.3% 1162.2 / 37.4 338intergenic (‑65/+39)eptA/adiClipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (143/177);  new base A (10/8);  total (153/185)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 345 total aligned reads displayed.

CAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAA  >  NC_000913/4335507‑4335803
                                                                                                                                                      |                                                                                                                                                    
 aGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTc                                                                                                                                                     <  3:2389153/149‑1 (MQ=255)
   gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCtt                                                                                                                                                   >  2:248528/1‑149 (MQ=255)
     aGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTaa                                                                                                                                                 <  3:82837/149‑1 (MQ=255)
      gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAg                                                                                                                                                <  2:683131/149‑1 (MQ=255)
       aTAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                               >  3:1247567/1‑149 (MQ=255)
         agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                             <  2:1878642/149‑1 (MQ=255)
         agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                             <  4:2015895/149‑1 (MQ=255)
            gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAACATAAAAGCCAACCTTAAGAACTTAAGGTTg                                                                                                                                            <  2:848688/149‑1 (MQ=38)
             atatAAAAGGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                                                                         >  1:1132658/1‑149 (MQ=255)
                   aaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCtttt                                                                                                                                      >  2:139633/1‑146 (MQ=255)
                    aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGt                                                                                                                                  >  3:1799415/1‑149 (MQ=255)
                     aGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGtt                                                                                                                                 <  3:515356/149‑1 (MQ=255)
                      gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGttt                                                                                                                                <  4:237390/149‑1 (MQ=255)
                        ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                              <  4:1862254/149‑1 (MQ=255)
                        ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGATGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                              >  3:767758/1‑149 (MQ=255)
                           caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTgg                                                                                                                           >  1:1512879/1‑149 (MQ=255)
                            aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGa                                                                                                                          >  1:2226320/1‑149 (MQ=255)
                              caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGAtt                                                                                                                        <  2:303544/149‑1 (MQ=255)
                                 tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                     <  3:810958/149‑1 (MQ=255)
                                  aGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgg                                                                                                                    <  1:2254364/149‑1 (MQ=255)
                                  aGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgg                                                                                                                    <  3:714431/149‑1 (MQ=255)
                                    ccAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgggg                                                                                                                  <  1:1403271/149‑1 (MQ=255)
                                      aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                                <  1:396143/149‑1 (MQ=255)
                                      aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                                >  1:1774928/1‑149 (MQ=255)
                                          gAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTaa                                                                                                            <  2:2001838/149‑1 (MQ=255)
                                            gTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATc                                                                                                          >  3:1321449/1‑149 (MQ=255)
                                             tAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCt                                                                                                         >  2:1200152/1‑149 (MQ=255)
                                               aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCttt                                                                                                       >  3:1572378/1‑149 (MQ=255)
                                                aTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                      >  4:2112380/1‑149 (MQ=255)
                                                 ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                     >  4:1950406/1‑149 (MQ=255)
                                                 ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCACAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                     >  2:1697801/1‑149 (MQ=255)
                                                    aaaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                     >  2:1453380/1‑146 (MQ=255)
                                                     aaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                                 <  2:2411338/149‑1 (MQ=255)
                                                      agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                                >  1:1963584/1‑149 (MQ=255)
                                                       gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGAt                                                                                                                         >  3:2004874/1‑123 (MQ=255)
                                                       gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                               <  3:1134754/149‑1 (MQ=255)
                                                         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                             <  3:602225/149‑1 (MQ=255)
                                                         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                             >  3:918012/1‑149 (MQ=255)
                                                                 ttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGTATCTaa                                                                                     <  1:1960252/149‑1 (MQ=255)
                                                                  tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                                    <  1:786568/149‑1 (MQ=255)
                                                                   agtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                                                                   <  4:792974/149‑1 (MQ=255)
                                                                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCGTATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                                 >  3:2230840/1‑134 (MQ=255)
                                                                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCGTATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                                 <  4:2230840/134‑1 (MQ=255)
                                                                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAggg                                                                                  <  1:2016590/149‑1 (MQ=255)
                                                                    gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAggg                                                                                  >  3:1464840/1‑149 (MQ=255)
                                                                        gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                     >  1:2402615/1‑126 (MQ=255)
                                                                        gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATAc                                                                              <  3:1673987/149‑1 (MQ=255)
                                                                          gcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                                                                            >  1:465178/1‑149 (MQ=255)
                                                                            ttCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATc                                                                                                          <  1:1206854/117‑1 (MQ=255)
                                                                              cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttt                                                                        <  2:1223780/149‑1 (MQ=255)
                                                                                 cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAGGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                     <  4:848449/149‑1 (MQ=255)
                                                                                 cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                     <  3:1038679/149‑1 (MQ=255)
                                                                                 cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                     <  4:1149713/149‑1 (MQ=255)
                                                                                  aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                    <  2:1963584/149‑1 (MQ=255)
                                                                                  aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                    <  2:2226169/149‑1 (MQ=255)
                                                                                     aCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa                                                                 <  4:1510082/149‑1 (MQ=255)
                                                                                       ggTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc                                                               <  1:2115270/149‑1 (MQ=255)
                                                                                           tttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                                <  4:1031915/96‑1 (MQ=255)
                                                                                            ttCCATCGAACAAAGTGCGCATATGTTCGCAAATCAAA‑AT‑TAAGCCAACCTTAAGTTCCTAAGGATGGCTTTTATGTTTGCTGGATTAAGGGAAATTAATCTGTGCGTATTGGTGCATCTAATGGATACAGGTTTTTATGTAGCCGGtt                                                          <  4:1676647/149‑1 (MQ=255)
                                                                                            ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGtt                                                          <  1:1473822/149‑1 (MQ=255)
                                                                                             tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                         <  2:621332/149‑1 (MQ=255)
                                                                                               cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTa                                                       >  1:622197/1‑149 (MQ=255)
                                                                                                aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                      >  3:2119750/1‑149 (MQ=255)
                                                                                                aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                      >  4:2291223/1‑149 (MQ=255)
                                                                                                  cGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                                    >  4:338018/1‑149 (MQ=255)
                                                                                                    aaCAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                               <  2:809952/104‑1 (MQ=255)
                                                                                                    aaCAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                               >  1:809952/1‑104 (MQ=255)
                                                                                                    aaCAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCa                                                  <  1:248528/149‑1 (MQ=255)
                                                                                                    aaCAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCa                                                  <  4:1321449/149‑1 (MQ=255)
                                                                                                    aaCAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGACAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTTGATTAAGGAGAATTACTCTTTGGTTATTGGTGCATCTAACGGATACGGGTTTTAATGTAGCCGGTTGTAATCCa                                                  <  4:1922350/149‑1 (MQ=255)
                                                                                                      cAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgc                                                <  3:338018/149‑1 (MQ=255)
                                                                                                       aaaGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgcg                                               >  1:354967/1‑149 (MQ=255)
                                                                                                          gTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCa                                                  <  4:1087071/143‑1 (MQ=255)
                                                                                                          gTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCAt                                            >  4:1374654/1‑149 (MQ=255)
                                                                                                              gcATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                    >  1:1689748/1‑121 (MQ=255)
                                                                                                              gcATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACAt                                        <  1:1370821/149‑1 (MQ=255)
                                                                                                                atatGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATgg                                      >  3:399667/1‑149 (MQ=255)
                                                                                                                      tCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                                <  3:1391866/149‑1 (MQ=255)
                                                                                                                      tCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                                <  2:465178/149‑1 (MQ=255)
                                                                                                                      tCGCAAAGCAAA‑AT‑TAAGCAAACCTTAAGTTCTTAAGGTTGGCTATTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                                                                           >  4:544076/1‑90 (MQ=255)
                                                                                                                        gcaaagcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                              >  4:690773/1‑149 (MQ=255)
                                                                                                                         caaagcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                             <  2:1783898/149‑1 (MQ=255)
                                                                                                                         caaagcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                             >  1:1896312/1‑149 (MQ=255)
                                                                                                                           aagcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAcc                           <  1:1107075/149‑1 (MQ=255)
                                                                                                                             gcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                      >  2:1331366/1‑72 (MQ=255)
                                                                                                                             gcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                      <  1:1331366/72‑1 (MQ=255)
                                                                                                                                    taaaaGCCAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAggg                    >  2:1123258/4‑149 (MQ=255)
                                                                                                                                    t‑tAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTTTGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCACCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt                   <  3:2037624/149‑1 (MQ=255)
                                                                                                                                    t‑tAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                               >  2:1005124/1‑73 (MQ=255)
                                                                                                                                    t‑tAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                               <  1:1005124/73‑1 (MQ=255)
                                                                                                                                    t‑tAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                         >  4:162432/1‑111 (MQ=255)
                                                                                                                                         gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa                                                                 <  3:1735643/99‑1 (MQ=255)
                                                                                                                                         gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa                                                                 >  4:1735643/1‑99 (MQ=255)
                                                                                                                                           cAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa             >  3:163336/1‑149 (MQ=255)
                                                                                                                                            aaCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAAt            >  2:903627/1‑149 (MQ=255)
                                                                                                                                             aCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                                     <  4:1396179/75‑1 (MQ=255)
                                                                                                                                               cTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAc         <  2:1502436/149‑1 (MQ=255)
                                                                                                                                                ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                                     >  1:866886/1‑72 (MQ=255)
                                                                                                                                                     ttCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATa   >  3:1913276/1‑149 (MQ=255)
                                                                                                                                                     ttCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATa   >  3:19363/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAA  >  NC_000913/4335507‑4335803

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑