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 4,064,564 A→G Y185H (TAT→CAT)  yihS ← sulfoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,064,5640AG99.1% 410.6 / ‑4.0 110Y185H (TAT→CAT) yihSsulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (42/67);  total (42/68)
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.42e-01

Only 100 of 111 total aligned reads displayed.

AATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGG  >  NZ_CP009273/4064477‑4064647
                                                                                       |                                                                                   
aaTGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATgaa                                                                                   <  1:1521959/90‑1 (MQ=255)
    attattTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAAt                                                                               >  2:165344/1‑90 (MQ=255)
       attTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCaa                                                                            >  1:2281736/1‑90 (MQ=255)
        tttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAg                                                                           >  2:1300329/1‑90 (MQ=255)
        tttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAg                                                                           >  1:2472313/1‑90 (MQ=255)
        tttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAg                                                                           >  1:1670734/1‑90 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:165344/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:304811/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:400501/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  2:2773834/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:56501/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  2:589348/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  2:691979/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:2880495/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  2:1015520/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  2:1010685/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:1051141/90‑1 (MQ=255)
         ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                          <  1:2887283/90‑1 (MQ=255)
           cTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaaa                                                                        <  2:468612/90‑1 (MQ=255)
           cTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaaa                                                                        <  2:1494669/90‑1 (MQ=255)
           cTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaaa                                                                        <  2:2333487/90‑1 (MQ=255)
            tGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAg                                                                       >  1:2608975/1‑90 (MQ=255)
            tGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAg                                                                       >  2:1974905/1‑90 (MQ=255)
             ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                                      >  1:2659718/1‑90 (MQ=255)
             ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                                      >  1:309300/1‑90 (MQ=255)
                 aCGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTcc                                                                  >  1:289088/1‑90 (MQ=255)
                 aCGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTcc                                                                  <  1:1415967/90‑1 (MQ=255)
                  cgtcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCa                                                                 >  2:2622759/1‑90 (MQ=255)
                  cgtcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCa                                                                 >  1:638377/1‑90 (MQ=255)
                  cgtcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCa                                                                 >  1:934729/1‑90 (MQ=255)
                   gtcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAc                                                                <  2:2436259/90‑1 (MQ=255)
                    tcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAcc                                                               <  2:2279468/90‑1 (MQ=255)
                     cgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACcg                                                              <  2:1967150/90‑1 (MQ=255)
                     cgtGGATAAACAAGGAAGCAAAGCGAAACGCGCAGCCCCGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAcgg                                                              <  1:518920/90‑3 (MQ=255)
                      gtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACcgc                                                             <  1:2965434/90‑1 (MQ=255)
                        ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                           <  1:495867/90‑1 (MQ=255)
                        ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                           <  1:1584157/90‑1 (MQ=255)
                        ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                           <  2:1649930/90‑1 (MQ=255)
                        ggATAATCACGGAAGCCACGCGAAGCGCGGGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                           <  1:2755131/90‑1 (MQ=255)
                         gATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTg                                                          <  1:2622759/90‑1 (MQ=255)
                          aTAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGc                                                         <  2:912043/90‑1 (MQ=255)
                          aTAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGc                                                         <  1:1185176/90‑1 (MQ=255)
                            aaTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCat                                                       >  1:2973810/1‑90 (MQ=255)
                            aaTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCat                                                       >  2:2592578/1‑90 (MQ=255)
                            aaTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCat                                                       >  1:1766964/1‑90 (MQ=255)
                            aaTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCat                                                       >  1:1572881/1‑90 (MQ=255)
                             aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                      >  2:1578854/1‑90 (MQ=255)
                             aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                      >  2:2226399/1‑90 (MQ=255)
                             aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                      >  2:2459116/1‑90 (MQ=255)
                             aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                      >  2:2474166/1‑90 (MQ=255)
                             aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                      >  2:1196884/1‑90 (MQ=255)
                             aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                      >  1:947110/1‑90 (MQ=255)
                              tCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCatat                                                     >  1:2542786/1‑90 (MQ=255)
                              tCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCatat                                                     <  2:706107/90‑1 (MQ=255)
                              tCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCatat                                                     <  1:1111064/90‑1 (MQ=255)
                                aCGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTg                                                   >  2:427926/1‑90 (MQ=255)
                                aCGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTg                                                   >  1:1012741/1‑90 (MQ=255)
                                 cGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTgg                                                  >  1:2181672/1‑90 (MQ=255)
                                  ggAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGc                                                 >  2:2812338/1‑90 (MQ=255)
                                  ggAAGCCACGCGAATCCCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGc                                                 >  1:2248833/1‑90 (MQ=255)
                                      gCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTg                                             <  2:115494/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  1:100799/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  2:1132941/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  1:435501/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  2:2472313/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  1:298004/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  2:1529190/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  1:2358973/90‑1 (MQ=255)
                                        cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  1:2473413/90‑1 (MQ=255)
                                               aTCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCgg                                    >  1:72054/1‑90 (MQ=255)
                                               aTCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCgg                                    <  2:246329/90‑1 (MQ=255)
                                               aTCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCgg                                    <  1:2155507/90‑1 (MQ=255)
                                                  gcgcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAc                                 <  1:570383/90‑1 (MQ=255)
                                                  gcgcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAc                                 <  1:2567078/90‑1 (MQ=255)
                                                  gcgcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAc                                 <  2:1927098/90‑1 (MQ=255)
                                                  gcgcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAc                                 <  2:272584/90‑1 (MQ=255)
                                                    gcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGc                                                         >  1:1813083/1‑64 (MQ=255)
                                                    gcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGc                                                         <  2:1813083/64‑1 (MQ=255)
                                                       aTCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTc                            <  2:289088/90‑1 (MQ=255)
                                                       aTCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTc                            <  1:1931179/90‑1 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           >  1:1329107/1‑65 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                           <  2:1329107/65‑1 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTc                            <  1:2192968/80‑1 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTc                            >  2:2192968/1‑80 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg                  >  1:2127951/1‑90 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg                  >  1:2883562/1‑90 (MQ=255)
                                                                 tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg                  >  2:530985/1‑90 (MQ=255)
                                                                      gTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGc             >  1:2147361/1‑90 (MQ=255)
                                                                      gTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGc             >  2:339917/1‑90 (MQ=255)
                                                                        cGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCtt           <  1:1196884/90‑1 (MQ=255)
                                                                         gtgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTc          <  1:204264/90‑1 (MQ=255)
                                                                         gtgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTc          <  1:2790974/90‑1 (MQ=255)
                                                                         gtgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTc          >  1:145691/1‑90 (MQ=255)
                                                                              tgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTccc     <  1:79123/90‑1 (MQ=255)
                                                                                aCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATCTTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAg   <  2:2281736/90‑1 (MQ=255)
                                                                                aCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAg   <  1:1888617/90‑1 (MQ=255)
                                                                                aCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAg   <  1:1366480/90‑1 (MQ=255)
                                                                                aCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAg   <  2:3006789/90‑1 (MQ=255)
                                                                                 cGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAgg  <  2:639688/90‑1 (MQ=255)
                                                                                    cATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAgg  <  2:921641/87‑1 (MQ=255)
                                                                                       |                                                                                   
AATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGG  >  NZ_CP009273/4064477‑4064647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ 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
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  NZ_CP009273/4064468‑4064658
                                                                                                |                                                                                              
ataagagacagTGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.12191/89‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.1000583/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.2203875/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.2473869/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.2069396/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.1965301/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.1459460/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.1253033/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.690813/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.831708/100‑1 (MQ=60)
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAAC                                                                                             <  SRR3722030.1249748/100‑1 (MQ=60)
      GATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAA                                                                                       >  SRR3722030.2323054/1‑100 (MQ=60)
       ATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAG                                                                                      >  SRR3722030.2517479/1‑100 (MQ=60)
       ATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAG                                                                                      >  SRR3722030.1697866/1‑100 (MQ=60)
         AATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAA                                                                                    <  SRR3722030.1546498/100‑1 (MQ=60)
           TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAG                                                                                  >  SRR3722030.2656882/1‑100 (MQ=60)
            GATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGC                                                                                 >  SRR3722030.313804/1‑100 (MQ=60)
            GATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGC                                                                                 >  SRR3722030.2708710/1‑100 (MQ=60)
                ATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCC                                                                             >  SRR3722030.293276/1‑100 (MQ=60)
                 TTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCA                                                                            >  SRR3722030.648284/1‑100 (MQ=60)
                 TTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCA                                                                            >  SRR3722030.949028/1‑100 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.167720/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.2488838/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.1066939/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.2933649/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.2940556/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.309254/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.404052/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.406526/100‑1 (MQ=60)
                  TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722030.57372/100‑1 (MQ=60)
                          ACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCA                                                                   <  SRR3722030.1438316/100‑1 (MQ=60)
                           CGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCAT                                                                  >  SRR3722030.1796385/1‑100 (MQ=60)
                           CGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCAT                                                                  >  SRR3722030.3028603/1‑100 (MQ=60)
                           CGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCAT                                                                  >  SRR3722030.1598342/1‑100 (MQ=60)
                            GTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATA                                                                 >  SRR3722030.961577/1‑100 (MQ=60)
                             TCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATAT                                                                >  SRR3722030.2589363/1‑100 (MQ=60)
                             TCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATAT                                                                >  SRR3722030.841488/1‑100 (MQ=60)
                              cgtggataaacaaggaagcaaagcgaaacgcgcagccccGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACGGCGTGCATATT                                                               <  SRR3722030.526890/61‑1 (MQ=60)
                               GTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTG                                                              <  SRR3722030.3020098/100‑1 (MQ=60)
                               GTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTG                                                              <  SRR3722030.2020557/100‑1 (MQ=60)
                               GTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTG                                                              >  SRR3722030.1028062/1‑100 (MQ=60)
                                TGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGG                                                             >  SRR3722030.2220854/1‑100 (MQ=60)
                                 GGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                            <  SRR3722030.503468/100‑1 (MQ=60)
                                 GGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                            <  SRR3722030.1609833/100‑1 (MQ=60)
                                 GGATAATCACGGAAGCCACGCGAATCCCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                            >  SRR3722030.2289431/1‑100 (MQ=60)
                                 GGATAATCACGGAAGCCACGCGAAGCGCGGGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                            <  SRR3722030.2806017/100‑1 (MQ=60)
                                  GATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCA                                                           <  SRR3722030.2671009/100‑1 (MQ=60)
                                   ATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCAT                                                          <  SRR3722030.1202736/100‑1 (MQ=60)
                                       TCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCC                                                      <  SRR3722030.1127649/100‑1 (MQ=60)
                                              AGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGG                                               >  SRR3722030.73157/1‑100 (MQ=60)
                                                 tgtcgtcggcagcgtcagatgtgtataagagacagGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722030.2106747/65‑1 (MQ=60)
                                                 CACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722030.442097/100‑1 (MQ=60)
                                                 CACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722030.302350/100‑1 (MQ=60)
                                                 CACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722030.102288/100‑1 (MQ=60)
                                                 CACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722030.2401862/100‑1 (MQ=60)
                                                 CACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722030.2518600/100‑1 (MQ=60)
                                                   CGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATctgtctcttatacaca                                          >  SRR3722030.1843722/1‑84 (MQ=60)
                                                        ATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGG                                     <  SRR3722030.2194123/100‑1 (MQ=60)
                                                           GCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTT                                  >  SRR3722030.1543758/1‑100 (MQ=60)
                                                           GCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTT                                  <  SRR3722030.579218/100‑1 (MQ=60)
                                                           GCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTT                                  <  SRR3722030.2614157/100‑1 (MQ=60)
                                                                ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACctgtctcttatacac                             >  SRR3722030.1349545/1‑85 (MQ=60)
                                                                ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             >  SRR3722030.2936770/1‑100 (MQ=60)
                                                                ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             <  SRR3722030.2864238/100‑1 (MQ=60)
                                                                ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             <  SRR3722030.2232380/100‑1 (MQ=60)
                                                                ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             >  SRR3722030.2165959/1‑100 (MQ=60)
                                                                ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             <  SRR3722030.1964774/100‑1 (MQ=60)
                                                                     GCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGC                        >  SRR3722030.2185772/1‑100 (MQ=60)
                                                                        ATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTC                     >  SRR3722030.147787/1‑100 (MQ=60)
                                                                                 CGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGA            <  SRR3722030.1214601/100‑1 (MQ=60)
                                                                                  GTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGAT           <  SRR3722030.207172/100‑1 (MQ=60)
                                                                                  GTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGAT           <  SRR3722030.2842512/100‑1 (MQ=60)
                                                                                   TGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATT          >  SRR3722030.1303885/1‑100 (MQ=60)
                                                                                   TGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATT          >  SRR3722030.935759/1‑100 (MQ=60)
                                                                                       TGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAG      <  SRR3722030.80346/100‑1 (MQ=60)
                                                                                         ACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGC    <  SRR3722030.1921140/100‑1 (MQ=60)
                                                                                         ACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGC    <  SRR3722030.1387816/100‑1 (MQ=60)
                                                                                           GTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  SRR3722030.1628790/1‑100 (MQ=60)
                                                                                                |                                                                                              
TAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  NZ_CP009273/4064468‑4064658

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