Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 262,225 G→T 10.5% V241V (GTG→GTT proA → gamma‑glutamylphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913262,2250GT10.5% 243.9 / 7.2 98V241V (GTG→GTTproAgamma‑glutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (37/47);  new base T (10/2);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGCAGCTGC  >  NC_000913/262078‑262371
                                                                                                                                                   |                                                                                                                                                  
aTACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGAt                                                                                                                                                  <  2:741921/150‑1 (MQ=255)
 tACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATc                                                                                                                                                 <  1:762143/150‑1 (MQ=255)
 tACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATAATTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATc                                                                                                                                                 >  2:599205/1‑150 (MQ=255)
   cATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAACCTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCga                                                                                                                                               >  2:1053840/1‑149 (MQ=255)
   cATCGACATGCTGATCCCGAGTGGAGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAAGCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTTAAAGTGATCGt                                                                                                                                               >  1:844606/1‑150 (MQ=255)
    aTCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGACGTGAACAGTCGACACTCCCGGTGATAACAGGTGGTATAGTCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTTATCTTc                                                                                                                                              >  2:1619621/1‑150 (MQ=255)
      cGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCaa                                                                                                                                            <  2:745010/150‑1 (MQ=255)
         cATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAAcgc                                                                                                                                         <  2:1694837/150‑1 (MQ=255)
         cATGCTGATCCCGCGTGGGGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTCCGTTGATGAAAGGGTAGAGATCGCTGAAGCATTAAAAGTGCTCGTCAAcgc                                                                                                                                         >  1:1044100/1‑150 (MQ=255)
          aTGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAAcgcg                                                                                                                                        <  2:619081/150‑1 (MQ=255)
           tGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGAAGAGATCGCTGAAGCAATAAAAGTGATCGTCAACGCGa                                                                                                                                       >  2:1629719/1‑150 (MQ=255)
           tGCTGATCCCGCGTGGTGGCGCTGGTTTGAATAAACTGTGCCGTGAACAGTCGACAATCCCGGTCAACACAGGTGGTATAGGGCTATGCAATATTTACGTTGATTAAAGTGTTGAGATCGCTGAAGCATTAAAAGTGAGCGTCAACGCGa                                                                                                                                       >  1:165761/1‑150 (MQ=255)
              tGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGaaaa                                                                                                                                    >  1:156332/1‑150 (MQ=255)
                  cccGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCa                                                                                                                                <  1:555094/150‑1 (MQ=255)
                  cccGCGTGGTGGCGCTGGTTTGCATAAACTATGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCa                                                                                                                                >  1:909611/1‑150 (MQ=255)
                   ccGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAg                                                                                                                               <  1:477865/150‑1 (MQ=255)
                   ccGCGTGGGGGCGCTGGGTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAg                                                                                                                               <  1:1441724/150‑1 (MQ=255)
                     gcgTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCg                                                                                                                             <  2:1685343/150‑1 (MQ=255)
                     gcgTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCACTGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCg                                                                                                                             <  1:1698042/150‑1 (MQ=255)
                                cTGGTTTGCATAAACTGTGCCGTTAAAAGTCGACAATCCCTGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACAt                                                                                                                  <  1:226941/150‑1 (MQ=255)
                                       gCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATAc                                                                                                           <  2:972378/150‑1 (MQ=255)
                                         aTAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACgg                                                                                                         <  1:636014/150‑1 (MQ=255)
                                          tAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGt                                                                                                        >  1:751742/1‑150 (MQ=255)
                                          tAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGt                                                                                                        <  2:1207522/150‑1 (MQ=255)
                                          tAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGt                                                                                                        >  2:690836/1‑150 (MQ=255)
                                          tAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGt                                                                                                        <  2:144013/150‑1 (MQ=255)
                                           aaaCTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGtt                                                                                                       <  2:1099463/150‑1 (MQ=255)
                                            aaCTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTg                                                                                                      <  1:48137/150‑1 (MQ=255)
                                            aaCTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTg                                                                                                      >  2:507793/1‑150 (MQ=255)
                                                 tgCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACg                                                                                                 <  2:897812/150‑1 (MQ=255)
                                                 tgCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTCATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACg                                                                                                 <  2:1249810/150‑1 (MQ=255)
                                                   ccGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGGCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGtt                                                                                               <  2:1196631/150‑1 (MQ=255)
                                                    cGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTg                                                                                              <  2:788099/150‑1 (MQ=255)
                                                          cAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTg                                                                                        <  2:1340420/150‑1 (MQ=255)
                                                          cAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTGGAAACGTTGCTGGTg                                                                                        <  2:1069660/150‑1 (MQ=255)
                                                          cAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGAAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCCTCCGAGCACATGTAATACGGTTGAAACGATGCTGGTg                                                                                        >  1:1561309/1‑150 (MQ=255)
                                                            gTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGaa                                                                                      <  1:48598/150‑1 (MQ=255)
                                                             tCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAAt                                                                                     <  2:1436364/150‑1 (MQ=255)
                                                              cGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATa                                                                                    >  1:1594678/1‑150 (MQ=255)
                                                              cGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGAACGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATa                                                                                    >  2:924549/1‑150 (MQ=255)
                                                              cGACAATCCCCGTGATCACAGCTGGTAAAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTCAAGCATTAAAAGTCATCGTCAACGCGAAAACTCACCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATa                                                                                    >  1:1731545/1‑150 (MQ=255)
                                                               gACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATaa                                                                                   >  1:1161615/1‑150 (MQ=255)
                                                                 cAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGATATCGCTGAAACAATAAAAATTATAATAAAATCTAAAACGAAGCGTCCGAGCACAGGTAATACGGTTGACACTTTGCTGGTGAATaaaa                                                                                 >  2:726669/1‑150 (MQ=255)
                                                                 cAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATaaaa                                                                                 <  2:631759/150‑1 (MQ=255)
                                                                 cAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATaaaa                                                                                 <  2:57677/150‑1 (MQ=255)
                                                                   aTCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGCAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCAGAGCACATATAATACGGTTGAAACGTTGCGGGTGAATaaaaaa                                                                               >  1:222123/1‑149 (MQ=255)
                                                                   aTCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAAc                                                                               >  1:220185/1‑150 (MQ=255)
                                                                   aTCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGATGGTGAATaaaaaa                                                                               >  1:1110295/1‑149 (MQ=255)
                                                                    tCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACa                                                                              <  2:909611/150‑1 (MQ=255)
                                                                    tCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACa                                                                              >  2:1172032/1‑150 (MQ=255)
                                                                     cccGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGATATCGCTGAACCAATTAAAATTAAAAAAAAAAAGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACAt                                                                             >  2:1003216/1‑150 (MQ=255)
                                                                     cccGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAATAGTTATCGTCAACGCGAAAACTCAGCTTACGAGAACATGTAATACGATTGAAACGTTGCTGGTGAATAAAAACAt                                                                             >  2:855703/1‑150 (MQ=255)
                                                                     cccGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACAt                                                                             <  2:1278098/150‑1 (MQ=255)
                                                                       cGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCg                                                                           <  2:944020/150‑1 (MQ=255)
                                                                           gATCACAGGTGGTTTAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGa                                                                       >  1:1051199/1‑150 (MQ=255)
                                                                            aTCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAACCAATAAAAATTATCTACAACACGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGAt                                                                      >  2:1344853/1‑150 (MQ=255)
                                                                            aTCACAGGGGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTTATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGAt                                                                      <  2:1236216/150‑1 (MQ=255)
                                                                                caGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGATATCGCTGAACCATTAAAAGTGAACTACAACTCGAAAACTCAGAGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCt                                                                  >  2:457508/1‑150 (MQ=255)
                                                                                  ggtggtATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTc                                                                >  2:1467754/1‑150 (MQ=255)
                                                                                   gtggtATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTcc                                                               <  1:690836/150‑1 (MQ=255)
                                                                                       tataGGCGTATGCCATATTTACGTTGATGAAAGTGTAGATATCGCTGAACCAATAAAAATTATAAAAAAAACGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCTCCGATAGCTTCCTGcc                                                           >  2:288644/1‑150 (MQ=255)
                                                                                       tataGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGcc                                                           <  1:413729/150‑1 (MQ=255)
                                                                                       tataGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGcc                                                           >  2:404811/1‑150 (MQ=255)
                                                                                         taGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCg                                                         <  2:413156/150‑1 (MQ=255)
                                                                                          aGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCGGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGc                                                        <  1:1499491/150‑1 (MQ=255)
                                                                                           ggCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCa                                                       >  1:892816/1‑150 (MQ=255)
                                                                                           ggCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCa                                                       >  1:850613/1‑150 (MQ=255)
                                                                                              ttATGCCATATTTACGTTGATGAAAGTGTATAGATGGCTGAAGCATTATAAGTGTTCGTCAACGCGAAAACTCAGCGTCCGAGAACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTa                                                    <  2:165186/149‑1 (MQ=255)
                                                                                                aTGCCATATTTACGATGATGTAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACTCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTTGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAg                                                  <  2:1395692/150‑1 (MQ=255)
                                                                                                 tGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGc                                                 <  1:513223/150‑1 (MQ=255)
                                                                                                  gCCATATTTACCTTGATGAATGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACTTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGca                                                <  1:687863/150‑1 (MQ=255)
                                                                                                   ccATATTTACGTTGATGAAAGTGTAGGTATCGCTGAAGCAATTAAAGTTAAAGTCAAAGAGAAAACTCAGAGTCCTAGAAAATGTAATACGGTTGAAACTATGGTGGCGAATAAAAATATCTCAGATAGCTTCCTGCCCGCATTAAGcaa                                               >  2:721022/1‑150 (MQ=255)
                                                                                                   ccATATTTACGTTGATGAAAGTGTAGATATCGCTGAAGCAATAAACTTATAATCAAAGAGAAAACTCAGCTTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCAGTAAGcaaa                                              >  2:1047357/1‑150 (MQ=255)
                                                                                                    cATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGcaaa                                              >  1:1274874/1‑150 (MQ=255)
                                                                                                           aCGTTGATGAATGTGTAGAGAGCGCTGAAGCATTAAAAGTGTGCGTCAACGCGAAAACTCAGCGGCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATgg                                       <  2:913294/150‑1 (MQ=255)
                                                                                                              ttGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCgg                                    <  2:992220/150‑1 (MQ=255)
                                                                                                                 atgaAAGTGTAGATATCGCTGAACCAATAAAAATTATAAAAAAAACGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCGTTAAGCAAACAAATGGCGGaac                                 >  2:238209/1‑149 (MQ=255)
                                                                                                                   gaAAGTGTAGAGCTCGCTGAAGCAATAAAAATTATAATAAAAGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACCTCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAgc                               >  2:96710/1‑150 (MQ=255)
                                                                                                                         gtAGAGATCGCTTTAGCATTAAAAGTTATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTg                         <  2:536630/150‑1 (MQ=255)
                                                                                                                          tAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGa                        >  2:630653/1‑150 (MQ=255)
                                                                                                                          tAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGa                        >  1:1342679/1‑150 (MQ=255)
                                                                                                                            gagaTCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACTTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACa                      <  1:1396110/150‑1 (MQ=255)
                                                                                                                             agaTCGCTGAAGCATTAAAAGGGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACAt                     <  2:207603/150‑1 (MQ=255)
                                                                                                                                tCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTac                  <  1:448704/150‑1 (MQ=255)
                                                                                                                                 cGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTaca                 >  1:810270/1‑150 (MQ=255)
                                                                                                                                  gCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTacac                >  1:809483/1‑150 (MQ=255)
                                                                                                                                  gCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTacac                >  1:906178/1‑150 (MQ=255)
                                                                                                                                   cTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCACGCATTAATCAAACAAATGGCGGAAAGCGGCGTGACATAACACg               >  1:297075/1‑150 (MQ=255)
                                                                                                                                    tGAAGCATTAAAAGTGATCGTCAACGCGAAACCTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGc              <  1:1053840/150‑1 (MQ=255)
                                                                                                                                     taaGCATTAAAAGTGTTCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCa             <  1:1468483/149‑1 (MQ=255)
                                                                                                                                     gAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCa             <  2:1428768/150‑1 (MQ=255)
                                                                                                                                      aaGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAg            <  2:220185/150‑1 (MQ=255)
                                                                                                                                          aTTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGc        >  1:630661/1‑150 (MQ=255)
                                                                                                                                          aTTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGc        >  2:588472/1‑150 (MQ=255)
                                                                                                                                           ttAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGCa       <  2:1371442/150‑1 (MQ=255)
                                                                                                                                             aaaaGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGCAGc     >  2:1403158/1‑150 (MQ=255)
                                                                                                                                              aaaGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGCAGCt    >  1:1241275/1‑150 (MQ=255)
                                                                                                                                                aGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGCAGCTGc  >  1:339356/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATACATCGACATGCTGATCCCGCGTGGTGGCGCTGGTTTGCATAAACTGTGCCGTGAACAGTCGACAATCCCGGTGATCACAGGTGGTATAGGCGTATGCCATATTTACGTTGATGAAAGTGTAGAGATCGCTGAAGCATTAAAAGTGATCGTCAACGCGAAAACTCAGCGTCCGAGCACATGTAATACGGTTGAAACGTTGCTGGTGAATAAAAACATCGCCGATAGCTTCCTGCCCGCATTAAGCAAACAAATGGCGGAAAGCGGCGTGACATTACACGCAGATGCAGCTGC  >  NC_000913/262078‑262371

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