Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,308,193 A→C 6.7% F96V (TTC→GTC) ‡ 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,1930AC 6.7% 257.6 / 8.8 87F96V (TTC→GTC) ‡alsED‑allulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base C (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.09e-01

GTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACG  >  NC_000913/4308046‑4308342
                                                                                                                                                   |                                                                                                                                                     
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:1010412/1‑150 (MQ=255)
    tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGTTCCGCCTg                                                                                                                                                 >  1:1004275/1‑150 (MQ=255)
    tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTg                                                                                                                                                 <  2:802678/150‑1 (MQ=255)
    tCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTg                                                                                                                                                 <  2:1147934/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                                  >  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: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)
                                                                                                                                                   aCGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACg  <  1:204064/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACG  >  NC_000913/4308046‑4308342

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: