Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,308,191 G→T 6.3% F96L (TTC→TTA) ‡ 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,1910GT 6.3% 294.7 / 5.1 88F96L (TTC→TTA) ‡alsED‑allulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base T (2/4);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

AACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAA  >  NC_000913/4308042‑4308337
                                                                                                                                                     |                                                                                                                                                  
aaCGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCgg                                                                                                                                                    <  1:741717/150‑1 (MQ=255)
aaCGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCgg                                                                                                                                                    >  2:154680/1‑150 (MQ=255)
 aCGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGa                                                                                                                                                   <  1:1144041/150‑1 (MQ=255)
    gTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACg                                                                                                                                                <  2:138312/150‑1 (MQ=255)
    gTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACg                                                                                                                                                <  2:586630/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:1010412/1‑150 (MQ=255)
        tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTg                                                                                                                                            >  1:1004275/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:473427/150‑1 (MQ=255)
                                                     cTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGc                                                                                               <  2:409943/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                             >  1:344460/1‑150 (MQ=255)
                                                                                                                       gtcatgACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGt                             <  2:1418356/150‑1 (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:1198443/1‑150 (MQ=255)
                                                                                                                                             atcaGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACtt       >  1:29022/1‑150 (MQ=255)
                                                                                                                                                 ggCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGGTTAGTTGCCAGTTTTTTAACCTGACTTACga   >  1:711181/1‑150 (MQ=255)
                                                                                                                                                  gCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACgaa  >  1:918168/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAA  >  NC_000913/4308042‑4308337

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: