Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,641 C→T 15.9% Q118Q (CAG→CAA tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,6410CT15.9% 1148.1 / 132.3 422Q118Q (CAG→CAAtktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (181/172);  new base T (34/35);  total (215/207)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 431 total aligned reads displayed.

GAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCAT  >  CP009273/3074492‑3074790
                                                                                                                                                     |                                                                                                                                                     
    aaaCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAc                                                                                                                                                   <  2:189936/150‑1 (MQ=255)
     aaCGTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCAAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAcc                                                                                                                                                  <  1:93447/150‑1 (MQ=255)
      gcTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAcct                                                                                                                                                 <  1:296482/149‑2 (MQ=255)
      gcTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAcc                                                                                                                                                  >  2:692078/2‑149 (MQ=255)
      aCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAcca                                                                                                                                                 <  2:287966/150‑2 (MQ=255)
          cGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAg                                                                                                                                               <  1:429527/148‑1 (MQ=255)
          cGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAg                                                                                                                                               <  2:371767/148‑1 (MQ=255)
           gTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCg                                                                                                                                             >  1:231363/1‑149 (MQ=255)
            tGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                                                                           >  2:867344/1‑150 (MQ=255)
            tGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                                                                           <  2:1024581/150‑1 (MQ=255)
             gcgagaTGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACggtc                                                                                                                                          >  2:565077/3‑148 (MQ=255)
              ggAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGAcc                                                                                                                                         <  1:739577/150‑1 (MQ=255)
               gagaTGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCg                                                                                                                                        <  1:835353/150‑1 (MQ=255)
                   tGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtgg                                                                                                                                    <  1:720603/150‑1 (MQ=255)
                     ccTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCggtcgtc                                                                                                                                  <  2:981380/150‑2 (MQ=255)
                      cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtgg                                                                                                                                 >  2:421277/1‑150 (MQ=255)
                      cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtgg                                                                                                                                 <  1:652650/150‑1 (MQ=255)
                          cattaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCCTGACCTAACAGGCCGGTCGTCGGTTc                                                                                                                             >  2:514843/7‑150 (MQ=255)
                           atcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTcc                                                                                                                            <  2:472747/150‑1 (MQ=255)
                           atcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTcc                                                                                                                            <  1:829541/150‑1 (MQ=255)
                              aagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgac                                                                                                                          <  1:320780/147‑1 (MQ=255)
                              aagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgac                                                                                                                          <  1:110087/147‑1 (MQ=255)
                               agCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGCGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgacgc                                                                                                                        >  1:569318/2‑148 (MQ=255)
                                 cAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGccc                                                                                                                      >  1:800489/1‑149 (MQ=255)
                                  aGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAg                                                                                                                     <  2:161321/150‑1 (MQ=255)
                                   gCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGc                                                                                                                    >  1:984546/1‑150 (MQ=255)
                                   gCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGc                                                                                                                    <  1:1010835/150‑1 (MQ=255)
                                   gCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGc                                                                                                                    >  1:172937/1‑150 (MQ=255)
                                     cGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCgg                                                                                                                  >  2:813124/1‑150 (MQ=255)
                                      gTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgt                                                                                                                 <  2:1036711/150‑1 (MQ=255)
                                      gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                                  >  1:890039/1‑144 (MQ=255)
                                      gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                                  >  2:318755/1‑144 (MQ=255)
                                      gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                                  >  2:426442/1‑144 (MQ=255)
                                      gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGGGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                                  >  1:129201/1‑144 (MQ=255)
                                         gCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcgca                                                                                                                >  2:828490/1‑141 (MQ=255)
                                         gCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcgca                                                                                                                >  2:592764/1‑141 (MQ=255)
                                          cccATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTaa                                                                                                             <  1:847448/150‑1 (MQ=255)
                                            cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                                           >  1:514677/1‑150 (MQ=255)
                                            cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                                           >  1:68217/1‑150 (MQ=255)
                                              tGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCa                                                                                                         <  1:608547/150‑1 (MQ=255)
                                               gAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGACCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCAc                                                                                                        >  2:130137/1‑150 (MQ=255)
                                                 aGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACtt                                                                                                      <  2:667180/150‑1 (MQ=255)
                                                 aGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACtt                                                                                                      <  1:48056/150‑1 (MQ=255)
                                                    cGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCg                                                                                                   >  1:1029073/1‑150 (MQ=255)
                                                       aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                                                >  2:304143/1‑150 (MQ=255)
                                                           gtAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                                                                            >  1:330466/1‑150 (MQ=255)
                                                               tGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGa                                                                                        <  2:38298/150‑1 (MQ=255)
                                                                    gACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGtttt                                                                                   <  1:999091/150‑1 (MQ=255)
                                                                       aaTGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATGCCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTaga                                                                                <  2:852168/150‑1 (MQ=255)
                                                                         tGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGt                                                                              <  2:651087/150‑1 (MQ=255)
                                                                             gTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCa                                                                           <  1:347845/149‑1 (MQ=255)
                                                                             gTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCa                                                                           <  2:231363/149‑1 (MQ=255)
                                                                             gTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAg                                                                          <  1:140684/150‑1 (MQ=255)
                                                                                  ccGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGAc                                                                     <  1:492073/150‑1 (MQ=255)
                                                                                     gACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGa                                                                  <  1:53644/150‑1 (MQ=255)
                                                                                       cGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAg                                                                <  2:963290/150‑1 (MQ=255)
                                                                                        ggTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt                                                               <  2:450637/150‑1 (MQ=255)
                                                                                         gTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtt                                                              >  1:273417/1‑150 (MQ=255)
                                                                                           tAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtttt                                                            >  1:634750/1‑150 (MQ=255)
                                                                                            aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                                           >  2:916751/1‑150 (MQ=255)
                                                                                            aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                                           >  1:416968/1‑150 (MQ=255)
                                                                                            aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                                           >  2:550234/1‑150 (MQ=255)
                                                                                            aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtttt                                                            <  1:1000151/149‑1 (MQ=255)
                                                                                             aaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                                           <  2:750295/149‑1 (MQ=255)
                                                                                             aaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                                           >  1:897655/1‑149 (MQ=255)
                                                                                                 gcgcCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAAAGTGCAGCGGACGGAAGTTTTTCAGtt                                                      >  2:196852/1‑150 (MQ=255)
                                                                                                     cgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttcttc                                                  >  2:712857/1‑150 (MQ=255)
                                                                                                       ccAGCGTTTTTTCTGCAATCGTCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                                                <  1:397356/150‑1 (MQ=255)
                                                                                                       ccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                                                <  2:270108/150‑1 (MQ=255)
                                                                                                       ccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                                                <  2:736922/150‑1 (MQ=255)
                                                                                                        cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                                               <  2:837178/150‑1 (MQ=255)
                                                                                                          gCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCg                                             <  1:558568/150‑1 (MQ=255)
                                                                                                          gCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCg                                             <  1:151830/150‑1 (MQ=255)
                                                                                                            gTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                            >  1:310573/1‑149 (MQ=255)
                                                                                                            gTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                            >  2:607017/1‑149 (MQ=255)
                                                                                                             ttttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                            >  1:183412/1‑148 (MQ=255)
                                                                                                             ttttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                            >  1:771251/1‑148 (MQ=255)
                                                                                                              tttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                            >  2:53236/1‑147 (MQ=255)
                                                                                                               ttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAg                                         <  1:376725/149‑1 (MQ=255)
                                                                                                                  tCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCg                                     >  1:957954/1‑150 (MQ=255)
                                                                                                                      cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGt                                    <  2:91054/147‑1 (MQ=255)
                                                                                                                      cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTa                                   <  1:874487/148‑1 (MQ=255)
                                                                                                                      cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAAc                                 >  1:755187/1‑150 (MQ=255)
                                                                                                                       aaTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAAcc                                >  1:451539/1‑150 (MQ=255)
                                                                                                                       aaTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAAcc                                <  2:887934/150‑1 (MQ=255)
                                                                                                                              aTACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAgg                         <  2:351390/150‑1 (MQ=255)
                                                                                                                                aCCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTg                       >  2:583754/1‑150 (MQ=255)
                                                                                                                                  cGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTgca                     <  1:453810/150‑1 (MQ=255)
                                                                                                                                   gACTGCGTTGGCAATACCCTGACCCAGCGGCCCGGTGGTGGTTTCCACACCAGCGGTGTAATACACCTCCGGGAGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTgcag                    >  2:933842/1‑150 (MQ=255)
                                                                                                                                    aCTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTgcag                    >  2:147671/1‑149 (MQ=255)
                                                                                                                                        cGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGc               >  2:173515/1‑150 (MQ=255)
                                                                                                                                          ttGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTg             <  2:630263/150‑1 (MQ=255)
                                                                                                                                            ggCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTg             >  2:279664/1‑148 (MQ=255)
                                                                                                                                              cAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTg             >  2:330977/1‑146 (MQ=255)
                                                                                                                                               aaTACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGAt        <  1:198199/150‑1 (MQ=255)
                                                                                                                                                aTACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATc       <  2:227886/150‑1 (MQ=255)
                                                                                                                                                 tACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGAt        >  2:39881/1‑148 (MQ=255)
                                                                                                                                                  aCCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAg     >  1:731384/1‑150 (MQ=255)
                                                                                                                                                   cccTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGc    <  2:504113/150‑1 (MQ=255)
                                                                                                                                                    ccTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCa   <  2:772343/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCAT  >  CP009273/3074492‑3074790

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