Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,667,358 G→A 15.9% W37* (TGG→TAG)  U712_08585 → Flagellar biosynthetic protein fliZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,667,3580GA15.9% 1020.0 / 118.7 344W37* (TGG→TAG) U712_08585Flagellar biosynthetic protein fliZ
Reads supporting (aligned to +/- strand):  ref base G (132/157);  new base A (26/29);  total (158/186)
Fisher's exact test for biased strand distribution p-value = 8.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 350 total aligned reads displayed.

AAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACTGTTTGTTA  >  CP006881/1667209‑1667507
                                                                                                                                                     |                                                                                                                                                     
   cAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGtt                                                                                                                                                    <  1:1751277/150‑1 (MQ=255)
   cAAACAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGtt                                                                                                                                                    <  1:100812/150‑1 (MQ=255)
         acaaGACATTTAAGTAAAGGCAGTACATCTGTGGAAAAAGTGATTATATTTTTTTCGTTTTATTTGTTTTTTCGTTTAATAAACAGTCCAACCATTTGCTGATGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCaaaa                                                                                                                                              <  2:1629363/150‑1 (MQ=255)
          caaAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAGTATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAaaag                                                                                                                                             <  2:780281/150‑1 (MQ=255)
          caaAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAaaag                                                                                                                                             >  1:1160361/1‑150 (MQ=255)
          caaAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAaaag                                                                                                                                             <  1:3473553/150‑1 (MQ=255)
           aaaaCATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAaaaga                                                                                                                                            <  1:2918000/150‑1 (MQ=255)
           aaaaCATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAaaaga                                                                                                                                            >  2:1973668/1‑150 (MQ=255)
            aaaCATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAaaagaa                                                                                                                                           >  2:1204150/1‑150 (MQ=255)
            aaaCATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAaaagaa                                                                                                                                           <  1:2028380/150‑1 (MQ=255)
               cATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAaagaaaga                                                                                                                                        <  1:4018151/150‑1 (MQ=255)
               cATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAaagaaaga                                                                                                                                        >  2:5451965/1‑150 (MQ=255)
                 ttaaattaaaGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATg                                                                                                                                      <  1:4121743/150‑1 (MQ=255)
                 ttaaattaaaGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATg                                                                                                                                      >  2:1777152/1‑150 (MQ=255)
                 ttaaattaaaGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATg                                                                                                                                      <  1:3802699/150‑1 (MQ=255)
                     attaaaGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGaaaa                                                                                                                                  >  2:2653928/1‑150 (MQ=255)
                     attaaaGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGaaaa                                                                                                                                  <  1:226307/150‑1 (MQ=255)
                      ttaaaGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGaaaaa                                                                                                                                 >  1:835550/1‑150 (MQ=255)
                         aaGGGTGTCCGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACt                                                                                                                              <  1:5576861/150‑1 (MQ=255)
                         aaGGGTGTACGACTGTTGAAAAAGAGCCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACt                                                                                                                              <  1:4561918/150‑1 (MQ=255)
                             gtgtACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTTCAACCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAg                                                                                                                          <  1:883130/150‑1 (MQ=255)
                             gtgtACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAg                                                                                                                          >  2:5416862/1‑150 (MQ=255)
                                tACGACTGTTGAAAAAGAGTCAATATTTTATTCTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTCCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGatc                                                                                                                       <  1:596524/150‑1 (MQ=255)
                                  cGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGatcaa                                                                                                                     >  2:1027518/1‑150 (MQ=255)
                                    aCTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGatcaatc                                                                                                                   >  2:2240394/1‑150 (MQ=255)
                                     cTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGatcaatca                                                                                                                  <  2:4847901/150‑1 (MQ=255)
                                      tGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCag                                                                                                                 >  2:3305062/1‑150 (MQ=255)
                                       gTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCaga                                                                                                                >  2:2484896/1‑150 (MQ=255)
                                        ttGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCagag                                                                                                               >  2:3616565/1‑150 (MQ=255)
                                         tGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGc                                                                                                              <  1:206559/150‑1 (MQ=255)
                                          gAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCa                                                                                                             <  2:4615349/150‑1 (MQ=255)
                                              aaGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAaaag                                                                                                         >  2:5167505/1‑150 (MQ=255)
                                              aaGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAaaag                                                                                                         <  1:429965/150‑1 (MQ=255)
                                               agagTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATTAATCAGAGCAaaaga                                                                                                        >  2:5692641/1‑150 (MQ=255)
                                                 agTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAaaagaaa                                                                                                      >  1:2419274/1‑150 (MQ=255)
                                                 agTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAaaagaaa                                                                                                      <  1:2770131/150‑1 (MQ=255)
                                                    ctatatTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAagaaagaa                                                                                                   <  2:5779455/148‑1 (MQ=255)
                                                    cAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAagaaagaa                                                                                                   >  2:4892806/1‑150 (MQ=255)
                                                      atatTTTATTGTTTTTATTTGTTTTTTGATTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGATAAAACTGCAGATCAATCAGAGCAAAAGAAAGaaaa                                                                                                 <  2:4528925/150‑1 (MQ=255)
                                                      atatTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGaaaa                                                                                                 <  2:142285/150‑1 (MQ=255)
                                                      atatTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGaaaa                                                                                                 >  1:1436577/1‑150 (MQ=255)
                                                         ttttATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAaac                                                                                              >  1:5243421/1‑150 (MQ=255)
                                                         ttttATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAaac                                                                                              >  2:3090913/1‑150 (MQ=255)
                                                             tttGTTTTTTTTTGTTTTTTCGTTTTATCCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATACAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAacaaca                                                                                          <  2:4112233/149‑1 (MQ=255)
                                                             aTTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAacaaca                                                                                          >  2:3456867/1‑150 (MQ=255)
                                                                  ttttATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAAc                                                                                     >  2:4978332/1‑150 (MQ=255)
                                                                    ttATTTGTTTTTTCGTTTTATTCAGTTTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAActg                                                                                   <  1:1057637/150‑1 (MQ=255)
                                                                     tATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAActgc                                                                                  <  1:1952290/150‑1 (MQ=255)
                                                                     tATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAActgc                                                                                  <  1:3954086/150‑1 (MQ=255)
                                                                      aTTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAActgct                                                                                 >  1:2984549/1‑150 (MQ=255)
                                                                      aTTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAActgct                                                                                 >  2:524942/1‑150 (MQ=255)
                                                                       tttGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAActgctg                                                                                <  1:4366186/150‑1 (MQ=255)
                                                                        ttGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCtga                                                                               >  2:5617593/1‑150 (MQ=255)
                                                                         tGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCtgat                                                                              <  2:4334237/150‑1 (MQ=255)
                                                                          gTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCtgatg                                                                             >  2:2151344/1‑150 (MQ=255)
                                                                          gTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCtgatg                                                                             >  2:4700408/1‑150 (MQ=255)
                                                                          gTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCtgatg                                                                             >  2:3134841/1‑150 (MQ=255)
                                                                           ttttttCGCTTTATTCATTGTTCATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGATAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCtgatga                                                                            <  1:5115426/150‑1 (MQ=255)
                                                                            tttttCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATgag                                                                           <  1:3879395/150‑1 (MQ=255)
                                                                             ttttCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAAATGAAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATGAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATgaga                                                                          >  1:5628605/1‑150 (MQ=255)
                                                                              tttCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGAc                                                                         <  1:2706144/150‑1 (MQ=255)
                                                                              tttCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGAc                                                                         >  2:4052630/1‑150 (MQ=255)
                                                                               ttCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGAc                                                                         >  2:916394/1‑149 (MQ=255)
                                                                                tCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACgg                                                                       >  1:2558897/1‑150 (MQ=255)
                                                                                    tttATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACggaggg                                                                   <  2:4264379/150‑1 (MQ=255)
                                                                                    tttATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACggaggg                                                                   <  1:2484896/150‑1 (MQ=255)
                                                                                     ttATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACggaggga                                                                  >  2:3863734/1‑150 (MQ=255)
                                                                                       aTTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGAGAGc                                                                >  2:1110901/1‑150 (MQ=255)
                                                                                          cAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGgc                                                             >  2:5111578/1‑150 (MQ=255)
                                                                                                aCATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAGAACTGGTGATGAAACGGAGCGAGCGGGTGATcc                                                       >  2:2065479/1‑150 (MQ=255)
                                                                                                aCATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTcc                                                       <  1:4457806/150‑1 (MQ=255)
                                                                                                   tCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTc                                                    <  1:677613/150‑1 (MQ=255)
                                                                                                   tCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTc                                                    <  1:3134841/150‑1 (MQ=255)
                                                                                                    ccGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCt                                                   >  1:3004050/1‑150 (MQ=255)
                                                                                                    ccGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCt                                                   <  1:534279/150‑1 (MQ=255)
                                                                                                      gATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGt                                                 >  1:2114903/1‑150 (MQ=255)
                                                                                                      gATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGt                                                 >  1:137088/1‑150 (MQ=255)
                                                                                                       aTTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGGTGCTCCTTCTGtc                                                >  2:678796/1‑150 (MQ=255)
                                                                                                          gctgctgcCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCa                                             <  1:3305062/150‑1 (MQ=255)
                                                                                                          gctgctgcCGCGGCAGACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCa                                             <  1:4892806/150‑1 (MQ=255)
                                                                                                                gcCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCtttt                                       <  2:2419274/150‑1 (MQ=255)
                                                                                                                 ccGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTg                                      >  2:2888945/1‑150 (MQ=255)
                                                                                                                    cgGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGAtt                                   >  1:1645609/1‑150 (MQ=255)
                                                                                                                     ggCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGAttt                                  <  1:2213308/150‑1 (MQ=255)
                                                                                                                     ggCAGACTCTGATAATTCAACTGTAAACGAATGGTGTCAAAAGAAAGATGCAAAAACTGCAGATGAATGAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGAttt                                  >  1:3651025/1‑150 (MQ=255)
                                                                                                                         gACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTa                              >  2:3770675/1‑150 (MQ=255)
                                                                                                                         gACTCTGATAATTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTa                              >  1:1307692/1‑150 (MQ=255)
                                                                                                                                aTAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAgatga                       <  1:4978332/150‑1 (MQ=255)
                                                                                                                                  aaTTCAACTGTAAACGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGAtt                     >  2:5543298/1‑150 (MQ=255)
                                                                                                                                     tCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTc                  <  1:4745620/150‑1 (MQ=255)
                                                                                                                                      cAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCg                 >  1:4261157/1‑150 (MQ=255)
                                                                                                                                            tAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTAc           <  2:356785/150‑1 (MQ=255)
                                                                                                                                              aaCGAATAGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtg         <  1:190746/150‑1 (MQ=255)
                                                                                                                                                cGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgtt       <  1:3090913/150‑1 (MQ=255)
                                                                                                                                                cGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgtt       <  2:4417138/150‑1 (MQ=255)
                                                                                                                                                 gAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgttt      <  1:3616565/150‑1 (MQ=255)
                                                                                                                                                  aaTGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgtttg     <  1:1815302/150‑1 (MQ=255)
                                                                                                                                                  aaTGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgtttg     >  2:1525447/1‑150 (MQ=255)
                                                                                                                                                   aTGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgtttgt    >  2:719701/1‑150 (MQ=255)
                                                                                                                                                    tGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACtgtttgtt   <  1:788395/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAATCAGAGCAAAAGAAAGAAAAAACAACAAAAACTGCTGATGAGACGGAGGGAGCGGCTGCTCCTTCTGTCTCAGCTTTTGATTTTGTAAAGATGATTTTCGCTTTACTGTTTGTTA  >  CP006881/1667209‑1667507

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