Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 826,969 G→T 6.3% A36D (GCT→GAT)  USA300HOU_0775 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730826,9690GT 6.3% 214.2 / 2.4 72A36D (GCT→GAT) USA300HOU_0775hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base T (3/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGGGACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAGATTATT  >  CP000730/826833‑827108
                                                                                                                                        |                                                                                                                                           
tGGGACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGcttc                                                                                                                                         <  2:701507/141‑1 (MQ=255)
  ggACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGcttcat                                                                                                                                       <  1:408481/141‑1 (MQ=255)
  ggACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGcttcat                                                                                                                                       <  1:470944/141‑1 (MQ=255)
   gACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGTACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTATcttcatt                                                                                                                                      >  1:513191/1‑140 (MQ=255)
   gACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTATcttcatc                                                                                                                                      >  1:513571/1‑141 (MQ=255)
     cacCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGcttcatctt                                                                                                                                    <  2:606731/141‑1 (MQ=255)
      acCGGCAGTTCCTGTCGGTTCACCATCATCCTTAGCCTTTTGAATGTTCAGTTAAGGTCCAATAGTATATGCAGAAGATTTATGAGTTGCATCTTTATGTACTCTTTTTATCGCAGCAATAAATGCTTTAGcttcatcttc                                                                                                                                   >  2:66464/1‑141 (MQ=255)
          gCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCTGCAATAAATGCTTTAGCTTCATCTTCAttt                                                                                                                               >  2:605059/1‑141 (MQ=255)
          ccAGTCCCTGTCGCTTCACCATCATCATTCGCTTTTTGAATACTCCGTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCAttt                                                                                                                               <  1:200259/140‑1 (MQ=255)
           cAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTGGCTTCATCTTCAtttt                                                                                                                              >  1:316569/1‑141 (MQ=255)
            ggTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCCGCACCAAATGCTTTAGCTTCATCTTCATTTTg                                                                                                                             <  1:755790/140‑1 (MQ=255)
            aGTTCCTGTTTGTTTAACATCATAATTCGCTTTATGAATATTCATTTCATGTCCAATAGTATATGCAGAAAAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTg                                                                                                                             <  2:725678/141‑1 (MQ=255)
               tCCTGTCGGTTCACCATCATCATTCGCTTTTTTTTTTTTTTTTTCAAAGCCAATAGTATATGCAGAATAATTATGAGTGGCATGTTTATGTTTTTTTTTTATTCCAGTATTAATTGTTTTATCTTCATATTCTTTTTTAAc                                                                                                                          >  2:323944/1‑141 (MQ=255)
               tCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAAc                                                                                                                          <  1:73563/141‑1 (MQ=255)
               tCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAAc                                                                                                                          <  1:489375/141‑1 (MQ=255)
                ccTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGGTTCATCTTCATTTTGAACa                                                                                                                         <  1:454499/141‑1 (MQ=255)
                   gTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTTCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACCGGt                                                                                                                      >  2:225447/1‑141 (MQ=255)
                   gTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATTTTCATTTTGAACacgg                                                                                                                      >  2:262792/1‑138 (MQ=255)
                   gTCGGTTCACCATCATAATTCGCTTTTTCAATATTCATGTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGt                                                                                                                      <  1:14943/141‑1 (MQ=255)
                     cGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGttt                                                                                                                    <  2:507620/141‑1 (MQ=255)
                        ttCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGat                                                                                                                 <  1:466393/141‑1 (MQ=255)
                                   catTTGCTTTTTGAATATTCAGTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGa                                                                                                      <  1:2114/141‑1 (MQ=255)
                                       cGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTTAACAGGTTTGATATTAGCCATGATTCt                                                                                                  >  2:684914/1‑141 (MQ=255)
                                       cGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTACGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCt                                                                                                  <  1:758962/141‑1 (MQ=255)
                                           ttttGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATACGAGCAATGAATCTTGAt                                                                                              <  2:746780/141‑1 (MQ=255)
                                                    ttCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATc                                                                                     <  2:567556/141‑1 (MQ=255)
                                                     tCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATca                                                                                    <  2:710850/141‑1 (MQ=255)
                                                            aGGTCCAATGGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTTATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTc                                                                             <  2:230468/141‑1 (MQ=255)
                                                              gTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACCTTTTCaa                                                                           >  2:619104/1‑141 (MQ=255)
                                                              ctcccacggTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTCATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCaa                                                                           <  2:410153/133‑1 (MQ=255)
                                                                 cAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCaataa                                                                        >  1:116255/1‑141 (MQ=255)
                                                                    tAGTATATGCAGAACAATTATGATTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAAtgt                                                                     <  1:36489/141‑1 (MQ=255)
                                                                    tAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTACAGCAATAAATGCTTTAGCTTCATCTTCATTTGGAACAAGGTTGATATGAGCAATGAATCTTCATTTACTTATCCCATCTTCaagaatat                                                                     >  2:512612/1‑139 (MQ=255)
                                                                      ttataTGCACAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTCCATTTTGAACAGGTTTGATATGATCAATGAATCTTGATTTACTAATCACATTTTCAATAATgtgt                                                                   <  1:602471/140‑1 (MQ=255)
                                                                       tataTGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGtt                                                                  >  2:616998/1‑141 (MQ=255)
                                                                          atGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCtt                                                               >  1:571359/1‑141 (MQ=255)
                                                                               gAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGACTCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTa                                                          >  2:509208/1‑141 (MQ=255)
                                                                                aaCAATTATGAGTGGCATCTTTATGTTCTTTTTTTTTTTTTTTATTAAACACTTTAGCTTCATCTTCATTTTCAACAGGTTTTATATGAGCAATGAATCTTTTTTTACTCATCACATTTTCAATAATGTGTTCTTTCTTaa                                                         >  2:268516/1‑141 (MQ=255)
                                                                                aaCAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTTCATCAATAAAAGCATTAGCTTCATCTTCAATTCGAACAGGTTTGATATGAGCAATGAATCTTGATCTACTAATCACACTTTCAATAATGTGTTCGTTTTTaa                                                         >  2:207262/1‑141 (MQ=255)
                                                                                    aTTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGt                                                     >  1:65964/1‑141 (MQ=255)
                                                                                     ttATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTa                                                    >  1:25180/1‑141 (MQ=255)
                                                                                     ttATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTa                                                    >  1:265552/1‑141 (MQ=255)
                                                                                     ttATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTa                                                    >  1:302312/1‑141 (MQ=255)
                                                                                      tATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTaa                                                   >  1:162411/1‑141 (MQ=255)
                                                                                      tATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTaa                                                   <  2:117248/141‑1 (MQ=255)
                                                                                       aTGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAAt                                                  >  1:84863/1‑141 (MQ=255)
                                                                                        tGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATACTGTGTTCTTTTTTAACAGTAATg                                                 <  2:736675/141‑1 (MQ=255)
                                                                                        tGAGTGGCATCTTGATGTTCTTTTTTTATGGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATg                                                 <  1:421155/141‑1 (MQ=255)
                                                                                          agcggCATTTTTATGTTCTTTTTTTTTTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGat                                               <  2:630199/138‑1 (MQ=255)
                                                                                          aGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGat                                               <  2:240969/141‑1 (MQ=255)
                                                                                                  cTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTtgtg                                       <  1:372294/141‑1 (MQ=255)
                                                                                                   tttATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTgtgt                                      <  2:652458/141‑1 (MQ=255)
                                                                                                       tGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCata                                  >  1:116761/1‑141 (MQ=255)
                                                                                                            tttttttATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCATATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAAc                             <  1:506769/141‑1 (MQ=255)
                                                                                                            ttttttcATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTTATTTACTTATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAAc                             <  1:614142/141‑1 (MQ=255)
                                                                                                              tttttATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTc                           >  1:444125/1‑141 (MQ=255)
                                                                                                              tttttATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTc                           <  2:573204/141‑1 (MQ=255)
                                                                                                                  tATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTTACATTATTTATATTTTGTTTCTTTATAACTCCTTa                       >  2:435519/1‑141 (MQ=255)
                                                                                                                  tATTGCAGCAATAAATGCTTTAGCTTACTCCGAATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATCATGTGTTCTTTTTTATCAGTGATGATATTTTGTTTAATAATAACCCCTTa                       >  2:196366/1‑141 (MQ=255)
                                                                                                                   aTTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTaa                      >  1:353637/1‑141 (MQ=255)
                                                                                                                    tttcagcaATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCAAATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAAt                     <  2:415706/138‑1 (MQ=255)
                                                                                                                    ttGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAAt                     <  1:277711/141‑1 (MQ=255)
                                                                                                                              aaaTGCTTTATTTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAACTCATAAGCtt           <  2:513571/141‑1 (MQ=255)
                                                                                                                              aaaTGCTTTATCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCtt           <  2:513191/141‑1 (MQ=255)
                                                                                                                               aaTGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTa          >  1:698070/1‑141 (MQ=255)
                                                                                                                                 tGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTaaa        >  1:313154/1‑140 (MQ=255)
                                                                                                                                 tGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAg        >  1:269392/1‑141 (MQ=255)
                                                                                                                                  gCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAGa       <  1:78823/141‑1 (MQ=255)
                                                                                                                                   cTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATACTACCTCCTTAATTCATAAGCTTAAGat      >  2:214691/1‑141 (MQ=255)
                                                                                                                                   cTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAGat      >  2:137072/1‑141 (MQ=255)
                                                                                                                                      tAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAGattat   <  2:116255/141‑1 (MQ=255)
                                                                                                                                       aGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAGattatt  >  1:144159/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                           
TGGGACACCGGCAGTTCCTGTCGGTTCACCATCATCATTCGCTTTTTGAATATTCATTTCAGGTCCAATAGTATATGCAGAACAATTATGAGTGGCATCTTTATGTTCTTTTTTTATTGCAGCAATAAATGCTTTAGCTTCATCTTCATTTTGAACAGGTTTGATATGAGCAATGAATCTTGATTTACTAATCACATTTTCAATAATGTGTTCTTTTTTAACAGTAATGATATTTTGTGTCATAATAACTCCTTAATTCATAAGCTTAAGATTATT  >  CP000730/826833‑827108

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: