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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927381,5620CA66.7% 249.4 / 260.2 225I236I (ATC→ATAleuOglobal transcription factor
Reads supporting (aligned to +/- strand):  ref base C (42/33);  new base A (88/62);  total (130/95)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

TGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAACAGCAG  >  CP009273/81414‑81708
                                                                                                                                                    |                                                                                                                                                  
tGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATcg                                                                                                                                                   >  1:4736628/1‑150 (MQ=255)
 gCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATagc                                                                                                                                                  <  2:3837085/150‑1 (MQ=255)
    ggTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGtc                                                                                                                                               >  2:3358567/1‑149 (MQ=255)
       aGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCa                                                                                                                                            <  2:1006142/150‑1 (MQ=255)
           aGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGc                                                                                                                                        >  2:4026768/1‑150 (MQ=255)
              aaaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATg                                                                                                                                     >  1:679337/1‑150 (MQ=255)
               aaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATgg                                                                                                                                    >  1:1027849/1‑150 (MQ=255)
               aaaaTCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATgg                                                                                                                                    <  1:3656000/150‑1 (MQ=255)
                   tcatcCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCCGTATAGCGTATCAGGGCATGGCAat                                                                                                                                >  1:3392708/1‑150 (MQ=255)
                       ccAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAatgatg                                                                                                                            >  2:1102834/1‑150 (MQ=255)
                         aacaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAg                                                                                                                          <  2:4419666/150‑1 (MQ=255)
                          acaaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGc                                                                                                                         <  1:3833186/150‑1 (MQ=255)
                           caaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTAATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCg                                                                                                                        <  1:4758077/150‑1 (MQ=255)
                            aaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGt                                                                                                                       >  2:4526788/1‑150 (MQ=255)
                            aaTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGt                                                                                                                       >  1:3820809/1‑150 (MQ=255)
                             aTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTATATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTa                                                                                                                      >  1:474953/1‑150 (MQ=255)
                             aTTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTa                                                                                                                      <  1:3559704/150‑1 (MQ=255)
                              ttcaGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTAc                                                                                                                     <  1:955514/147‑1 (MQ=255)
                              ttAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTAc                                                                                                                     <  1:955672/150‑1 (MQ=255)
                                aaGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACtt                                                                                                                   >  1:1488006/1‑150 (MQ=255)
                                     cccGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCgt                                                                                                              >  1:4864968/1‑150 (MQ=255)
                                        gTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCgtggt                                                                                                           >  2:3548089/1‑150 (MQ=255)
                                        gTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCgtggt                                                                                                           <  1:2233464/150‑1 (MQ=255)
                                           aCTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTc                                                                                                        <  2:855434/150‑1 (MQ=255)
                                            cTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGATCGTACTTAGCGTGGTGTCg                                                                                                       >  2:4117555/1‑150 (MQ=255)
                                              gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCa                                                                                                     >  2:2691482/1‑150 (MQ=255)
                                              gAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCa                                                                                                     >  2:3267722/1‑150 (MQ=255)
                                                 aCATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAAc                                                                                                  <  1:1249185/150‑1 (MQ=255)
                                                  cATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTATCGCGGCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCCGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACg                                                                                                 >  1:3360512/1‑150 (MQ=255)
                                                     gatgTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCACTGTTGAGCGTACTTAGCGTGGTGTCGCAATCGCAt                                                                                              >  2:3040511/1‑150 (MQ=255)
                                                     gatgTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGCCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCAt                                                                                              >  2:2068136/1‑150 (MQ=255)
                                                       tgTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCAttt                                                                                            <  2:1805365/150‑1 (MQ=255)
                                                         tttATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTgg                                                                                          >  2:507014/1‑150 (MQ=255)
                                                         tttATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTgg                                                                                          >  2:2726043/1‑150 (MQ=255)
                                                              aaCGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCt                                                                                     >  1:1849893/1‑150 (MQ=255)
                                                                cGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTAt                                                                                   >  2:1183598/1‑150 (MQ=255)
                                                                 gAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTAtt                                                                                  <  1:4676382/150‑1 (MQ=255)
                                                                 gAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTAtt                                                                                  <  1:2941274/150‑1 (MQ=255)
                                                                  aacaacATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAATCAAGCCACTATCGCGTATCAGGGCATGGCAATTATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTg                                                                                 >  1:420739/1‑150 (MQ=255)
                                                                  aacaacATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTg                                                                                 >  2:255410/1‑150 (MQ=255)
                                                                      acaTGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGcgcc                                                                             <  1:4592118/150‑1 (MQ=255)
                                                                      acaTGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGcgcc                                                                             <  2:1551009/150‑1 (MQ=255)
                                                                      acaTGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGcgcc                                                                             <  1:54672/150‑1 (MQ=255)
                                                                         tGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATTACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCcgcg                                                                          <  2:2675547/150‑1 (MQ=255)
                                                                           cggcggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGtt                                                                        <  2:804119/150‑1 (MQ=255)
                                                                           cggcggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGtt                                                                        >  2:1255647/1‑150 (MQ=255)
                                                                            ggcggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtg                                                                       >  2:1711379/1‑150 (MQ=255)
                                                                            ggcggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtg                                                                       <  2:3139655/150‑1 (MQ=255)
                                                                              cggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggc                                                                     >  1:747078/1‑150 (MQ=255)
                                                                               ggTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggct                                                                    >  1:3275436/1‑150 (MQ=255)
                                                                                gTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTtggctg                                                                   <  1:2010214/150‑1 (MQ=255)
                                                                                   tCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTggctggct                                                                <  1:1888584/150‑1 (MQ=255)
                                                                                     gCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGa                                                              >  1:4015388/1‑150 (MQ=255)
                                                                                     gCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGa                                                              >  1:62168/1‑150 (MQ=255)
                                                                                      cTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGaa                                                             >  2:2192876/1‑150 (MQ=255)
                                                                                       tCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAag                                                            >  2:316849/1‑150 (MQ=255)
                                                                                            cGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAACGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTc                                                       >  2:847887/1‑150 (MQ=255)
                                                                                             gTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCg                                                      >  1:1705157/1‑150 (MQ=255)
                                                                                               ttCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCt                                                    <  2:435873/150‑1 (MQ=255)
                                                                                               ttCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCt                                                    <  2:4765338/150‑1 (MQ=255)
                                                                                               ttCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCt                                                    <  1:3076328/150‑1 (MQ=255)
                                                                                                 cgcgTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGa                                                  >  1:2341586/1‑150 (MQ=255)
                                                                                                  gcgTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGaa                                                 >  2:2503609/1‑150 (MQ=255)
                                                                                                   cgTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAAt                                                <  1:2279370/150‑1 (MQ=255)
                                                                                                      cATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCt                                             >  1:1141592/1‑150 (MQ=255)
                                                                                                       aTTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCtt                                            <  2:2514736/150‑1 (MQ=255)
                                                                                                         ttAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAg                                          <  2:2805101/150‑1 (MQ=255)
                                                                                                         ttAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAg                                          <  1:3733720/150‑1 (MQ=255)
                                                                                                            gTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAAt                                       <  2:1431620/150‑1 (MQ=255)
                                                                                                             tCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAAtt                                      <  1:4253381/150‑1 (MQ=255)
                                                                                                             tCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAAtt                                      <  1:3374040/150‑1 (MQ=255)
                                                                                                              cAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTa                                     >  2:1085183/1‑150 (MQ=255)
                                                                                                               aaCCTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTAc                                    >  2:318242/1‑150 (MQ=255)
                                                                                                                  cTTGGTATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAgg                                 <  2:3473345/150‑1 (MQ=255)
                                                                                                                      gTATGACACGGTAGATACGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTAtt                             <  1:3310982/150‑1 (MQ=255)
                                                                                                                       tATGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTa                            >  1:2453517/1‑150 (MQ=255)
                                                                                                                         tGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAcc                          >  2:4860347/1‑150 (MQ=255)
                                                                                                                         tGACACGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTAcc                          <  2:4703792/150‑1 (MQ=255)
                                                                                                                          gACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCg                         >  1:3837227/1‑150 (MQ=255)
                                                                                                                          gACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCg                         >  1:3837285/1‑150 (MQ=255)
                                                                                                                           acacGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGc                        >  1:3254237/1‑150 (MQ=255)
                                                                                                                            cacGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCt                       <  1:3622441/150‑1 (MQ=255)
                                                                                                                              cGGTAGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGc                     >  1:3422429/1‑150 (MQ=255)
                                                                                                                                  aGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGtt                 <  2:4864968/150‑1 (MQ=255)
                                                                                                                                  aGATAAGCAAGCCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGtt                 <  2:3681322/150‑1 (MQ=255)
                                                                                                                                    aTAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACAGCTGCCGTTaa               >  2:2032533/1‑150 (MQ=255)
                                                                                                                                       agcaagcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaac            <  2:4504470/150‑1 (MQ=255)
                                                                                                                                        gcaagcCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaaca           <  2:423948/150‑1 (MQ=255)
                                                                                                                                        gcaagcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaaca           <  1:1109730/150‑1 (MQ=255)
                                                                                                                                        gcaagcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaaca           >  2:1328859/1‑150 (MQ=255)
                                                                                                                                        gcaagcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaaca           <  1:3394367/150‑1 (MQ=255)
                                                                                                                                         caagcCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaacaa          <  1:140392/150‑1 (MQ=255)
                                                                                                                                          aagcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaacaaa         <  1:295552/150‑1 (MQ=255)
                                                                                                                                           agcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaacaaaa        <  1:4476538/150‑1 (MQ=255)
                                                                                                                                            gcCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaacaaaac       <  1:683600/150‑1 (MQ=255)
                                                                                                                                            gcCAGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTaaaacaaaac       >  1:87978/1‑150 (MQ=255)
                                                                                                                                               aGTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcagc    >  2:1920027/1‑150 (MQ=255)
                                                                                                                                                gTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcagca   >  2:3994426/1‑150 (MQ=255)
                                                                                                                                                gTATAGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcagca   <  1:3526820/150‑1 (MQ=255)
                                                                                                                                                 tataGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAAcagcag  >  1:3995113/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCTGGTAGCCAGCAAAAATCATCCAACAATTAAGGGCCCGTTACTGAAACATGATGTTTATAACGAACAACATGCGGCGGTTTCGCTCGATCGTTTCGCGTCATTTAGTCAACCTTGGTATGACACGGTAGATAAGCAAGCCAGTATCGCGTATCAGGGCATGGCAATGATGAGCGTACTTAGCGTGGTGTCGCAAACGCATTTGGTCGCTATTGCGCCGCGTTGGCTGGCTGAAGAGTTCGCTGAATCCTTAGAATTACAGGTATTACCGCTGCCGTTAAAACAAAACAGCAG  >  CP009273/81414‑81708

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: