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,840,691 T→G N432H (AAC→CAC)  uhpC ← MFS transporter family glucose‑6‑phosphate receptor UhpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,840,6910TG98.1% 384.7 / ‑1.3 103N432H (AAC→CAC) uhpCMFS transporter family glucose‑6‑phosphate receptor UhpC
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (39/62);  total (40/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

Only 100 of 103 total aligned reads displayed.

CCAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTTCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC  >  NZ_CP009273/3840604‑3840778
                                                                                       |                                                                                       
ccAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCa                                                                                       >  2:573573/1‑90 (MQ=255)
  aGGCAAAACCTGGAAAGTTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaa                                                                                     >  1:2247401/1‑90 (MQ=255)
  aGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaa                                                                                     >  2:2606861/1‑90 (MQ=255)
   ggCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaaa                                                                                    >  2:373461/1‑90 (MQ=255)
   ggCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaaa                                                                                    <  1:2148074/90‑1 (MQ=255)
    gCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaaaa                                                                                   <  2:1770743/90‑1 (MQ=255)
       aaaCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAgggg                                                                                <  1:1466139/90‑1 (MQ=255)
       aaaCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAggg                                                                                <  2:836176/90‑1 (MQ=255)
       aaaCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAggg                                                                                <  1:1208536/90‑1 (MQ=255)
       aaaCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAggg                                                                                <  2:2964437/90‑1 (MQ=255)
       aaaCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAggg                                                                                <  1:402677/90‑1 (MQ=255)
        aaCCTGGATAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGc                                                                               >  2:253370/1‑90 (MQ=255)
        aaCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGc                                                                               >  2:447309/1‑90 (MQ=255)
         aCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTTCAAAAAGGGCa                                                                              <  2:1160877/90‑1 (MQ=255)
         aCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  1:2228576/90‑1 (MQ=255)
         aCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  2:2017502/90‑1 (MQ=255)
         aCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  1:2977884/90‑1 (MQ=255)
         aCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  2:1324782/90‑1 (MQ=255)
         aCCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  1:1440687/90‑1 (MQ=255)
            tGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTa                                                                           <  2:2220892/90‑1 (MQ=255)
            tGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTa                                                                           <  2:689089/90‑1 (MQ=255)
            tGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACACTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTa                                                                           <  1:1855163/90‑1 (MQ=255)
             ggAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTaa                                                                          >  2:684587/1‑90 (MQ=255)
               aaaaTTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAAca                                                                        <  1:1305190/90‑1 (MQ=255)
                 aaTTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagc                                                                      <  1:1895399/90‑1 (MQ=255)
                 aaTTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagc                                                                      <  1:62336/90‑1 (MQ=255)
                   tttCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagcag                                                                    <  1:1056575/90‑1 (MQ=255)
                    ttCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGt                                                                   <  1:826204/90‑1 (MQ=255)
                          gTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  2:2118013/73‑1 (MQ=255)
                          gTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              >  1:2118013/1‑73 (MQ=255)
                             ttGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGt                                                                            <  1:931772/72‑1 (MQ=255)
                                  ggATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggc                                                     <  2:113059/90‑1 (MQ=255)
                                  ggATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggc                                                     <  2:3000072/90‑1 (MQ=255)
                                   gATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggcg                                                    <  1:842009/90‑1 (MQ=255)
                                    atatGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              >  2:923846/1‑63 (MQ=255)
                                    atatGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCa                                                                              <  1:923846/63‑1 (MQ=255)
                                     tatGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCggcggc                                                  <  2:474535/90‑1 (MQ=255)
                                       tGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGa                                                <  1:2363256/90‑1 (MQ=255)
                                         aaGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATa                                              >  1:491545/1‑90 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaaa                                                                                    >  1:1596043/1‑48 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCaaaa                                                                                    <  2:1596043/48‑1 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggcg                                                    <  1:388906/80‑1 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggcg                                                    >  2:388906/1‑80 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAgaga                                          >  1:221373/1‑90 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAgaga                                          >  2:1782573/1‑90 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAgaga                                          >  2:612456/1‑90 (MQ=255)
                                             gAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAgaga                                          >  2:580150/1‑90 (MQ=255)
                                              aaaaGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGAt                                         <  1:1319307/90‑1 (MQ=255)
                                               aaaGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATa                                        <  1:253370/90‑1 (MQ=255)
                                                  ggTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAAcc                                     <  2:2582469/90‑1 (MQ=255)
                                                  ggTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAAcc                                     >  1:319261/1‑90 (MQ=255)
                                                          gCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACaaaaaa                             >  1:2973815/1‑90 (MQ=255)
                                                           cATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAAt                            >  2:2368373/1‑90 (MQ=255)
                                                            aTCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATc                           >  1:1963221/1‑90 (MQ=255)
                                                            aTCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATc                           >  1:923772/1‑90 (MQ=255)
                                                            aTCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATc                           >  1:2904878/1‑90 (MQ=255)
                                                             tCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATcc                          <  1:373461/90‑1 (MQ=255)
                                                              cACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCg                         <  1:684587/90‑1 (MQ=255)
                                                               aCGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagca                                                                     >  2:1060946/1‑45 (MQ=255)
                                                               aCGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagca                                                                     <  1:1060946/45‑1 (MQ=255)
                                                               aCGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGc                        >  2:2079175/1‑90 (MQ=255)
                                                                 gCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTc                      >  2:212167/1‑90 (MQ=255)
                                                                 gCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTc                      <  1:2197309/90‑1 (MQ=255)
                                                                  cTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGt                                                                   <  2:1682395/44‑1 (MQ=255)
                                                                  cTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGt                                                                   >  1:1682395/1‑44 (MQ=255)
                                                                  cTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTcc                     <  1:2685251/90‑1 (MQ=255)
                                                                  cTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTcc                     <  1:307117/90‑1 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATag                                             <  2:3003980/64‑1 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATag                                             >  1:3003980/1‑64 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   >  1:1843362/1‑90 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   <  1:2368373/90‑1 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   <  1:164207/90‑1 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   >  1:1280002/1‑90 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   <  2:710926/90‑1 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   >  1:1107203/1‑90 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   <  2:1107203/90‑1 (MQ=255)
                                                                    tCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAg                   >  1:694068/1‑90 (MQ=255)
                                                                     cgcgcgGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagc                                                                      >  1:1975189/1‑38 (MQ=38)
                                                                     cgcgcgGTGTCTGGGCGTGCAAAAAGGGCAGTAAcagc                                                                      <  2:1975189/38‑1 (MQ=38)
                                                                     cgcgcgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGt                  >  1:629532/1‑90 (MQ=255)
                                                                      gcgcgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTg                 >  2:1711316/1‑90 (MQ=255)
                                                                        gcgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGcc               <  1:2735047/90‑1 (MQ=255)
                                                                        gcgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGcc               <  1:612456/90‑1 (MQ=255)
                                                                        gcgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGcc               <  1:580150/90‑1 (MQ=255)
                                                                        gcgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGcc               <  1:212167/90‑1 (MQ=255)
                                                                         cgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCa              <  2:2797172/90‑1 (MQ=255)
                                                                         cgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCa              <  2:1992798/90‑1 (MQ=255)
                                                                         cgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCa              <  1:1711316/90‑1 (MQ=255)
                                                                         cgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCa              <  2:2820312/90‑1 (MQ=255)
                                                                         cgGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCAAGTGCCa              <  1:2905662/90‑1 (MQ=255)
                                                                           gtgtCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAgg            <  1:1652623/90‑1 (MQ=255)
                                                                                  ggCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAg     <  2:433665/90‑1 (MQ=255)
                                                                                  ggCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAg     <  1:573573/90‑1 (MQ=255)
                                                                                  ggCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAg     <  1:586937/90‑1 (MQ=255)
                                                                                  ggCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAg     <  1:2838818/90‑1 (MQ=255)
                                                                                   gCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggc                                                     >  2:549705/1‑41 (MQ=255)
                                                                                   gCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCcggc                                                     <  1:549705/41‑1 (MQ=255)
                                                                                   gCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGt    >  2:2560291/1‑90 (MQ=255)
                                                                                    cGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTa   <  2:2088367/90‑1 (MQ=255)
                                                                                     gTTCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAc  >  2:2056674/1‑90 (MQ=255)
                                                                                       |                                                                                       
CCAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTTCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC  >  NZ_CP009273/3840604‑3840778

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATAGCGTCCAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTTCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTTTCGCCAGCG  >  NZ_CP009273/3840596‑3840789
                                                                                               |                                                                                                  
GATAGCGTCCAGGCAAAACCTGGAAAGTTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAA                                                                                                >  SRR3722030.2287971/1‑100 (MQ=60)
    GCGTCCAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGG                                                                                            <  SRR3722030.1386011/100‑1 (MQ=60)
     CGTCCAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGG                                                                                           >  SRR3722030.3019082/1‑100 (MQ=60)
           GGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAA                                                                                     <  SRR3722030.2186501/100‑1 (MQ=60)
               AAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAGGGGCAGTAACAGC                                                                                 <  SRR3722030.1489518/100‑1 (MQ=60)
               AAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGC                                                                                 <  SRR3722030.408742/100‑1 (MQ=60)
               AAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGC                                                                                 <  SRR3722030.1226397/100‑1 (MQ=60)
                 atgtgtataagagacagGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAG                                                                               <  SRR3722030.937996/83‑1 (MQ=60)
                 ACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAG                                                                               <  SRR3722030.3032746/100‑1 (MQ=60)
                 ACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAG                                                                               <  SRR3722030.1463549/100‑1 (MQ=60)
                 ACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAG                                                                               <  SRR3722030.2268741/100‑1 (MQ=60)
                   agagacagTTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTG                                                                             <  SRR3722030.946033/92‑1 (MQ=60)
                    TGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACACTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGC                                                                            <  SRR3722030.1886849/100‑1 (MQ=60)
                       AAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGA                                                                         <  SRR3722030.1325051/100‑1 (MQ=60)
                        AAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGctgtctc                                                                        >  SRR3722030.2155802/1‑93 (MQ=60)
                         AATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAA                                                                       <  SRR3722030.63293/100‑1 (MQ=60)
                         AATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAA                                                                       <  SRR3722030.1928099/100‑1 (MQ=60)
                          tgtcgtcggcagcgtcagatgtgtataagagacagGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAAT                                                                      <  SRR3722030.1076877/65‑1 (MQ=60)
                           TTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATC                                                                     <  SRR3722030.1072449/100‑1 (MQ=60)
                            TTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCC                                                                    <  SRR3722030.838892/100‑1 (MQ=60)
                                       GCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATA                                                         >  SRR3722030.499077/1‑100 (MQ=60)
                                          gccatgtcgtcggcagcgtcagatgtgtataagagacaGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAG                                                      <  SRR3722030.558219/62‑1 (MQ=60)
                                           GATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACtgtctcttatacacatctgacgctgccgacg                                                     >  SRR3722030.1621933/1‑69 (MQ=60)
                                           GATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGA                                                     <  SRR3722030.394728/100‑1 (MQ=60)
                                           GATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGA                                                     >  SRR3722030.224534/1‑100 (MQ=60)
                                           GATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGA                                                     <  SRR3722030.854905/100‑1 (MQ=60)
                                             ggccatgtcgtcggcagcgtcagatgtgtataagagacagGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATA                                                   <  SRR3722030.379120/60‑1 (MQ=60)
                                               TGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAAC                                                 <  SRR3722030.2406228/100‑1 (MQ=60)
                                                GAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACC                                                >  SRR3722030.323913/1‑100 (MQ=60)
                                                      AAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAA                                          <  SRR3722030.1339506/100‑1 (MQ=60)
                                                       AAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAA                                         <  SRR3722030.257024/100‑1 (MQ=60)
                                                        AAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAA                                        >  SRR3722030.3028609/1‑100 (MQ=60)
                                                          GGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATC                                      >  SRR3722030.2958459/1‑100 (MQ=60)
                                                          GGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATC                                      >  SRR3722030.1997539/1‑100 (MQ=60)
                                                          GGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATC                                      >  SRR3722030.937921/1‑100 (MQ=60)
                                                                ATGCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCtgtctcttatacacatctgacgctgccgacgacat                                >  SRR3722030.1709713/1‑65 (MQ=60)
                                                                  GCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACctgtctcttatacaca                              >  SRR3722030.3059325/1‑84 (MQ=60)
                                                                  GCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAG                              >  SRR3722030.1123730/1‑100 (MQ=60)
                                                                  GCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAG                              >  SRR3722030.1299266/1‑100 (MQ=60)
                                                                  GCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAG                              >  SRR3722030.704809/1‑100 (MQ=60)
                                                                  GCATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAG                              >  SRR3722030.1874758/1‑100 (MQ=60)
                                                                   CATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAActgtctcttatacacatctgacgctgccgacgacatggccag                             >  SRR3722030.2009801/1‑58 (MQ=60)
                                                                   CATCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGT                             >  SRR3722030.639298/1‑100 (MQ=60)
                                                                     tgtataagagacaGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGC                           <  SRR3722030.2472556/87‑1 (MQ=60)
                                                                     TCACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGC                           <  SRR3722030.379030/100‑1 (MQ=60)
                                                                      CACGCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCC                          <  SRR3722030.695176/100‑1 (MQ=60)
                                                                         GCTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGG                       <  SRR3722030.2236816/100‑1 (MQ=60)
                                                                          CTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGT                      <  SRR3722030.311588/100‑1 (MQ=60)
                                                                          CTTCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGT                      <  SRR3722030.2734806/100‑1 (MQ=60)
                                                                            TCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTAT                    <  SRR3722030.166568/100‑1 (MQ=60)
                                                                            TCGCGCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTAT                    <  SRR3722030.2411450/100‑1 (MQ=60)
                                                                              gacaGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCG                  <  SRR3722030.2573427/96‑1 (MQ=60)
                                                                                GCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAG                <  SRR3722030.2785567/100‑1 (MQ=60)
                                                                                GCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAG                <  SRR3722030.589141/100‑1 (MQ=60)
                                                                                GCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAG                <  SRR3722030.621945/100‑1 (MQ=60)
                                                                                GCGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAG                <  SRR3722030.215193/100‑1 (MQ=60)
                                                                                 acagGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGT               <  SRR3722030.2088866/96‑1 (MQ=60)
                                                                                 CGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGT               <  SRR3722030.1739268/100‑1 (MQ=60)
                                                                                 CGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCAAGTGCCAGGTATCGAGT               <  SRR3722030.2959255/100‑1 (MQ=60)
                                                                                 AGGTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGT               <  SRR3722030.1710283/100‑1 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACctgtctcttatacacatctgacgctgccgacga             >  SRR3722030.2205746/1‑67 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTctgtctctt             >  SRR3722030.2535966/1‑91 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTct             >  SRR3722030.1966996/1‑98 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTct             >  SRR3722030.934955/1‑98 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.2476088/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.956807/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.905775/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.1952377/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.1380108/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.831860/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.824513/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.78006/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.778353/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.2353779/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.2446481/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.309931/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.2090156/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.1595282/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.1643525/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.409160/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             >  SRR3722030.2862950/1‑100 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTAC             <  SRR3722030.1679449/100‑1 (MQ=60)
                                                                                   GTGTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGGGCCAGGTATCGAGTAC             >  SRR3722030.1844730/1‑100 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCtgtctcttatacacatctgacgctgccgacgacatggccagtggagatct           >  SRR3722030.150009/1‑50 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTctgt           >  SRR3722030.824087/1‑96 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTctgt           >  SRR3722030.1800274/1‑96 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTT           >  SRR3722030.3049946/1‑100 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTT           >  SRR3722030.1244712/1‑100 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTT           >  SRR3722030.932558/1‑100 (MQ=60)
                                                                                     GTCTGGGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTT           >  SRR3722030.516824/1‑100 (MQ=60)
                                                                                          GGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTTTCGCC      <  SRR3722030.2891233/100‑1 (MQ=60)
                                                                                          GGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTTTCGCC      <  SRR3722030.582467/100‑1 (MQ=60)
                                                                                          GGCGTGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTTTCGCC      <  SRR3722030.596041/100‑1 (MQ=60)
                                                                                              TGCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTTTCGCCAGCG  >  SRR3722030.571092/1‑100 (MQ=60)
                                                                                               |                                                                                                  
GATAGCGTCCAGGCAAAACCTGGAAAATTTCTTAGTTTTGCCGGATATGAAGTGAAAAGGTGAGATGCATCACGCTTCGCGCGGTGTCTGGGCGTTCAAAAAGGGCAGTAACAGCAGTGCGGAAATCCCGGCGGCGATAGAGATAACCACAAAAAATCCGCTCCAGTGCCAGGTATCGAGTACTTTCGCCAGCG  >  NZ_CP009273/3840596‑3840789

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