Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A2 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 443,081 C→T L27L (CTG→CTA cyoE ← heme o synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273443,0810CT98.4% 231.2 / ‑7.2 62L27L (CTG→CTAcyoEheme o synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/30);  minor base A (0/1);  total (31/31)
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.92e-01

GTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCA  >  NZ_CP009273/442950‑443205
                                                                                                                                    |                                                                                                                             
gTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAAt                                                                                                                         <  1:922567/138‑1 (MQ=255)
     ttCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATccccc                                                                                                                    >  1:416761/1‑138 (MQ=255)
       ccATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCaa                                                                                                                  <  1:889673/138‑1 (MQ=255)
         aTCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCa                                                                                                               <  2:16520/139‑1 (MQ=255)
              tCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCAcc                                                                                                             >  2:1147370/1‑136 (MQ=255)
                     atatCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGt                                                                                                   >  1:1258088/1‑139 (MQ=255)
                            tGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGt                                                                                                   >  1:1083889/1‑132 (MQ=255)
                            tGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGt                                                                                                   <  2:1083889/132‑1 (MQ=255)
                                gATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatg                                                                                        <  1:120405/139‑1 (MQ=255)
                                   gTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGatgatg                                                                                     >  2:510731/1‑139 (MQ=255)
                                      ggtgttAAACACACACCACGACGCCACATCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGAGAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCCATCAACGCGATCATGTTGCCAAAGATGATGCCt                                                                                 >  2:284321/3‑139 (MQ=255)
                                        tgttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                                                                               >  2:1021518/1‑139 (MQ=255)
                                        tgttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                                                                               >  1:494491/1‑139 (MQ=255)
                                            aaaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                                                           >  2:183955/1‑139 (MQ=255)
                                                acacaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                                                       <  1:737598/139‑1 (MQ=255)
                                                  acaACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                                                     <  2:305593/139‑1 (MQ=255)
                                                    aaCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTg                                                                   >  2:338116/1‑139 (MQ=255)
                                                     aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                                                  <  2:494491/139‑1 (MQ=255)
                                                         gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt                                                                        >  2:805017/1‑129 (MQ=255)
                                                         gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGtt                                                                        <  1:805017/129‑1 (MQ=255)
                                                         gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                                                              <  2:33386/139‑1 (MQ=255)
                                                         gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                                                              <  2:927716/139‑1 (MQ=255)
                                                         gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                                                              <  1:1147370/139‑1 (MQ=255)
                                                         gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                                                              <  2:711804/139‑1 (MQ=255)
                                                          aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                             >  1:124932/1‑139 (MQ=255)
                                                          aCGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                             <  1:1185268/139‑1 (MQ=255)
                                                              cacaACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGc                                                         >  2:360471/1‑139 (MQ=255)
                                                                caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                                                              >  1:1098799/1‑132 (MQ=255)
                                                                caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGt                                                              <  2:1098799/132‑1 (MQ=255)
                                                                caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                             <  1:752455/133‑1 (MQ=255)
                                                                caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                             >  2:752455/1‑133 (MQ=255)
                                                                 aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                                                      <  2:1174944/139‑1 (MQ=255)
                                                                 aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                                                      <  1:238610/139‑1 (MQ=255)
                                                                 aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                                                      <  1:518570/139‑1 (MQ=255)
                                                                 aaCCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                                                      <  2:805609/139‑1 (MQ=255)
                                                                   ccAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaaa                                                    <  2:117811/139‑1 (MQ=255)
                                                                      gTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcat                                                 >  1:345031/1‑139 (MQ=255)
                                                                           acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatg                                                                                        >  2:1098454/1‑96 (MQ=255)
                                                                           acCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatg                                                                                        <  1:1098454/96‑1 (MQ=255)
                                                                                      cGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGctct                                  >  2:507329/1‑138 (MQ=255)
                                                                                       gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCcagc                                                                                                                               <  1:1039277/45‑1 (MQ=255)
                                                                                       gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCcagc                                                                                                                               >  2:1039277/1‑45 (MQ=255)
                                                                                       gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGctct                                  >  1:627824/1‑137 (MQ=255)
                                                                                       gtgtAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGctct                                  >  2:494319/1‑137 (MQ=255)
                                                                                          tAGAAAAAAAGGGATAAATCATGGCTGCCTTTTGAGCCAGGTAGGACCCCCCCAATCCCCGAGATCGGGTTGCAAAAGATGTTGCTTGTTTTCGTTCTTTGCAGTTTTTTTTTAAACATCATAACCGCCGCTCTTAGTg                             <  2:1258088/139‑1 (MQ=255)
                                                                                                    ggggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                                                  >  2:218039/1‑92 (MQ=255)
                                                                                                    ggggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGc                                                                  <  1:218039/92‑1 (MQ=255)
                                                                                                    ggggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtg                   >  2:897346/1‑139 (MQ=255)
                                                                                                    ggggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtg                   >  1:466784/1‑139 (MQ=255)
                                                                                                    ggggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtg                   >  1:170162/1‑139 (MQ=255)
                                                                                                     gggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtgt                  >  2:783645/1‑139 (MQ=255)
                                                                                                     gggATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtgt                  >  2:1105956/1‑139 (MQ=255)
                                                                                                     gggAAAAACAATGCCTCCCTTTGGGGGCCAGAAGGAACCCCCAAACCACCGAAATCGGGTTCCAAAAGATGTGCCTTGTTTTCTTTCTTTGAAGTTTTTTTTAAAACATAATAACCGCCGCTCTTAGTGCATCATCAtg                   <  2:940469/139‑1 (MQ=255)
                                                                                                       gATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                             >  1:836773/1‑94 (MQ=255)
                                                                                                       gATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                             <  2:836773/94‑1 (MQ=255)
                                                                                                       gATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtgttg                >  1:642243/1‑139 (MQ=255)
                                                                                                          aaTCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAg             >  1:619529/1‑139 (MQ=255)
                                                                                                             catggctccCCTTGGAGGCCAGTAGGACCCCCCCAACCCCCGATATCAGGTGGCAAAAGATGATGCCTGTTTTCTTTCTTGGCAGTTTTTTCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTg          <  1:783645/135‑1 (MQ=255)
                                                                                                             cAATGCTCCCTTTGGAGGCCAGTAGGAACCCCCCAATCCCCGAGATCGGGTTGCAAAAGATGATCCTTGTTTTCTTTTTTGGCGGGTTTTGCTAAAACAACATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTg          <  1:183955/139‑1 (MQ=255)
                                                                                                               aTGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                                                      >  2:386526/1‑93 (MQ=255)
                                                                                                               aTGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTa                                                      <  1:386526/93‑1 (MQ=255)
                                                                                                                  ctgcCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCACCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGt      <  1:943265/138‑1 (MQ=255)
                                                                                                                  ctgcCCTTTGAGGCCAGTAGGAACCCCCCAACCCCCGAGATCAGTTTGCAAAAGATGATGCCTGTTTCCTTTACTTGCGGGTATTGTTAAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGt      <  1:1137314/138‑1 (MQ=255)
                                                                                                                     cccTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCca  >  1:287134/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                             
GTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCA  >  NZ_CP009273/442950‑443205

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAA  >  NZ_CP009273/442933‑443228
                                                                                                                                                    |                                                                                                                                                   
CCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGT                                                                                                                                                     >  SRR3722129.309964/1‑149 (MQ=60)
   GCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAA                                                                                                                                                <  SRR3722129.510220/151‑1 (MQ=60)
   GCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAA                                                                                                                                                <  SRR3722129.1288579/151‑1 (MQ=60)
      CCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCC                                                                                                                                             <  SRR3722129.1069426/151‑1 (MQ=60)
          ttgttttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCC                                                                                                                                          <  SRR3722129.1053891/67‑1 (MQ=60)
           TTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCA                                                                                                                                         >  SRR3722129.422841/1‑150 (MQ=60)
                CGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCAC                                                                                                                                    <  SRR3722129.935686/150‑1 (MQ=60)
                       TCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATC                                                                                                                             <  SRR3722129.902324/150‑1 (MQ=60)
                           TCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTT                                                                                                                        >  SRR3722129.1275583/1‑151 (MQ=60)
                                  GTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAA                                                                                                                   >  SRR3722129.1099064/1‑149 (MQ=60)
                                              GTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGT                                                                                                     >  SRR3722129.501644/1‑151 (MQ=60)
                                                CGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTT                                                                                                   <  SRR3722129.122051/151‑1 (MQ=60)
                                                   cgtgggctcggagatgtgtataagagacagCAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTT                                                                                                   <  SRR3722129.1113783/118‑1 (MQ=60)
                                                               tctttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacaGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGG                                                                                      <  SRR3722129.707579/69‑1 (MQ=60)
                                                               ACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTA                                                                                    <  SRR3722129.816626/151‑1 (MQ=60)
                                                                CACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAT                                                                                   >  SRR3722129.126637/1‑151 (MQ=60)
                                                                CACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAT                                                                                   <  SRR3722129.748398/151‑1 (MQ=60)
                                                                 tttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTAT                                                                                   <  SRR3722129.479182/73‑1 (MQ=60)
                                                                      tagtctcgtgggctcggagatgtgtataagagacagGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT                                                                              <  SRR3722129.221077/114‑1 (MQ=60)
                                                                      ACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT                                                                              >  SRR3722129.1114132/1‑150 (MQ=60)
                                                                         CGACGCCACAACCAGTGACACACCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAC                                                                          <  SRR3722129.1163311/151‑1 (MQ=60)
                                                                          GACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACA                                                                         <  SRR3722129.1201757/151‑1 (MQ=60)
                                                                          GACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACA                                                                         <  SRR3722129.763428/151‑1 (MQ=60)
                                                                            CGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATC                                                                       >  SRR3722129.350006/1‑151 (MQ=60)
                                                                                 CAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAC                                                                  <  SRR3722129.526143/151‑1 (MQ=60)
                                                                                 CAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAC                                                                  <  SRR3722129.242003/151‑1 (MQ=60)
                                                                                  agtctcgtgggctcggagatgtgtataagagacagGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAC                                                                  <  SRR3722129.392176/115‑1 (MQ=60)
                                                                                             CCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTT                                                        >  SRR3722129.637158/1‑149 (MQ=60)
                                                                                                          GTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGT                                         >  SRR3722129.473574/1‑151 (MQ=60)
                                                                                                          GTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGT                                         >  SRR3722129.172502/1‑151 (MQ=60)
                                                                                                             GATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACActgtctcttatacacatctccgagcccacgagact                                      >  SRR3722129.848783/1‑116 (MQ=60)
                                                                                                             GATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGT                                      >  SRR3722129.651822/1‑151 (MQ=60)
                                                                                                                AAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGT                                   >  SRR3722129.628738/1‑151 (MQ=60)
                                                                                                                           AATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCAC                        >  SRR3722129.291253/1‑151 (MQ=60)
                                                                                                                             acaatgctccctttgGAGGCCAGTAGGAACCCCCCAATCCCCGAGATCGGGTTGCAAAAGATGATCCTTGTTTTCTTTTTTGGCGGGTTTTGCTAAAACAACATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACAT                      <  SRR3722129.186544/136‑1 (MQ=50)
                                                                                                                                  GCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCACCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATC                  <  SRR3722129.956615/150‑1 (MQ=60)
                                                                                                                                  GCTGCCCTTTGAGGCCAGTAGGAACCCCCCAACCCCCGAGATCAGTTTGCAAAAGATGATGCCTGTTTCCTTTACTTGCGGGTATTGTTAAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATC                  <  SRR3722129.1153137/150‑1 (MQ=60)
                                                                                                                                                 CAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCtgtctcttatacacatctccgagcccacgagactaaggcgaat  >  SRR3722129.372683/1‑108 (MQ=60)
                                                                                                                                                    |                                                                                                                                                   
CCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAA  >  NZ_CP009273/442933‑443228

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