Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 81,562 C→A 100% I236I (ATC→ATA leuO → global transcription factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927381,5620CA96.5% 847.5 / 19.5 258I236I (ATC→ATAleuOglobal transcription factor
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (129/120);  total (133/125)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 260 total aligned reads displayed.

CTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAACAGCAGA  >  CP009273/81416‑81709
                                                                                                                                                  |                                                                                                                                                   
    tAGCCAGCAAAAATCATCCAACAATTAAGGGCACTTTACTTAAACATGATGTTTATAACGAACAACATGCGGCGTTTTAGCTCGATCGTTTCTCGTCATTTAGTAAAAATTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATc                                                                                                                                              >  1:2234150/1‑150 (MQ=255)
     aGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCa                                                                                                                                             <  1:4578602/150‑1 (MQ=255)
     aGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCa                                                                                                                                             <  2:4021366/150‑1 (MQ=255)
      gCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAg                                                                                                                                            <  1:432490/150‑1 (MQ=255)
         aGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGc                                                                                                                                         >  1:89316/1‑150 (MQ=255)
           cAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCAt                                                                                                                                       >  1:4165588/1‑150 (MQ=255)
            aaaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATg                                                                                                                                      <  2:4423216/150‑1 (MQ=255)
             aaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATgg                                                                                                                                     <  2:3642783/150‑1 (MQ=255)
              aaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGc                                                                                                                                    <  1:4388210/150‑1 (MQ=255)
                atcatcCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCaa                                                                                                                                  <  1:4092121/150‑1 (MQ=255)
                atcatcCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCaa                                                                                                                                  <  1:1128212/150‑1 (MQ=255)
                  catcCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAatg                                                                                                                                <  1:3211806/150‑1 (MQ=255)
                   atcCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAatga                                                                                                                               <  2:4483053/150‑1 (MQ=255)
                   atcCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAatga                                                                                                                               >  2:84644/1‑150 (MQ=255)
                     ccAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAatgatg                                                                                                                             <  2:3296165/150‑1 (MQ=255)
                     ccAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAatgatg                                                                                                                             <  2:1491308/150‑1 (MQ=255)
                      caacaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAAtgatga                                                                                                                            <  2:4035268/150‑1 (MQ=255)
                      caacaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAAtgatga                                                                                                                            <  2:1701804/150‑1 (MQ=255)
                       aacaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAg                                                                                                                           >  1:3999903/1‑150 (MQ=255)
                       aacaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAg                                                                                                                           >  2:2180727/1‑150 (MQ=255)
                        acaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGc                                                                                                                          >  2:1205564/1‑150 (MQ=255)
                          aaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGTGCATGGCAATGATGAGCGt                                                                                                                        >  2:3500704/1‑150 (MQ=255)
                             tAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACt                                                                                                                     >  2:2407633/1‑150 (MQ=255)
                                 ggCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGc                                                                                                                 <  1:4462968/150‑1 (MQ=255)
                                  gCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCg                                                                                                                >  2:18969/1‑150 (MQ=255)
                                     cGTTACTGAAACATGCTGTTTATAACGAACAACCTGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCgtgg                                                                                                             <  2:470409/150‑1 (MQ=255)
                                     cGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTGGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCgtgg                                                                                                             <  2:3688768/150‑1 (MQ=255)
                                        tACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGgtgt                                                                                                          >  2:2942957/1‑150 (MQ=255)
                                         aCTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTc                                                                                                         <  2:4177360/150‑1 (MQ=255)
                                             aaccATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCaa                                                                                                     <  1:2509240/147‑1 (MQ=255)
                                              aaCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCaaa                                                                                                    >  2:2705387/1‑150 (MQ=255)
                                                  tgatgTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCa                                                                                                >  1:3782034/1‑150 (MQ=255)
                                                   gatgTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCAt                                                                                               >  1:3899964/1‑150 (MQ=255)
                                                     tgTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCAttt                                                                                             <  1:3770129/150‑1 (MQ=255)
                                                      gTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTg                                                                                            <  1:4337432/150‑1 (MQ=255)
                                                       tttATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTgg                                                                                           <  2:2226363/150‑1 (MQ=255)
                                                          ataACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCg                                                                                        <  1:2454105/150‑1 (MQ=255)
                                                            aaCGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCt                                                                                      <  2:3227717/150‑1 (MQ=255)
                                                             aCGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTa                                                                                     <  2:2387568/150‑1 (MQ=255)
                                                              cGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTAt                                                                                    >  2:3211928/1‑150 (MQ=255)
                                                               gAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTAtt                                                                                   >  1:650773/1‑150 (MQ=255)
                                                                aacaacATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTg                                                                                  >  1:4442679/1‑150 (MQ=255)
                                                                aacaacATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTg                                                                                  >  1:2111314/1‑150 (MQ=255)
                                                                aacaacATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTg                                                                                  >  1:1901195/1‑150 (MQ=255)
                                                                aacaacATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTg                                                                                  >  2:1328750/1‑150 (MQ=255)
                                                                      aTGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGcgccgc                                                                            >  1:2803801/1‑150 (MQ=255)
                                                                       tGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCcgcg                                                                           <  2:4083102/150‑1 (MQ=255)
                                                                             ggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggct                                                                     <  2:3943280/150‑1 (MQ=255)
                                                                             ggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggct                                                                     >  1:2978474/1‑150 (MQ=255)
                                                                               tttCGCTCGATCTTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggctgg                                                                   <  2:4118692/150‑1 (MQ=255)
                                                                                ttCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggctggc                                                                  >  1:960626/1‑150 (MQ=255)
                                                                                   gCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGa                                                               >  2:3021340/1‑150 (MQ=255)
                                                                                    cTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGaa                                                              >  2:1357403/1‑150 (MQ=255)
                                                                                         tCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGtt                                                         <  1:3632269/150‑1 (MQ=255)
                                                                                          cGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTc                                                        >  2:577603/1‑150 (MQ=255)
                                                                                          cGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGATGGCTGAAGAGTTc                                                        <  1:775603/150‑1 (MQ=255)
                                                                                          cGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGATATTGCGCCGCGTTGGCTGGCTGAAGAGTTc                                                        >  1:1034889/1‑150 (MQ=255)
                                                                                           gTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCg                                                       >  2:4605172/1‑150 (MQ=255)
                                                                                             ttCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCt                                                     >  1:372416/1‑150 (MQ=255)
                                                                                             ttCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCt                                                     >  2:1975589/1‑150 (MQ=255)
                                                                                                gcgTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGaa                                                  >  2:581628/1‑150 (MQ=255)
                                                                                                 cgTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAAt                                                 >  1:1421367/1‑150 (MQ=255)
                                                                                                  gTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATc                                                >  1:1960787/1‑150 (MQ=255)
                                                                                                  gTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATc                                                <  1:2875568/150‑1 (MQ=255)
                                                                                                  gTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATc                                                <  1:557200/150‑1 (MQ=255)
                                                                                                   tCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATcc                                               >  2:2273466/1‑150 (MQ=255)
                                                                                                    cATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCt                                              >  2:3793519/1‑150 (MQ=255)
                                                                                                    cATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTAGCTGAATCCt                                              >  1:1323763/1‑150 (MQ=255)
                                                                                                      tttAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTa                                            <  2:958536/150‑1 (MQ=255)
                                                                                                       ttAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAg                                           <  1:2668391/150‑1 (MQ=255)
                                                                                                         aGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGaa                                         >  2:2809308/1‑150 (MQ=255)
                                                                                                         aGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGaa                                         >  1:1773477/1‑150 (MQ=255)
                                                                                                          gTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAAt                                        >  1:3322246/1‑150 (MQ=255)
                                                                                                           tCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAAtt                                       <  2:274788/150‑1 (MQ=255)
                                                                                                               ccTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAg                                   >  1:1291724/1‑150 (MQ=255)
                                                                                                                 ttGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGt                                 <  2:1818157/150‑1 (MQ=255)
                                                                                                                  tGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGTAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTa                                <  2:286849/150‑1 (MQ=255)
                                                                                                                     tATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTa                             <  2:306172/150‑1 (MQ=255)
                                                                                                                      aTGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAc                            <  1:108066/150‑1 (MQ=255)
                                                                                                                      aTGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAc                            <  1:3591802/150‑1 (MQ=255)
                                                                                                                       tGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAcc                           <  1:62399/150‑1 (MQ=255)
                                                                                                                       tGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAcc                           <  1:454481/150‑1 (MQ=255)
                                                                                                                       tGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAcc                           <  1:3269463/150‑1 (MQ=255)
                                                                                                                          cacGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCt                        <  2:1845941/150‑1 (MQ=255)
                                                                                                                          cacGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCt                        <  2:67319/150‑1 (MQ=255)
                                                                                                                           acGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTg                       <  2:3563229/150‑1 (MQ=255)
                                                                                                                             ggTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGcc                     <  2:3650279/150‑1 (MQ=255)
                                                                                                                               tAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGt                   >  1:4009353/1‑150 (MQ=255)
                                                                                                                                aGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGtt                  <  1:1199598/150‑1 (MQ=255)
                                                                                                                                 gATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTa                 >  2:697274/1‑150 (MQ=255)
                                                                                                                                  aTAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaa                <  2:3709437/150‑1 (MQ=255)
                                                                                                                                  aTAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaa                >  1:4207067/1‑150 (MQ=255)
                                                                                                                                  aTAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaa                <  2:3709446/150‑1 (MQ=255)
                                                                                                                                   tAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaa               <  1:18969/150‑1 (MQ=255)
                                                                                                                                      gcaagcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaaca            <  2:1182539/150‑1 (MQ=255)
                                                                                                                                           ccAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAaaacaaaaca       >  1:1110675/1‑150 (MQ=255)
                                                                                                                                            cAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcag      >  1:3102370/1‑150 (MQ=255)
                                                                                                                                              gTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcagca    <  1:1615251/150‑1 (MQ=255)
                                                                                                                                              gTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCAGTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcagca    >  2:868463/1‑150 (MQ=255)
                                                                                                                                                aTCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAACAGCAGa  >  1:3008519/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAACAGCAGA  >  CP009273/81416‑81709

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