Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,308,192 A→T 6.9% F96Y (TTC→TAC) ‡ alsE ← D‑allulose‑6‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,308,1920AT 6.9% 271.8 / 7.8 87F96Y (TTC→TAC) ‡alsED‑allulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base T (2/4);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAAC  >  NC_000913/4308043‑4308341
                                                                                                                                                     |                                                                                                                                                     
aCGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGa                                                                                                                                                       <  1:1144041/150‑1 (MQ=255)
   gTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACg                                                                                                                                                    <  2:586630/150‑1 (MQ=255)
   gTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACg                                                                                                                                                    <  2:138312/150‑1 (MQ=255)
     tGTCCGGCAAAGCCGGGATCGACTGTCATGACCGTAATTTTATTGGGCTTATGGATTTAGTATTTCATGGGCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCGACTTTCATGTCATGTCGGCGGATTTCATCAATCAGGCGGAACGcc                                                                                                                                                  <  1:1571281/150‑1 (MQ=255)
     tGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATTGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGcc                                                                                                                                                  <  1:1577218/150‑1 (MQ=255)
      gTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCt                                                                                                                                                 >  1:568294/1‑150 (MQ=255)
       tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGTCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTg                                                                                                                                                <  2:1128886/150‑1 (MQ=255)
       tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTg                                                                                                                                                >  1:1004275/1‑150 (MQ=255)
       tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTg                                                                                                                                                >  1:1010412/1‑150 (MQ=255)
       tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTg                                                                                                                                                <  2:1147934/150‑1 (MQ=255)
       tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTg                                                                                                                                                <  2:802678/150‑1 (MQ=255)
        ccGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCAGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTgg                                                                                                                                               >  1:265207/1‑150 (MQ=255)
         cGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGc                                                                                                                                              <  1:1034775/150‑1 (MQ=255)
         cGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCAATTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGc                                                                                                                                              >  1:277166/1‑150 (MQ=255)
          ggCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTgatt                                                                                                                                             <  1:775760/150‑4 (MQ=255)
            cAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTgattga                                                                                                                                           <  2:1010412/150‑6 (MQ=255)
            cAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTgattga                                                                                                                                           <  2:1004275/150‑6 (MQ=255)
              aaGGCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATTTAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTgattgatg                                                                                                                                         <  2:619283/150‑8 (MQ=255)
              aaGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTg                                                                                                                                         <  2:294703/150‑1 (MQ=255)
                    ggATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGtt                                                                                                                                   <  2:741518/150‑1 (MQ=255)
                        cGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCATTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCg                                                                                                                               <  1:806013/150‑1 (MQ=255)
                        cGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCATCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCg                                                                                                                               >  2:1179749/1‑150 (MQ=255)
                         gACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCgg                                                                                                                              >  2:982303/1‑150 (MQ=255)
                           cAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGAt                                                                                                                            >  1:641496/1‑150 (MQ=255)
                             gTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGc                                                                                                                          <  2:961786/150‑1 (MQ=255)
                               cATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCag                                                                                                                        <  2:625143/150‑1 (MQ=255)
                               cATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCag                                                                                                                        >  2:1383519/1‑150 (MQ=255)
                                    ccGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGtga                                                                                                                   <  2:129407/150‑1 (MQ=255)
                                      gTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGtgatg                                                                                                                 <  1:1260689/150‑1 (MQ=255)
                                      gTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGtgatg                                                                                                                 <  1:1383519/150‑1 (MQ=255)
                                         aTTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGaaa                                                                                                              >  2:1602453/1‑150 (MQ=255)
                                          ttttATCGGCCTTATGGATATAGTATTTCATGGCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCACACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATc                                                                                                            >  2:355219/1‑150 (MQ=255)
                                           tttATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATc                                                                                                            >  1:978802/1‑150 (MQ=255)
                                            ttATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGGCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCt                                                                                                           <  2:1364141/150‑1 (MQ=255)
                                               tCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCt                                                                                                        <  2:622118/150‑1 (MQ=255)
                                               tCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCt                                                                                                        <  2:1064865/150‑1 (MQ=255)
                                                 ggCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTcc                                                                                                      <  1:154680/150‑1 (MQ=255)
                                                 ggCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTcc                                                                                                      <  1:982303/150‑1 (MQ=255)
                                                  gCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTccc                                                                                                     >  2:110315/1‑150 (MQ=255)
                                                  gCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTccc                                                                                                     >  2:1076620/1‑150 (MQ=255)
                                                    cTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGc                                                                                                   <  2:409943/150‑1 (MQ=255)
                                                    cTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGc                                                                                                   <  2:473427/150‑1 (MQ=255)
                                                      tATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCAc                                                                                                 >  1:809459/1‑150 (MQ=255)
                                                       aTGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCAcg                                                                                                >  2:75176/1‑150 (MQ=255)
                                                          gATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACgcgc                                                                                             >  1:765873/1‑150 (MQ=255)
                                                             ataGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAg                                                                                          <  2:1233058/150‑1 (MQ=255)
                                                                  aTTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAg                                                                                     >  1:1543525/1‑150 (MQ=255)
                                                                  aTTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAg                                                                                     >  2:1461095/1‑150 (MQ=255)
                                                                   tttCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGc                                                                                    <  2:377278/150‑1 (MQ=255)
                                                                     tCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCaa                                                                                  <  1:687590/150‑1 (MQ=255)
                                                                       aTGGGCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATg                                                                                <  1:949002/150‑1 (MQ=255)
                                                                        tGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGt                                                                               <  1:359423/150‑1 (MQ=255)
                                                                           ccTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAAt                                                                            >  1:904675/1‑150 (MQ=255)
                                                                           ccTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAAt                                                                            >  2:138512/1‑150 (MQ=255)
                                                                               aaCTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTg                                                                        <  1:77886/150‑1 (MQ=255)
                                                                                 cTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCg                                                                      >  2:1521854/1‑150 (MQ=255)
                                                                                   ggCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCCCGCGCCAGTTGAGCAATGTAATCCTGCGGc                                                                    >  2:1318452/1‑150 (MQ=255)
                                                                                        ctcCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGt                                                               >  2:1465418/1‑150 (MQ=255)
                                                                                          ccGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCa                                                             >  2:681669/1‑150 (MQ=255)
                                                                                           cGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCAc                                                            >  2:1016563/1‑150 (MQ=255)
                                                                                              gTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGAAATGTAATCCTGCGGCCGCGTCACCAt                                                         >  2:587216/1‑150 (MQ=255)
                                                                                                tAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCa                                                       <  1:1461095/150‑1 (MQ=255)
                                                                                                   ggATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGAt                                                    <  1:337768/150‑1 (MQ=255)
                                                                                                    gATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATg                                                   >  2:198604/1‑150 (MQ=255)
                                                                                                     aTCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGa                                                  >  2:560078/1‑150 (MQ=255)
                                                                                                       cAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACa                                                <  1:1179749/150‑1 (MQ=255)
                                                                                                        aGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAg                                               <  1:1342184/150‑1 (MQ=255)
                                                                                                        aGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAg                                               >  2:322300/1‑150 (MQ=255)
                                                                                                           ccaaCTTTAATGTCATGACGGCTGATTTCATCAATCAGGCGGAACACCTGGCCGTTGATGGTTTCCAGATGAAGAGTGATGAAATCTGCTCCCGAAAGCGCCAGTTGAGAAATGTAATCCTGCAGCCGCGAAACCATCAGAAGAAAgaaa                                            >  2:710243/4‑147 (MQ=255)
                                                                                                               cTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCg                                        <  1:1283137/150‑1 (MQ=255)
                                                                                                                 ttCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGt                                      <  2:385202/150‑1 (MQ=255)
                                                                                                                  tcatgtcatgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGtt                                     <  1:1050502/150‑1 (MQ=255)
                                                                                                                   catgtcatgACGGCGGATTTCATCAATCAGGCGGAACGCCTGCCCGTAGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGttt                                    >  1:1141990/1‑150 (MQ=255)
                                                                                                                     tgtcatgACGGCGGATTTCATCAATCAGGCCGAACCCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGAAGAAATCTGCTCCCGCAAGCGCCAGTTCCGCAATGTAATCCTGCGGCCGCGACACCATCAGATGACAGTCGAGCGCTTAAg                                  >  2:947306/1‑150 (MQ=255)
                                                                                                                      gtcatgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGt                                 <  2:1418356/150‑1 (MQ=255)
                                                                                                                      gtcatgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGt                                 >  1:344460/1‑150 (MQ=255)
                                                                                                                      gtcatgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGt                                 >  1:1342300/1‑150 (MQ=255)
                                                                                                                          tgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGcc                             >  1:1221202/1‑150 (MQ=255)
                                                                                                                                cggATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGtttt                       <  2:1558118/150‑1 (MQ=255)
                                                                                                                                 ggATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATTTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGttttt                      <  1:1016563/150‑1 (MQ=255)
                                                                                                                                 ggATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGttttt                      <  1:322300/150‑1 (MQ=255)
                                                                                                                                       catcaATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCt                >  2:1054008/1‑150 (MQ=255)
                                                                                                                                         tcaatcaGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGa              >  1:1561071/1‑150 (MQ=255)
                                                                                                                                           aatcaGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACt            >  2:255403/1‑150 (MQ=255)
                                                                                                                                            atcaGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACtt           >  1:29022/1‑150 (MQ=255)
                                                                                                                                            atcaGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACtt           >  1:1198443/1‑150 (MQ=255)
                                                                                                                                                ggCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGGTTAGTTGCCAGTTTTTTAACCTGACTTACga       >  1:711181/1‑150 (MQ=255)
                                                                                                                                                 gCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACgaa      >  1:918168/1‑150 (MQ=255)
                                                                                                                                                     aaCGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAAc  >  1:974501/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAAC  >  NC_000913/4308043‑4308341

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