Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,466,283 A→T 5.7% N344Y (AAC→TAC)  yjgL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,466,2830AT 5.7% 446.7 / 13.6 123N344Y (AAC→TAC) yjgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (52/64);  new base T (3/4);  total (55/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 123 total aligned reads displayed.

TTAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAATA  >  CP009273/4466136‑4466431
                                                                                                                                                   |                                                                                                                                                    
ttAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAAc                                                                                                                                                    >  2:990574/1‑150 (MQ=255)
ttAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAAc                                                                                                                                                    >  2:2138257/1‑150 (MQ=255)
 tAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAAcc                                                                                                                                                   >  1:3780968/1‑150 (MQ=255)
   aTGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTg                                                                                                                                                 <  1:892995/150‑1 (MQ=255)
      ccTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGtct                                                                                                                                              >  1:3965810/1‑150 (MQ=255)
        tCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCa                                                                                                                                            >  2:897070/1‑150 (MQ=255)
               tgtgATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCtt                                                                                                                                     >  2:2382526/1‑150 (MQ=255)
               tgtgATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCtt                                                                                                                                     <  1:356301/150‑1 (MQ=255)
                gtgATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATTACATCAAAGTTAACCTGTCTCATATTGTTa                                                                                                                                    >  2:2577595/1‑150 (MQ=255)
                   atgaATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTg                                                                                                                                 <  2:1769198/150‑1 (MQ=255)
                    tgaATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGa                                                                                                                                <  1:2882481/150‑1 (MQ=255)
                     gaATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGAt                                                                                                                               >  1:4029054/1‑150 (MQ=255)
                      aaTCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGAtt                                                                                                                              >  1:4695264/1‑150 (MQ=255)
                       aTCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTc                                                                                                                             <  2:2213737/150‑1 (MQ=255)
                          gATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCtt                                                                                                                          <  2:3390756/150‑1 (MQ=255)
                           aTTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTa                                                                                                                         <  2:1066674/150‑1 (MQ=255)
                             ttCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAga                                                                                                                       <  1:2938635/150‑1 (MQ=255)
                             ttCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAga                                                                                                                       <  1:2939374/150‑1 (MQ=255)
                                 tacttactGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAgatgat                                                                                                                   <  2:2475572/150‑1 (MQ=255)
                                 tacttactGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAgatgat                                                                                                                   <  2:1371654/150‑1 (MQ=255)
                                      actGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCaaa                                                                                                              <  2:2913055/150‑1 (MQ=255)
                                      actGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCaaa                                                                                                              <  2:2748170/150‑1 (MQ=255)
                                      actGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCaaa                                                                                                              <  2:2913780/150‑1 (MQ=255)
                                              aaaaaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAAcc                                                                                                      >  2:1082738/1‑150 (MQ=255)
                                              aaaaaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAAcc                                                                                                      >  2:1082375/1‑150 (MQ=255)
                                               aaaaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCt                                                                                                     >  2:2207192/1‑150 (MQ=255)
                                                 aaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTa                                                                                                   >  2:4274263/1‑150 (MQ=255)
                                                 aaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTa                                                                                                   <  1:3903488/150‑1 (MQ=255)
                                                 aaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTa                                                                                                   >  1:455306/1‑150 (MQ=255)
                                                 aaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTa                                                                                                   >  1:456151/1‑150 (MQ=255)
                                                 aaGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTa                                                                                                   <  1:4749477/150‑1 (MQ=255)
                                                          ttGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCa                                                                                          <  2:4505796/150‑1 (MQ=255)
                                                           tGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAt                                                                                         >  1:3645163/1‑150 (MQ=255)
                                                           tGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAt                                                                                         >  2:648410/1‑150 (MQ=255)
                                                            gaagaaATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtc                                                                                        <  1:1455664/150‑1 (MQ=255)
                                                            gaagaaATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtc                                                                                        <  2:3235961/150‑1 (MQ=255)
                                                             aagaaATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcg                                                                                       <  2:4107954/150‑1 (MQ=255)
                                                                 aaTAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCAtcgattc                                                                                   <  2:3925382/150‑1 (MQ=255)
                                                                      aTAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCgattcgattc                                                                              <  1:2856776/150‑1 (MQ=255)
                                                                       tAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCt                                                                             >  1:275556/1‑150 (MQ=255)
                                                                        aaGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTg                                                                            <  1:3583797/150‑1 (MQ=255)
                                                                        aaGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTg                                                                            <  1:2606130/150‑1 (MQ=255)
                                                                          gCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGga                                                                          <  2:1384834/150‑1 (MQ=255)
                                                                             ctACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGAc                                                                       >  1:63957/1‑150 (MQ=255)
                                                                             ctACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGAc                                                                       >  1:2801918/1‑150 (MQ=255)
                                                                              tACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACt                                                                      >  2:4535521/1‑150 (MQ=255)
                                                                              tACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACt                                                                      >  2:4535462/1‑150 (MQ=255)
                                                                               aCGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACtt                                                                     <  2:1639228/150‑1 (MQ=255)
                                                                                cGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACttt                                                                    <  1:3631499/150‑1 (MQ=255)
                                                                                   aTTAAATGATCACGGAATTGATAATGACACAGGAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTc                                                                 <  2:3867816/150‑1 (MQ=255)
                                                                                   aTTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTc                                                                 >  2:2652888/1‑150 (MQ=255)
                                                                                   aTTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTc                                                                 >  2:2652978/1‑150 (MQ=255)
                                                                                     tAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCaa                                                               >  2:3633650/1‑150 (MQ=255)
                                                                                        aTGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAAt                                                            <  2:3780481/150‑1 (MQ=255)
                                                                                        aTGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAAt                                                            <  2:3780968/150‑1 (MQ=255)
                                                                                          gATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATcc                                                          >  1:127945/1‑150 (MQ=255)
                                                                                             cACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCa                                                       <  2:456151/150‑1 (MQ=255)
                                                                                             cACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCa                                                       <  2:455306/150‑1 (MQ=255)
                                                                                               cGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCata                                                     >  2:3213518/1‑150 (MQ=255)
                                                                                                   aTTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATAtt                                                 <  2:4020940/150‑1 (MQ=255)
                                                                                                   aTTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCATTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATAtt                                                 >  1:1806430/1‑150 (MQ=255)
                                                                                                    ttGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTa                                                >  2:1162362/1‑150 (MQ=255)
                                                                                                     tGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTaa                                               >  2:2906689/1‑150 (MQ=255)
                                                                                                       aTAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatg                                             >  1:485574/1‑150 (MQ=255)
                                                                                                        tAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatga                                            <  1:65157/150‑1 (MQ=255)
                                                                                                          aTGACAAAGTAGGTAACAATGAAATTAATGATATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatgatg                                          >  2:3971542/1‑150 (MQ=255)
                                                                                                          aTGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatgatg                                          <  2:3147283/150‑1 (MQ=255)
                                                                                                          aTGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAatgatg                                          >  2:3971533/1‑150 (MQ=255)
                                                                                                           tGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGt                                         <  1:4244830/150‑1 (MQ=255)
                                                                                                            gACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTa                                        <  2:2338629/150‑1 (MQ=255)
                                                                                                            gACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTa                                        <  2:2338147/150‑1 (MQ=255)
                                                                                                             aCAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTaa                                       >  2:3724431/1‑150 (MQ=255)
                                                                                                              cAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAg                                      >  2:2192963/1‑150 (MQ=255)
                                                                                                                 aGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAAt                                   <  2:2913921/150‑1 (MQ=255)
                                                                                                                 aGTAGGTAACAATGAAATTAATGACATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAAt                                   <  2:540977/150‑1 (MQ=255)
                                                                                                                       tAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGa                             <  1:4595777/150‑1 (MQ=255)
                                                                                                                        aaCAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGAt                            <  1:2531929/150‑1 (MQ=255)
                                                                                                                           aaTGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTgg                         >  1:1048010/1‑150 (MQ=255)
                                                                                                                            aTGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGgt                        >  1:3184374/1‑150 (MQ=255)
                                                                                                                            aTGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGgt                        >  1:3184705/1‑150 (MQ=255)
                                                                                                                             tGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGgtg                       >  1:937833/1‑150 (MQ=255)
                                                                                                                               aaaTTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGtt                     <  2:3965810/150‑1 (MQ=255)
                                                                                                                                 aTTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGtttt                   <  1:1203398/150‑1 (MQ=255)
                                                                                                                                   tAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaa                 >  1:3170187/1‑150 (MQ=255)
                                                                                                                                   tAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaa                 >  1:3170088/1‑150 (MQ=255)
                                                                                                                                    aacgACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaat                <  1:2482244/147‑1 (MQ=255)
                                                                                                                                     aTGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaata               <  1:2053465/150‑1 (MQ=255)
                                                                                                                                     aTGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaata               <  1:4434431/150‑1 (MQ=255)
                                                                                                                                     aTGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaata               <  1:44576/150‑1 (MQ=255)
                                                                                                                                     aTGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaata               <  1:44744/150‑1 (MQ=255)
                                                                                                                                      tGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTtaataa              >  2:4371812/1‑150 (MQ=255)
                                                                                                                                       gACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATaaa             <  2:932796/150‑1 (MQ=255)
                                                                                                                                         cATCAAAGTTTACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAg           >  2:3730449/1‑150 (MQ=255)
                                                                                                                                         cATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTACTGATGTAAGAATACTGGATTGGTGTTTTTATAAAAg           >  1:1787844/1‑150 (MQ=255)
                                                                                                                                          aTCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGc          >  2:4052895/1‑150 (MQ=255)
                                                                                                                                            cAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCAt        >  1:375147/1‑150 (MQ=255)
                                                                                                                                            cAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCAt        <  1:2207192/150‑1 (MQ=255)
                                                                                                                                             aaaGTTAAACTGTCTCATATTATTATTGATTACTTTGATGATGCAAATGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATg       <  1:1569606/150‑1 (MQ=255)
                                                                                                                                                  tAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAata  <  2:2706913/150‑1 (MQ=255)
                                                                                                                                                  tAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAata  >  2:54653/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTAATGCCTCATGTATGTGATGAATCGATTTCATACTTACTGGGCGAAAAAGGTATACTTGAAGAAATAGATAAGCTCTACGCATTAAATGATCACGGAATTGATAATGACAAAGTAGGTAACAATGAAATTAATGACATCAAAGTTAACCTGTCTCATATTCTTATTGATTCCTTAGATGATGCAAAGGTTAACCTTACACCGGTCATCGATTCGATTCTGGAGACTTTTTCAAAATCCCCATATATTAATGATGTAAGAATACTGGATTGGTGTTTTAATAAAAGCATGCAATA  >  CP009273/4466136‑4466431

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: