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
A1 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,786,991 T→C D18G (GAT→GGT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,9910TC100.0% 334.4 / NA 87D18G (GAT→GGT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/61);  total (26/61)

TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATT  >  NZ_CP009273/3786904‑3787076
                                                                                       |                                                                                     
tCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAAccc                                                                                     <  2:2407541/90‑1 (MQ=255)
  aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                   <  2:1557422/90‑1 (MQ=255)
  aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                   >  1:336761/1‑90 (MQ=255)
  aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                   <  2:133337/90‑1 (MQ=255)
   aGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCct                                                                                  <  1:2773150/90‑1 (MQ=255)
    gTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctc                                                                                 <  1:2406724/90‑1 (MQ=255)
     tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                                >  2:1865484/1‑90 (MQ=255)
     tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                                >  2:1963272/1‑90 (MQ=255)
     tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                                >  2:525908/1‑90 (MQ=255)
        ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                   >  2:1524560/1‑84 (MQ=255)
        ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                   <  1:1524560/84‑1 (MQ=255)
         gCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCa                                                                            <  1:1287234/90‑1 (MQ=255)
             gACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGt                                                                        <  1:2822183/90‑1 (MQ=255)
             gACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGt                                                                        >  1:1752374/1‑90 (MQ=255)
               cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                      <  1:635653/90‑1 (MQ=255)
               cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                      <  2:157953/90‑1 (MQ=255)
               cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                      <  2:1604547/90‑1 (MQ=255)
                gCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCaa                                                                     <  1:2598069/90‑1 (MQ=255)
                gCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCaa                                                                     <  1:1325378/90‑1 (MQ=255)
                  tCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCaaaa                                                                   <  2:838742/90‑1 (MQ=255)
                  tCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCaaaa                                                                   <  1:2460198/90‑1 (MQ=255)
                  tCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCaaaa                                                                   <  1:1577486/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  2:2494809/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:762174/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:821289/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:893115/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:947296/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:315598/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:1572397/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:2780623/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  2:767912/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  2:1830151/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  2:2946508/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  2:336761/90‑1 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctc                                                                                 >  2:2329623/1‑69 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctc                                                                                 <  1:2329623/69‑1 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                            <  2:2139021/90‑1 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                            >  2:1703668/1‑90 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                            >  1:1444580/1‑90 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                            <  2:749531/90‑1 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                            >  2:3019074/1‑90 (MQ=255)
                           cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                          <  1:1557257/90‑1 (MQ=255)
                           cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                          <  1:3019074/90‑1 (MQ=255)
                             aaaaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGAt                                                        >  2:610835/1‑90 (MQ=255)
                               aaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATaa                                                      >  2:2356685/1‑90 (MQ=255)
                                aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                     <  1:739267/90‑1 (MQ=255)
                                aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                     <  1:1703668/90‑1 (MQ=255)
                                  aaaCATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATtg                                                   <  1:1963272/90‑1 (MQ=255)
                                     cATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtgg                                                >  1:2665794/1‑90 (MQ=255)
                                      aTCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtggt                                               <  2:2531973/90‑1 (MQ=255)
                                          aCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgat                                           <  2:1752374/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          <  2:2928139/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          <  2:3004239/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          <  2:1750964/90‑1 (MQ=255)
                                              aaCGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGAtt                                       >  1:1089446/1‑90 (MQ=255)
                                               aCGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGAttt                                      >  2:1314809/1‑90 (MQ=255)
                                                cGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTc                                     <  2:2665794/90‑1 (MQ=255)
                                                  aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                   >  2:856804/1‑90 (MQ=255)
                                                  aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                   >  1:2364648/1‑90 (MQ=255)
                                                   gCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAtt                                  <  1:1052635/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                                 <  2:1089446/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                                 <  2:529587/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                                 <  1:1999203/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                                 <  1:2356685/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                                 <  2:2909875/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                                 <  1:2766636/90‑1 (MQ=255)
                                                     tAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAc                                <  1:1791896/90‑1 (MQ=255)
                                                       gTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             >  2:2054822/1‑59 (MQ=255)
                                                       gTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:2054822/59‑1 (MQ=255)
                                                       gTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACta                              >  2:1904986/1‑90 (MQ=255)
                                                          tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                           <  2:2491062/90‑1 (MQ=255)
                                                          tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                           <  1:1314809/90‑1 (MQ=255)
                                                          tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                           <  1:2856646/90‑1 (MQ=255)
                                                            tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAg                         <  2:701944/90‑1 (MQ=255)
                                                             cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGa                        <  1:2626344/90‑1 (MQ=255)
                                                             cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGa                        <  2:1444580/90‑1 (MQ=255)
                                                                 gAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             <  1:1465872/49‑1 (MQ=255)
                                                                 gAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                             >  2:1465872/1‑49 (MQ=255)
                                                                      cTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          >  1:254811/1‑63 (MQ=255)
                                                                      cTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          <  2:254811/63‑1 (MQ=255)
                                                                              gCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtgg                                                >  1:795821/1‑49 (MQ=255)
                                                                              gCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtgg                                                <  2:795821/49‑1 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          <  1:2471379/53‑1 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                          >  2:2471379/1‑53 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa     >  1:1382943/1‑90 (MQ=255)
                                                                                  tCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAAt   >  2:2410118/1‑90 (MQ=255)
                                                                                  tCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAAt   >  2:23492/1‑90 (MQ=255)
                                                                                   ccAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAAtt  <  2:13647/90‑1 (MQ=255)
                                                                                       |                                                                                     
TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATT  >  NZ_CP009273/3786904‑3787076

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATAAATATT  >  NZ_CP009273/3786894‑3787090
                                                                                                 |                                                                                                   
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCC                                                                                                   <  SRR3722030.2171437/100‑1 (MQ=60)
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCC                                                                                                   <  SRR3722030.1534920/100‑1 (MQ=60)
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCC                                                                                                   <  SRR3722030.1453724/100‑1 (MQ=60)
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCC                                                                                                   <  SRR3722030.678569/100‑1 (MQ=60)
  TCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCC                                                                                                 >  SRR3722030.341720/1‑100 (MQ=60)
     CATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCT                                                                                              >  SRR3722030.501577/1‑100 (MQ=60)
       TTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAC                                                                                            <  SRR3722030.545689/100‑1 (MQ=60)
       TTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAC                                                                                            <  SRR3722030.246682/100‑1 (MQ=60)
       TTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAC                                                                                            <  SRR3722030.1520997/100‑1 (MQ=60)
        TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722030.1577999/100‑1 (MQ=60)
        TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722030.866572/100‑1 (MQ=60)
        TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722030.1595498/100‑1 (MQ=60)
        TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722030.2129499/100‑1 (MQ=60)
        TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722030.2100995/100‑1 (MQ=60)
         agagacaGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCA                                                                                          <  SRR3722030.1129026/93‑1 (MQ=60)
         TTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCA                                                                                          <  SRR3722030.1809545/100‑1 (MQ=60)
            AAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATG                                                                                       <  SRR3722030.1549153/100‑1 (MQ=60)
             AGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGT                                                                                      <  SRR3722030.2824373/100‑1 (MQ=60)
             AGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGT                                                                                      >  SRR3722030.1781419/1‑100 (MQ=60)
              ataagagacagCGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTC                                                                                     <  SRR3722030.2371920/89‑1 (MQ=60)
              GTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTC                                                                                     <  SRR3722030.2450575/100‑1 (MQ=60)
                   GCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAAT                                                                                <  SRR3722030.1306682/100‑1 (MQ=60)
                       GACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGC                                                                            <  SRR3722030.2874270/100‑1 (MQ=60)
                         CGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAG                                                                          >  SRR3722030.1467528/1‑100 (MQ=60)
                         CGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAG                                                                          <  SRR3722030.645529/100‑1 (MQ=60)
                          GCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGT                                                                         <  SRR3722030.2645765/100‑1 (MQ=60)
                          GCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGT                                                                         <  SRR3722030.1345732/100‑1 (MQ=60)
                            TCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGA                                                                       <  SRR3722030.2505159/100‑1 (MQ=60)
                            TCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGA                                                                       <  SRR3722030.1603040/100‑1 (MQ=60)
                                  tcagatgtgtataagagacagGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.2091251/79‑1 (MQ=60)
                                  gtcggcagcgtcagatgtgtataagagacagGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.1489247/69‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.320191/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.2831993/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.961764/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.773916/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.833904/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.1597852/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722030.906814/100‑1 (MQ=60)
                                     CTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGG                                                              >  SRR3722030.2714921/1‑100 (MQ=60)
                                     CTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGG                                                              <  SRR3722030.3074672/100‑1 (MQ=60)
                                     CTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGG                                                              <  SRR3722030.1582445/100‑1 (MQ=60)
                                          ATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGAT                                                         <  SRR3722030.750661/100‑1 (MQ=60)
                                          ATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGAT                                                         <  SRR3722030.1731435/100‑1 (MQ=60)
                                            AAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGA                                                       <  SRR3722030.1997591/100‑1 (MQ=60)
                                              ACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATT                                                     >  SRR3722030.1105729/1‑100 (MQ=60)
                                                  CAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAT                                                 >  SRR3722030.2407648/1‑100 (MQ=60)
                                                     gcagcgtcagatgtgtataagagacagCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAC                                              <  SRR3722030.2516535/73‑1 (MQ=60)
                                                             GCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGA                                      <  SRR3722030.1068452/100‑1 (MQ=60)
                                                              CTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC                                     <  SRR3722030.2817730/100‑1 (MQ=60)
                                                              CTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC                                     <  SRR3722030.2399533/100‑1 (MQ=60)
                                                              CTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC                                     <  SRR3722030.2034364/100‑1 (MQ=60)
                                                               TAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACC                                    <  SRR3722030.1821981/100‑1 (MQ=60)
                                                                    TTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAA                               <  SRR3722030.2909401/100‑1 (MQ=60)
                                                                    TTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAA                               <  SRR3722030.1334912/100‑1 (MQ=60)
                                                                      TCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACctgtctcttatacacat                             >  SRR3722030.258485/1‑83 (MQ=60)
                                                                       CACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATAT                            <  SRR3722030.2674649/100‑1 (MQ=60)
                                                                              CCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCtgtctcttatacacatctgacgctgccgac                     >  SRR3722030.808059/1‑70 (MQ=60)
                                                                                CTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA                   >  SRR3722030.1404621/1‑100 (MQ=60)
                                                                                  TTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGctgtctcttatacacatctgacgctgccgacgacatggccagtgtaga                 >  SRR3722030.2398391/1‑52 (MQ=60)
                                                                                                 CCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATAAATATT  >  SRR3722030.2274343/1‑100 (MQ=60)
                                                                                                 |                                                                                                   
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATAAATATT  >  NZ_CP009273/3786894‑3787090

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