Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,456,705 C→A 6.0% S72Y (TCC→TAC)  yfdF → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,456,7050CA 6.0% 829.4 / 20.4 216S72Y (TCC→TAC) yfdFuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (110/93);  new base A (6/7);  total (116/100)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 219 total aligned reads displayed.

AATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTC  >  CP009273/2456557‑2456854
                                                                                                                                                    |                                                                                                                                                     
 atgatgAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAAAAATTCCt                                                                                                                                                     >  2:1075507/1‑150 (MQ=255)
  tgatgaAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTc                                                                                                                                                    >  2:3583592/1‑150 (MQ=255)
    atgaAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAg                                                                                                                                                  <  1:1769732/150‑1 (MQ=255)
      gaAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGtt                                                                                                                                                <  2:2400663/150‑1 (MQ=255)
       aaGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTa                                                                                                                                               <  1:3784619/150‑1 (MQ=255)
       aaGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTa                                                                                                                                               >  1:485061/1‑150 (MQ=255)
       aaGTTAATGGATTAGTACAAGAGTTAAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTa                                                                                                                                               <  1:487783/150‑1 (MQ=255)
        aGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAATCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTACTCAGTTAc                                                                                                                                              >  2:1230131/1‑150 (MQ=255)
        aGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTAc                                                                                                                                              >  2:2040753/1‑150 (MQ=255)
          ttAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACtt                                                                                                                                            >  1:3352444/1‑150 (MQ=255)
           tAATGGATTAGTACAATAGTTCAAAAACCTTTTTAATGGTAAGGATGGAATAAGCACCTGTATTAAACATCTACGTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTACTATTACTAATTCCTCAGTTACTta                                                                                                                                           <  2:4106604/150‑1 (MQ=255)
           tAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTta                                                                                                                                           >  1:642567/1‑150 (MQ=255)
              tGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatat                                                                                                                                        >  2:2615326/1‑150 (MQ=255)
              tGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatat                                                                                                                                        <  1:2029328/150‑1 (MQ=255)
              tGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatat                                                                                                                                        >  2:1015086/1‑150 (MQ=255)
              tGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatat                                                                                                                                        <  2:2783106/150‑1 (MQ=255)
               ggATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGACGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAAATCCTCAGTAACTTatatt                                                                                                                                       >  2:856044/1‑150 (MQ=255)
                gATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatattg                                                                                                                                      >  2:4371915/1‑150 (MQ=255)
                         aaGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACt                                                                                                                             >  2:1937507/1‑150 (MQ=255)
                         aaGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACt                                                                                                                             <  2:3105042/150‑1 (MQ=255)
                         aaGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCATATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGAca                                                                                                                             >  2:3476679/1‑149 (MQ=255)
                          agagTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTc                                                                                                                            <  1:2311099/150‑1 (MQ=255)
                          agagTTCAAAAACCTTTTTAATGGTAAGGAAGGAAAAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTACTCAGTTACTTATATTGATATTGACTc                                                                                                                            <  1:3016724/150‑1 (MQ=255)
                                 aaaaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGAt                                                                                                                     <  1:1209399/150‑1 (MQ=255)
                                 aaaaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGAt                                                                                                                     >  2:3159721/1‑150 (MQ=255)
                                 aaaaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGAt                                                                                                                     >  2:272490/1‑150 (MQ=255)
                                  aaaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATa                                                                                                                    >  1:1807031/1‑150 (MQ=255)
                                   aaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATAc                                                                                                                   >  1:346245/1‑150 (MQ=255)
                                    aaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACa                                                                                                                  >  2:2776931/1‑150 (MQ=255)
                                     aCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAg                                                                                                                 >  2:4523334/1‑150 (MQ=255)
                                     aCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAg                                                                                                                 <  2:127938/150‑1 (MQ=255)
                                     aCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAAGTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAg                                                                                                                 <  2:127736/150‑1 (MQ=255)
                                       cTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGAc                                                                                                               <  1:1429557/150‑1 (MQ=255)
                                       cTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGAc                                                                                                               >  2:1470696/1‑150 (MQ=255)
                                            tAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCatat                                                                                                          >  2:4635480/1‑150 (MQ=255)
                                              aTGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTa                                                                                                        <  1:3706827/150‑1 (MQ=255)
                                                       aaggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTa                                                                                               >  1:4299627/1‑150 (MQ=255)
                                                        aggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTAt                                                                                              <  1:211036/150‑1 (MQ=255)
                                                        aggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTAt                                                                                              >  2:2880727/1‑150 (MQ=255)
                                                         ggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATc                                                                                             >  2:1402088/1‑150 (MQ=255)
                                                          gaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGAATTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCg                                                                                            >  1:612559/1‑150 (MQ=255)
                                                            aTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCAGACGAGTAAACGCCGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGAc                                                                                          <  2:4305031/150‑1 (MQ=255)
                                                                        aTTAAACATCTACTTGATCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCa                                                                              <  1:2478862/150‑1 (MQ=255)
                                                                        aTTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCa                                                                              >  2:2879603/1‑150 (MQ=255)
                                                                        aTTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCa                                                                              >  1:3077953/1‑150 (MQ=255)
                                                                         ttAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCaa                                                                             <  2:2369964/150‑1 (MQ=255)
                                                                          tAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAAt                                                                            <  2:3412447/150‑1 (MQ=255)
                                                                           aaaCATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAaccaataaccaata                                                                           >  1:1090679/1‑150 (MQ=255)
                                                                           aaaCATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATa                                                                           >  2:1587388/1‑150 (MQ=255)
                                                                            aaCATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAg                                                                          <  1:773398/150‑1 (MQ=255)
                                                                            aaCATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAg                                                                          >  2:3805433/1‑150 (MQ=255)
                                                                               aTCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAAt                                                                       >  1:1240367/1‑150 (MQ=255)
                                                                               aTCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAAt                                                                       >  1:1785418/1‑150 (MQ=255)
                                                                               aTCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAAt                                                                       >  2:1701193/1‑150 (MQ=255)
                                                                                tCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAAtt                                                                      >  1:2767524/1‑150 (MQ=255)
                                                                                 cTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTa                                                                     <  1:186401/150‑1 (MQ=255)
                                                                                 cTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTa                                                                     >  1:2470168/1‑150 (MQ=255)
                                                                                     ttGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTg                                                                 <  1:1060926/150‑1 (MQ=255)
                                                                                     ttGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTg                                                                 >  2:2383446/1‑150 (MQ=255)
                                                                                        aGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGa                                                              >  1:1187537/1‑150 (MQ=255)
                                                                                        aGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGa                                                              >  1:725782/1‑150 (MQ=255)
                                                                                        aGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGa                                                              <  1:2196854/150‑1 (MQ=255)
                                                                                          cTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAAc                                                            <  2:1526066/150‑1 (MQ=255)
                                                                                          cTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAAc                                                            >  2:1116236/1‑150 (MQ=255)
                                                                                          cTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAAc                                                            >  2:2211633/1‑150 (MQ=255)
                                                                                           ttATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACt                                                           <  2:4463314/150‑1 (MQ=255)
                                                                                            tataAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACta                                                          <  1:1571200/150‑1 (MQ=255)
                                                                                                aaaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATaaa                                                      <  2:2350787/150‑1 (MQ=255)
                                                                                                aaaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATaaa                                                      >  2:163358/1‑150 (MQ=255)
                                                                                                 aaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAg                                                     >  1:2794942/1‑150 (MQ=255)
                                                                                                 aaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAg                                                     >  1:1927394/1‑150 (MQ=255)
                                                                                                 aaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAg                                                     >  2:819310/1‑150 (MQ=255)
                                                                                                      gCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACaaa                                                <  1:4052425/150‑1 (MQ=255)
                                                                                                        cATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTACTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGa                                              <  2:3686850/150‑1 (MQ=255)
                                                                                                         aTACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGaa                                             <  1:4128895/150‑1 (MQ=255)
                                                                                                          tACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAAc                                            <  1:4349035/150‑1 (MQ=255)
                                                                                                                tAAACGAAGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTcc                                      <  2:1437770/150‑1 (MQ=255)
                                                                                                                   acgacgATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTa                                   <  2:3771626/150‑1 (MQ=255)
                                                                                                                     gacgaTCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACt                                 >  1:3764283/1‑150 (MQ=255)
                                                                                                                      acgaTCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTa                                <  2:1839732/150‑1 (MQ=255)
                                                                                                                       cgaTCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTAt                               >  1:3875383/1‑150 (MQ=255)
                                                                                                                          tCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCAt                            <  1:4468964/150‑1 (MQ=255)
                                                                                                                            cTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTa                          <  2:2013273/150‑1 (MQ=255)
                                                                                                                              tataGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAAt                        <  1:2152558/150‑1 (MQ=255)
                                                                                                                              tataGATTTAATATTAATAATTACTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAAt                        <  2:1953306/150‑1 (MQ=255)
                                                                                                                               ataGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAAGAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATg                       >  2:3613671/1‑150 (MQ=255)
                                                                                                                                 aGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGac                     <  2:1582029/150‑1 (MQ=255)
                                                                                                                                 aGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGac                     >  1:131619/1‑150 (MQ=255)
                                                                                                                                  gATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGaca                    >  1:365702/1‑150 (MQ=255)
                                                                                                                                  gATTTAATATTAATAATTACTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGaca                    >  2:879/1‑150 (MQ=255)
                                                                                                                                     ttaatattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGacaaca                 <  2:3281424/150‑1 (MQ=255)
                                                                                                                                        atattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAg              <  1:1118056/150‑1 (MQ=255)
                                                                                                                                            taataaTTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTga          >  2:1504401/1‑150 (MQ=255)
                                                                                                                                             aataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgag         <  2:2941325/150‑1 (MQ=255)
                                                                                                                                              ataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgaga        <  2:81313/150‑1 (MQ=255)
                                                                                                                                              ataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgaga        <  2:81465/150‑1 (MQ=255)
                                                                                                                                                aatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGagaag      <  1:2908539/150‑1 (MQ=255)
                                                                                                                                                aatTACTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGagaag      <  2:638175/150‑1 (MQ=255)
                                                                                                                                                 atTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTCTCATTAATGACAACATAGTTGagaaga     >  2:2662342/1‑150 (MQ=255)
                                                                                                                                                 atTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGagaaga     >  2:4026836/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTC  >  CP009273/2456557‑2456854

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: