Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 959,667 A→C 5.3% T286P (ACC→CCC)  aroA → 5‑enolpyruvylshikimate‑3‑phosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913959,6670AC 5.3% 340.5 / 6.6 112T286P (ACC→CCC) aroA5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base A (49/57);  new base C (3/3);  total (52/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 114 total aligned reads displayed.

CTTATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATC  >  NC_000913/959518‑959813
                                                                                                                                                     |                                                                                                                                                  
cTTATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGa                                                                                                                                                    >  1:1069876/1‑150 (MQ=255)
cTTATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAACTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGa                                                                                                                                                    <  1:282897/150‑1 (MQ=255)
   aTTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCa                                                                                                                                                 >  1:1734570/1‑150 (MQ=255)
   aTTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCa                                                                                                                                                 >  1:681175/1‑150 (MQ=255)
    tttGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCAt                                                                                                                                                >  1:1097164/1‑150 (MQ=255)
     ttGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCAtt                                                                                                                                               <  2:1513442/150‑1 (MQ=255)
       ggTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTg                                                                                                                                             >  2:1132394/1‑150 (MQ=255)
        gTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGc                                                                                                                                            <  1:728600/150‑1 (MQ=255)
          cGAAGGCGTAGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCACAGGCGGCACTGTATAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGTCGCGCCCAttttttc                                                                                                                                          <  1:930702/150‑5 (MQ=255)
          cGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTg                                                                                                                                          <  1:506952/150‑1 (MQ=255)
           gAAGGCGATGCATCTTCGGCTTCTTACTTTCTGTCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGTCGCGCCCAttttttcc                                                                                                                                         <  1:396816/150‑6 (MQ=255)
            aatgCGATGCATCTTCGGCTTATTACTTTCTGGCAGCAGCAGAAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTTTGCATGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTggg                                                                                                                                        <  2:1792266/147‑1 (MQ=255)
              ggCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGc                                                                                                                                      >  1:500125/1‑150 (MQ=255)
                cGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCga                                                                                                                                    <  1:1456166/150‑1 (MQ=255)
                    gCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCgatgat                                                                                                                                <  2:965666/150‑1 (MQ=255)
                       tCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATtat                                                                                                                             <  1:1586077/150‑1 (MQ=255)
                        cTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATtata                                                                                                                            >  2:55462/1‑150 (MQ=255)
                         ttCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGTACGTAACAGTATGCAGGGTGATATTAGCTTTGCTGATGTGCTGGAAAAAATGGGAGCGACCATTTGCTGGGGCGATGATTatat                                                                                                                           >  2:414110/1‑150 (MQ=255)
                           cGGCTTCTTACTTTCTGGCAGCATCAGCACTCAAAGGCGGCACAGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTCCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATAttt                                                                                                                         >  2:706295/1‑150 (MQ=255)
                             gCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTcc                                                                                                                       <  1:1138131/150‑1 (MQ=255)
                             gCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTcc                                                                                                                       <  1:1020309/150‑1 (MQ=255)
                              cttcttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGGGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCt                                                                                                                      <  1:1025058/150‑1 (MQ=255)
                              attcatAGTTTCTGGCAGCAGCAGCAATCAAAGGCGTCACTGTAACAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATAGTCGCTTTGCTGATGTGCTTGAAAAAATGGGCGCGACCAATTGCTGGGGCGATGATTATATTTCCt                                                                                                                      <  2:310293/149‑1 (MQ=255)
                                tcttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGc                                                                                                                    <  2:1345504/150‑1 (MQ=255)
                                tcttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGc                                                                                                                    <  2:1352894/150‑1 (MQ=255)
                                 cttACTTTCTGTCATCAGCAGCAATCAAAGGCTGAACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTTATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCa                                                                                                                   <  1:1169200/150‑1 (MQ=255)
                                   tACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGCTGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCAcg                                                                                                                 >  2:539090/1‑150 (MQ=255)
                                   tACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCAcg                                                                                                                 <  2:1065194/150‑1 (MQ=255)
                                         cTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggt                                                                                                           >  1:612663/1‑150 (MQ=255)
                                         cTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggt                                                                                                           <  2:1548989/150‑1 (MQ=255)
                                          tGGCCGCCGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGCCGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggtg                                                                                                          <  1:334582/150‑1 (MQ=255)
                                          tGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggtg                                                                                                          <  1:1087921/150‑1 (MQ=255)
                                            gcagcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaa                                                                                                        >  2:1794537/1‑150 (MQ=255)
                                             cagcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGTTGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaac                                                                                                       <  1:602572/150‑1 (MQ=255)
                                              agcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTCCACGCGTGGtgaact                                                                                                      <  2:1240442/150‑1 (MQ=255)
                                               gcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaactg                                                                                                     <  2:376786/150‑1 (MQ=255)
                                                cagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTCACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaactga                                                                                                    >  1:255474/1‑150 (MQ=255)
                                                   cagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACg                                                                                                 <  2:674220/150‑1 (MQ=255)
                                                     gcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCt                                                                                               <  2:1440583/150‑1 (MQ=255)
                                                     gcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCt                                                                                               >  1:769993/1‑150 (MQ=255)
                                                      caATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTa                                                                                              >  2:62140/1‑150 (MQ=255)
                                                         tCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTg                                                                                           <  2:1700520/150‑1 (MQ=255)
                                                          cAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTga                                                                                          <  1:1184092/150‑1 (MQ=255)
                                                          cAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACCTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTga                                                                                          <  2:691878/150‑1 (MQ=255)
                                                            aaGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTgata                                                                                        >  2:1624705/1‑150 (MQ=255)
                                                               gcggcACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTgatatgg                                                                                     >  1:1823360/1‑150 (MQ=255)
                                                                gggcACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTgatatgga                                                                                    <  1:1466063/149‑1 (MQ=255)
                                                                    aCTGTAAAAGTGACAGGTATTGGACGTAACAGTAAGCATGGTGAAATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTgatatggatatg                                                                                <  1:1262762/150‑1 (MQ=255)
                                                                     cTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAGCTGAACGCTATTGATATGGATATGa                                                                               >  2:86053/1‑150 (MQ=255)
                                                                        tAAAAGTGACCGGTATTGGACTTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAAcc                                                                            <  1:888149/150‑1 (MQ=255)
                                                                        tAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAAcc                                                                            >  1:1691839/1‑150 (MQ=255)
                                                                           aaGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCata                                                                         <  2:714913/150‑1 (MQ=255)
                                                                             gTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATAtt                                                                       >  1:6531/1‑150 (MQ=255)
                                                                               gACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTcc                                                                     <  2:1734570/150‑1 (MQ=255)
                                                                               gACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTcc                                                                     <  2:448156/150‑1 (MQ=255)
                                                                                aCCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCt                                                                    >  2:594417/1‑150 (MQ=255)
                                                                                      aTTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATgcg                                                              <  1:1132394/150‑1 (MQ=255)
                                                                                        tGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGATAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGGTATTGATATGGATATGAACCATATTCCTGCTgcggc                                                            >  1:1183455/1‑150 (MQ=255)
                                                                                        tGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATgcggc                                                            <  1:1624705/150‑1 (MQ=255)
                                                                                        tGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATgcggc                                                            >  2:1552864/1‑150 (MQ=255)
                                                                                         ggACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATgcggcg                                                           <  1:478128/150‑1 (MQ=255)
                                                                                          gACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGATCCATATTCCTGATGCGGCGa                                                          >  2:618024/1‑150 (MQ=255)
                                                                                            cGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATg                                                        <  1:1264496/150‑1 (MQ=255)
                                                                                            cGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCAAGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATg                                                        >  1:633748/1‑150 (MQ=255)
                                                                                              tAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAc                                                      >  1:596942/1‑150 (MQ=255)
                                                                                              tAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAc                                                      >  1:1396219/1‑150 (MQ=255)
                                                                                              tAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAc                                                      >  1:983008/1‑150 (MQ=255)
                                                                                              tAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAc                                                      >  2:1041253/1‑150 (MQ=255)
                                                                                               aaCAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAcc                                                     >  1:516932/1‑150 (MQ=255)
                                                                                               aaCAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTAGCTGGGGCGATGATTATATTTCCTGCACGCCTGGTGATCTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAcc                                                     >  1:99824/1‑150 (MQ=255)
                                                                                                   gTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCAttt                                                 >  1:215078/1‑149 (MQ=255)
                                                                                                   gTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTg                                                 <  2:1726200/150‑1 (MQ=255)
                                                                                                   gTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTg                                                 <  2:548818/150‑1 (MQ=255)
                                                                                                    tATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGc                                                <  1:1794671/150‑1 (MQ=255)
                                                                                                    tATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGc                                                <  2:1743510/150‑1 (MQ=255)
                                                                                                    tATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGc                                                <  2:462384/150‑1 (MQ=255)
                                                                                                     aTGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGcc                                               >  1:149490/1‑150 (MQ=255)
                                                                                                        cAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCAcg                                            >  2:1375362/1‑150 (MQ=255)
                                                                                                              gaaaatCGCTTTTGTGATGTGCTGGAAAAAATTGGCGCGACCATTTGCCGGGGCGATGATTATATTTCCTGCACGCGTGGGGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGgcggcg                                      <  1:1229039/145‑1 (MQ=255)
                                                                                                              gATATTCGCTTTGCGGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGAGGCGATGACCATTGCCACGgcggcg                                      >  2:262272/1‑150 (MQ=255)
                                                                                                               atatTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGCCCATTTGCTGGGGCGATGATTATATTTCCTGCACTCGAGGTGATCTGAACGCTATTGATATGGATATGAACAATATTCGTGATGCGGCGATGACCATTGCCACggcggcct                                     >  1:286599/1‑148 (MQ=255)
                                                                                                               atatTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGt                                     <  2:1551984/150‑1 (MQ=255)
                                                                                                               atatTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGt                                     <  1:1052425/150‑1 (MQ=255)
                                                                                                                tatTCGCTTTGCTGATGGGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGtt                                    <  2:1684445/150‑1 (MQ=255)
                                                                                                                 atTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTa                                   >  1:1400035/1‑150 (MQ=255)
                                                                                                                    cGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTAttt                                <  1:1041253/150‑1 (MQ=255)
                                                                                                                    cGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCAACGGCGGCGTTAttt                                >  1:1328965/1‑150 (MQ=255)
                                                                                                                     gCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTg                               <  2:500125/150‑1 (MQ=255)
                                                                                                                                 tgCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTACTTGCAAAAGGcacca                   >  1:1550837/1‑150 (MQ=255)
                                                                                                                                  gCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccac                  >  1:898668/1‑150 (MQ=255)
                                                                                                                                    tGGAAAAAATGGGCGCGACGATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccacca                <  2:681175/150‑1 (MQ=255)
                                                                                                                                    tGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGACCCATAATCCAGATGCGGCGATGCCCATTGCCACGGCGGCGTTATTTGCAAAAGGcacaaccc                >  1:916182/1‑149 (MQ=255)
                                                                                                                                     ggAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGCCCATTGCCCCGGCGGCGTTATTTGCAAAAGGcaccaccac               >  1:1580784/1‑150 (MQ=255)
                                                                                                                                     ggAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccaccac               >  2:401890/1‑150 (MQ=255)
                                                                                                                                     ggAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccaccac               >  1:147853/1‑150 (MQ=255)
                                                                                                                                      cagcaaaTGCTCGCGCCCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg              >  1:1632648/5‑150 (MQ=255)
                                                                                                                                        aaaaaTGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCt            <  2:571346/150‑1 (MQ=255)
                                                                                                                                                  gcgACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATc  <  2:225810/150‑1 (MQ=255)
                                                                                                                                                  gcgACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATc  >  2:1420576/1‑150 (MQ=255)
                                                                                                                                                  gcgACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATc  <  1:55462/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTTATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATC  >  NC_000913/959518‑959813

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: