Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,333,789 T→G 13.7% Y64D (TAT→GAT)  USA300HOU_1237 → GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,333,7890TG13.7% 143.9 / 3.2 73Y64D (TAT→GAT) USA300HOU_1237GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (23/36);  new base G (11/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 4.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAATTATT  >  CP000730/1333649‑1333926
                                                                                                                                            |                                                                                                                                           
tCTTAGTCGTGGTACATGCTCAGGATTATAAGCAATTTAATTGTATGTCCACAATTTTGGAATTATCATCTTTATCAGAGACTTGCCAGCTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAat                                                                                                                                             <  2:83268/141‑1 (MQ=255)
  ttAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAatat                                                                                                                                           >  1:135839/1‑141 (MQ=255)
  ttAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAatat                                                                                                                                           <  1:33680/141‑1 (MQ=255)
  ttAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAatat                                                                                                                                           <  1:187098/141‑1 (MQ=255)
        ggTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGtt                                                                                                                                   >  1:165244/1‑141 (MQ=255)
         gtgtACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                                                                                                                  >  1:305219/1‑141 (MQ=255)
          tgtACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTgg                                                                                                                                 <  1:434041/141‑1 (MQ=255)
          tgtACATGCTCAAGATGATAAGCAATTTAATTTTGAGCTTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTgg                                                                                                                                 <  2:47306/141‑1 (MQ=255)
           gtACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggt                                                                                                                                >  2:659964/1‑141 (MQ=255)
                 gCTCAAGATGATAAGCAATTTAATTTTGAGTCTACCATGGAGGAATTATCATCTTTATCAGAGACCAGGATACATTAAGAGTTAGGTGTTGGTGCGCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaaggt                                                                                                                          <  1:370475/141‑1 (MQ=255)
                 gCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaaggt                                                                                                                          <  1:533294/141‑1 (MQ=255)
                 gCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaaggt                                                                                                                          >  1:456883/1‑141 (MQ=255)
                          gatAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGAAAAGGTAAAATTgaa                                                                                                                 >  1:305393/1‑141 (MQ=255)
                            tAAGCAATTTAATTTTTAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTgaaga                                                                                                               <  1:770555/141‑1 (MQ=255)
                             aaGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATTTTTATCATAAACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTgaagaa                                                                                                              <  2:232398/141‑1 (MQ=255)
                               gCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTTATGGTGGCATTGGAAAAAATGAAGAAAt                                                                                                            >  2:389514/1‑141 (MQ=255)
                                cAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAAtt                                                                                                           >  1:621237/1‑141 (MQ=255)
                                cAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACATGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAAtt                                                                                                           >  1:466265/1‑141 (MQ=255)
                                 aatttaatttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCGAATGTCACGTTGGAAATGGTAAAATTGAAGAAAtta                                                                                                          >  2:229243/1‑140 (MQ=255)
                                 aatttaatttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGGTGGTAAAGGTAAAATTGAAGAAAtta                                                                                                          >  2:546965/1‑140 (MQ=255)
                                    ttaatttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCTAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAg                                                                                                       <  2:527954/141‑1 (MQ=255)
                                     taatttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGc                                                                                                      <  2:758982/141‑1 (MQ=255)
                                       atttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCAt                                                                                                    >  2:414151/1‑141 (MQ=255)
                                       atttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCCTGCTTATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCAt                                                                                                    >  1:404317/1‑141 (MQ=255)
                                           tGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAGGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGCACACGTCAAATATCATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTAt                                                                                                >  1:175089/1‑141 (MQ=255)
                                            gAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTAtt                                                                                               <  1:511996/141‑1 (MQ=255)
                                              gTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGa                                                                                             <  2:530052/141‑1 (MQ=255)
                                               tCTACAATGGAGGAATTATCATCTTTATCAGAAACTTGCCAACTTGAAGTGTTGGGTCAAATAACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAg                                                                                            >  1:555298/1‑141 (MQ=255)
                                                 tACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCACAGGTTACGTTACTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTAAGtt                                                                                          >  2:345305/1‑141 (MQ=255)
                                                 tACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGtt                                                                                          >  1:280954/1‑141 (MQ=255)
                                                  aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATCTTATGTTGGTAAAGGTAAAATTAAAGAAATTCAAGCATTTATTAAGtta                                                                                         >  1:29023/1‑140 (MQ=255)
                                                   cAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGATTGTGTTGTTATAGAAAAAATTGAAGAAATTCAAACATTTATTGAGttaa                                                                                        >  2:203676/1‑139 (MQ=255)
                                                     aTGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCaaa                                                                                      >  1:669971/1‑141 (MQ=255)
                                                      tGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAg                                                                                     <  2:581978/141‑1 (MQ=255)
                                                        gaggaATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGat                                                                                   <  2:285210/141‑1 (MQ=255)
                                                        gaggaATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGGTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGat                                                                                   >  2:205618/1‑141 (MQ=255)
                                                          ggaATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCCGTATTATAAGAAAAAAAGAAGAAATTCAGGCATTTATTGTGTTCAAAGatat                                                                                 >  2:606006/1‑141 (MQ=255)
                                                            aaTTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTg                                                                               >  1:136614/1‑141 (MQ=255)
                                                             aTTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGa                                                                              <  2:61691/141‑1 (MQ=255)
                                                              ttATCATCTGTATCAGAGACTTGCCAACTTGAAGGGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATTGTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGAt                                                                             <  2:602884/141‑1 (MQ=255)
                                                              ttATCAGCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGAt                                                                             <  2:260097/141‑1 (MQ=255)
                                                               tATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAGATGTCGCCGTAGCATTGAAAAAATAAAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATg                                                                            >  2:108240/1‑141 (MQ=255)
                                                                        tATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTCCTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGtcatca                                                                   <  1:722998/141‑1 (MQ=255)
                                                                           cAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACaa                                                                <  1:659964/141‑1 (MQ=255)
                                                                           cAGAGACATGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCCGTAGCATAAAAAAAAAAAAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCGCaa                                                                >  2:544495/1‑141 (MQ=255)
                                                                               gaCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCGGTAGCAATAGAAAAAATGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAatga                                                            >  2:596779/1‑141 (MQ=255)
                                                                               gaCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGAGATCACCGTGGCAAAGAAAAAAATGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAatga                                                            >  2:116270/1‑141 (MQ=255)
                                                                               gaCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATATAGTCATCACAAatga                                                            <  1:385083/141‑1 (MQ=255)
                                                                                aCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCCGTATCATTAGAAAAATAGAAGAAGTTGAAGGATTTGTTGAGTTCGGAGATATTGATGTAGACATAACAAatgag                                                           >  2:731143/1‑140 (MQ=255)
                                                                                   tGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCCGTGTCATTAGAAAAAAAAAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGGTGTAGTCATGACAAACGATGaa                                                        >  2:178339/1‑141 (MQ=255)
                                                                                    gCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAAt                                                       <  2:305219/141‑1 (MQ=255)
                                                                                     ccAACTTGAAGTGTTGGGTCAAATTACTTAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAACGATATTGATGTTGTCCACACAAATGTTGAAtt                                                      <  2:767794/141‑1 (MQ=255)
                                                                                     ccAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAAtt                                                      >  1:316755/1‑141 (MQ=255)
                                                                                           tGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGAc                                                >  1:564414/1‑141 (MQ=255)
                                                                                              agggtTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGc                                             <  2:34903/138‑1 (MQ=255)
                                                                                              aGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAATATATTGATGTAGTCATCACAAATGATGAATTAACGACTGc                                             <  1:56656/141‑1 (MQ=255)
                                                                                              aGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGc                                             <  2:416642/141‑1 (MQ=255)
                                                                                              aGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAAGGATGAATTAACGACTGc                                             >  1:356517/1‑141 (MQ=255)
                                                                                                 gtTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCGGTGGTCGACGTATCATAAA‑‑AAAAATAAAAAAAATACAAGCATTTATTAAATACAAAGATATCGAGGTAATTATCACAAAAGATGGATTCAcaccaacaac                                          >  2:646036/1‑132 (MQ=255)
                                                                                                    gggTCAAATTACTCAAAACAGGGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAAGTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATc                                       <  2:83327/141‑1 (MQ=255)
                                                                                                           aTTGCTGAAAACGGATATCTTGTATATCTCAAAGATTATGTGGGTAAAGGTAAAAGTGAAGAAATTCAAGCATTTATGGGGTTTAAAGATATTGATGTAGTCATCACAAAGGATGAAGTAACGACTGCACAATCCAAATCa                                <  1:157010/141‑1 (MQ=255)
                                                                                                            ttACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCAc                               <  1:282724/141‑1 (MQ=255)
                                                                                                             tACCCAAACCAGAGATCGTGCAGATCGCAAAAAATATTTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACt                              <  1:272373/141‑1 (MQ=255)
                                                                                                              aCTCAAAACAGAGATCGTGTAGATCGCGGTGGTTGCGGTGGCATTAGTAAAAATGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTa                             >  2:753295/1‑141 (MQ=255)
                                                                                                                  aaaaCAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAAGCCAAATCACTAAATg                         <  1:571925/141‑1 (MQ=255)
                                                                                                                   aaaCAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGa                        <  2:564414/141‑1 (MQ=255)
                                                                                                                     aCAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAg                      <  2:631508/141‑1 (MQ=255)
                                                                                                                        gaaaTCGTGTAGGACGCAAATATTCCGCCTTTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCtt                   <  2:578190/138‑1 (MQ=255)
                                                                                                                        ataagcgagtatctCGCAAAGACTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCACCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCtt                   <  2:582265/128‑1 (MQ=255)
                                                                                                                         agaTCGTGTAGATCGCAAATATTATGTTGGTCCCGGTGTTATGGAAGAACTTCATTCTTTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCttt                  >  1:765074/1‑141 (MQ=255)
                                                                                                                           aTCGTTTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAg                <  1:619957/141‑1 (MQ=255)
                                                                                                                           aTCGTGTAGACCGCAAATAGTATGTTGGTAGAGGTGAAATTGAAGAAATTCCAGCATTTATTGAGTTCAAAGATATTGAGGTAGTGATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAg                <  1:210818/141‑1 (MQ=255)
                                                                                                                               tgtAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGgtgt            <  1:424604/141‑1 (MQ=255)
                                                                                                                                ctAGATCGAAAATACTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTa           <  1:1389/140‑1 (MQ=255)
                                                                                                                                   gATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCATAATCCAAATCACTAAATGAAGCTTTAGGTGTaaaa        <  1:414151/141‑1 (MQ=255)
                                                                                                                                     tCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAat      <  1:205618/141‑1 (MQ=255)
                                                                                                                                         aaaTATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAattatt  <  2:280954/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
TCTTAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAAAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAATTATT  >  CP000730/1333649‑1333926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: