Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,470,017 G→A 93.6% M157I (ATG→ATA yjgN → DUF898 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,470,0170GA93.6% 652.7 / 22.8 202M157I (ATG→ATAyjgNDUF898 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base A (78/111);  total (86/116)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 202 total aligned reads displayed.

TTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTTC  >  CP009273/4469872‑4470162
                                                                                                                                                 |                                                                                                                                                 
ttttGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGa                                                                                                                                               >  2:1045220/1‑150 (MQ=255)
  ttGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGCAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACt                                                                                                                                             <  2:532852/150‑1 (MQ=255)
   tGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTg                                                                                                                                            >  1:2314528/1‑150 (MQ=255)
    gCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGt                                                                                                                                           >  1:4231893/1‑150 (MQ=255)
    gCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGt                                                                                                                                           >  1:5136175/1‑150 (MQ=255)
      ttATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGttt                                                                                                                                         <  1:2924738/150‑1 (MQ=255)
      ttATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGttt                                                                                                                                         <  1:5160677/150‑1 (MQ=255)
       tataTTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGtttt                                                                                                                                        <  1:1649895/150‑1 (MQ=255)
       tataTTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCAGTGGGATAGGGACTGtttt                                                                                                                                        >  2:2768870/1‑150 (MQ=255)
        atatTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTc                                                                                                                                       <  1:496477/150‑1 (MQ=255)
        atatTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTc                                                                                                                                       <  1:496450/150‑1 (MQ=255)
           tttATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCttt                                                                                                                                    >  1:3313759/1‑150 (MQ=255)
           tttATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCttt                                                                                                                                    <  2:3909920/150‑1 (MQ=255)
            ttATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTa                                                                                                                                   >  1:4655852/1‑150 (MQ=255)
              aTGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAATGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTAtc                                                                                                                                 <  1:1607787/150‑1 (MQ=255)
                ggCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTAtctc                                                                                                                               <  1:1364790/150‑1 (MQ=255)
                 gcagcaAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCg                                                                                                                              <  1:2720575/150‑1 (MQ=255)
                         aGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGAt                                                                                                                      <  1:1246899/150‑1 (MQ=255)
                         aGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGAt                                                                                                                      <  1:1246856/150‑1 (MQ=255)
                             cTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTa                                                                                                                  <  2:3271564/150‑1 (MQ=255)
                               gCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGTTAGGGACTGTTTTCTTTATCTCGACAAAGATGCTAcc                                                                                                                <  2:3156946/150‑1 (MQ=255)
                               gCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTAcc                                                                                                                <  1:1139771/150‑1 (MQ=255)
                               gCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTAcc                                                                                                                <  1:2257818/150‑1 (MQ=255)
                                 gTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTg                                                                                                              <  2:4692224/150‑1 (MQ=255)
                                 gTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTg                                                                                                              >  2:4134335/1‑150 (MQ=255)
                                 gTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTg                                                                                                              >  2:5052550/1‑150 (MQ=255)
                                   catcaGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGcc                                                                                                            >  2:1016793/1‑150 (MQ=255)
                                    atcaGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                           <  1:3602105/150‑1 (MQ=255)
                                    atcaGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                           >  1:4695448/1‑150 (MQ=255)
                                      caGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                         >  1:621766/1‑150 (MQ=255)
                                      caGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                         >  1:4370872/1‑150 (MQ=255)
                                       aGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATa                                                                                                        >  2:1535917/1‑150 (MQ=255)
                                          ccTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGtt                                                                                                     >  2:2490431/1‑150 (MQ=255)
                                           cTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTc                                                                                                    >  2:3219769/1‑150 (MQ=255)
                                               aTGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagt                                                                                                >  2:5185418/1‑150 (MQ=255)
                                               aTGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagt                                                                                                <  1:392194/150‑1 (MQ=255)
                                                 gACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagtag                                                                                              <  2:1887916/150‑1 (MQ=255)
                                                    cAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACACAGATGCTACCTGCCAATAGTTCAAGTAgtgt                                                                                           >  2:2844238/1‑150 (MQ=255)
                                                     aGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGtt                                                                                          >  1:866186/1‑150 (MQ=255)
                                                      gTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGtta                                                                                         <  2:1767245/150‑1 (MQ=255)
                                                        ctcAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttatt                                                                                       >  2:4965668/1‑150 (MQ=255)
                                                          cAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTtattat                                                                                     <  1:2803671/150‑1 (MQ=255)
                                                          cAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTtattat                                                                                     >  2:4817641/1‑150 (MQ=255)
                                                             cGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATc                                                                                  >  2:4371275/1‑150 (MQ=255)
                                                               gCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCa                                                                                <  1:4089015/150‑1 (MQ=255)
                                                                cGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCAt                                                                               >  1:2969242/1‑150 (MQ=255)
                                                                cGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCAt                                                                               >  2:79530/1‑150 (MQ=255)
                                                                 gTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATg                                                                              >  1:2995455/1‑150 (MQ=255)
                                                                  tAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATgg                                                                             <  2:3224548/150‑1 (MQ=255)
                                                                   aaGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCACTGGt                                                                            >  1:2856260/1‑150 (MQ=255)
                                                                   aaGATTTAGTTTTAATTGCTCGATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGt                                                                            >  1:908654/1‑150 (MQ=255)
                                                                      aTTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGTCCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGTGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                         <  2:3029399/150‑1 (MQ=255)
                                                                      aTTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                         <  1:3908557/150‑1 (MQ=255)
                                                                       tttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTg                                                                        <  2:1614813/150‑1 (MQ=255)
                                                                          aGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATg                                                                     <  2:641426/150‑1 (MQ=255)
                                                                                 aTTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAg                                                              >  1:4519296/1‑150 (MQ=255)
                                                                                 aTTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAg                                                              >  1:1958746/1‑150 (MQ=255)
                                                                                     ctctATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTgg                                                          >  2:3336611/1‑150 (MQ=255)
                                                                                     ctctATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTgg                                                          <  1:1179753/150‑1 (MQ=255)
                                                                                         aTGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTAtt                                                      <  2:3210993/150‑1 (MQ=255)
                                                                                           gAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGt                                                    <  2:4551367/150‑1 (MQ=255)
                                                                                            aaaGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGtt                                                   <  1:1422433/150‑1 (MQ=255)
                                                                                               tggTTCTGCAGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATAGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGGGAGATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTcc                                                >  2:37126/2‑150 (MQ=255)
                                                                                               gggTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTcc                                                <  2:2322669/150‑1 (MQ=255)
                                                                                                ggTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCa                                               <  2:1834038/150‑1 (MQ=255)
                                                                                                 gTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCAt                                              <  1:2465902/150‑1 (MQ=255)
                                                                                                    cTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTgg                                           <  2:1097541/150‑1 (MQ=255)
                                                                                                    cTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTgg                                           >  1:1874842/1‑150 (MQ=255)
                                                                                                     tggtggGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGt                                          >  1:3675929/1‑150 (MQ=255)
                                                                                                        tgggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAtt                                       >  1:2173818/1‑150 (MQ=255)
                                                                                                        tgggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAtt                                       >  1:4501727/1‑150 (MQ=255)
                                                                                                        tgggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAtt                                       >  1:3767345/1‑150 (MQ=255)
                                                                                                         gggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAttt                                      <  2:4105987/150‑1 (MQ=255)
                                                                                                           gtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAttttt                                    <  2:290704/150‑1 (MQ=255)
                                                                                                            tgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTa                                   >  2:384788/1‑150 (MQ=255)
                                                                                                              aCCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAAt                                 <  1:4580118/150‑1 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                                >  2:628159/1‑150 (MQ=255)
                                                                                                                cTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATgg                               <  2:3049237/150‑1 (MQ=255)
                                                                                                                  ttttCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTa                             >  2:2250789/1‑150 (MQ=255)
                                                                                                                  ttttCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTa                             >  2:2250823/1‑150 (MQ=255)
                                                                                                                      cTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACttt                         <  2:2860443/150‑1 (MQ=255)
                                                                                                                      cTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACttt                         >  1:4833457/1‑150 (MQ=255)
                                                                                                                       ttGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTta                        <  2:4296603/150‑1 (MQ=255)
                                                                                                                       ttGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTta                        <  2:1955586/150‑1 (MQ=255)
                                                                                                                        tGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTtat                       <  1:758370/150‑1 (MQ=255)
                                                                                                                         gCCGATTTTAATGGCCATTGGGAGAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTtata                      <  2:3128244/150‑1 (MQ=255)
                                                                                                                              ttttAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTc                 <  2:3310278/150‑1 (MQ=255)
                                                                                                                               tttAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCt                <  1:1519285/150‑1 (MQ=255)
                                                                                                                               tttAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCt                <  2:4193728/150‑1 (MQ=255)
                                                                                                                               tttAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCt                <  2:2593272/150‑1 (MQ=255)
                                                                                                                                ttAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTg               <  2:866186/150‑1 (MQ=255)
                                                                                                                                   aTGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTa            <  2:493149/150‑1 (MQ=255)
                                                                                                                                        cATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAg       <  2:448685/150‑1 (MQ=255)
                                                                                                                                         aTTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGt      >  1:2760592/1‑150 (MQ=255)
                                                                                                                                         aTTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGt      >  2:307798/1‑150 (MQ=255)
                                                                                                                                         aTTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGt      >  1:5189260/1‑150 (MQ=255)
                                                                                                                                           tGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGttt    <  2:3415177/150‑1 (MQ=255)
                                                                                                                                             ggATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTtc  <  1:3264592/150‑1 (MQ=255)
                                                                                                                                             ggATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTtc  <  1:4694433/150‑1 (MQ=255)
                                                                                                                                             ggATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTtc  <  1:1282997/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTTC  >  CP009273/4469872‑4470162

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: