Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A2 F1 I1 R1 404 265.1 10478538 96.8% 10143224 150.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,704,352 (G)12→11 intergenic (‑292/+172) RPPX_25370 ← / ← RPPX_25375 GntR family transcriptional regulator/porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,704,3410G.92.9% 147.2 / 10.0 56intergenic (‑281/+183)RPPX_25370/RPPX_25375GntR family transcriptional regulator/porin
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base . (40/12);  total (41/15)
Fisher's exact test for biased strand distribution p-value = 5.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGGTGAACCCGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGCGCGC  >  CP009974/5704191‑5704491
                                                                                                                                                       |                                                                                                                                                      
gggTGAACCCGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACg                                                                                                                                                        <  2:1191541/151‑1 (MQ=255)
   tGAACCCGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGGACCCGCCCCCACgggg                                                                                                                                                     >  1:2466527/1‑151 (MQ=255)
   tGAACCCGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACgggg                                                                                                                                                     <  2:1578737/151‑1 (MQ=255)
       cccGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACgggggggg                                                                                                                                                 <  1:3699892/151‑1 (MQ=255)
        ccGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTGGGAGCGGGTTCACCCGCGAAGCAGTCAACGCGGGGCAAGGCCCCGGCGAC‑GCCGGTGTTCGCGGGCGAACCCGCTCCGACgggggcggg                                                                                                                                                >  2:1361921/1‑151 (MQ=14)
           cTCCTGCGGGGGGGCCCGAAAAAACACCAACCGGGCTCAAGACGGGGGCTGTAAGGGTGGGAGCGGGTTTCCCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACgggggggggggc                                                                                                                                             <  1:3805471/151‑2 (MQ=14)
             ccTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACgggggggggggcgg                                                                                                                                           >  1:5155270/1‑148 (MQ=255)
              cTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAT‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACgggggggggggcgcc                                                                                                                                          >  1:1022263/1‑147 (MQ=21)
                    ggggggggGCCAAGAAACCACCCGGGCACAGGACGGGCGCTGTAAGTGTAGGGGGGGGTTTACCCGCGAAGCAGGCACCGCGGGGCAGGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGa                                                                                                                                   <  1:3467792/151‑1 (MQ=11)
                    gggggCGCGCCAAGAACCCAACCTGGCTCAAGGCAGGCGCTGTAAAGGAAGGAGCGGGTTTCCCCGCGAAGCAGGCAACGCGGGGCATGGCCCCGGCTGC‑GCCGGTGTCCGCGGGTGACCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGa                                                                                                                                   <  1:3590799/151‑1 (MQ=11)
                     ggggCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGGGGCG‑CAAGaa                                                                                                                                  >  1:4434678/1‑151 (MQ=255)
                       ggCGCGGCAAGAGACCACCCTGGCTCGAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGAAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAaacc                                                                                                                               <  2:4026142/151‑1 (MQ=11)
                              cAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCCCCGGGCAC‑GCCGGTGTTCGCGGGAAAACCCGCCCCTAC‑‑‑GGGGGGGGGCGCGCCAAGAaaacaacccagc                                                                                                                       >  1:3012382/1‑147 (MQ=11)
                                    accaaccTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAAAAACAAACATGGCGCaaa                                                                                                                  >  1:4303586/1‑150 (MQ=11)
                                     ccaaccTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCAAGAAACCAACCGGgcgacag                                                                                                                  >  1:4026154/1‑147 (MQ=255)
                                      caaccTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGGGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGGCAAGAACACCACCCGGCGCAAGa                                                                                                                 >  2:2040733/1‑151 (MQ=14)
                                          cgggttCAAGACATGCGCTGGAAATGTAGGAGCGGGTTTCCCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTA‑GGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATg                                                                                                             <  1:3389715/146‑1 (MQ=14)
                                               tCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAAACAACCGGGCGCAAGACATGGGCTgg                                                                                                       >  1:3694973/1‑150 (MQ=11)
                                               gcaggACATGCGCGGAAAAGGTGGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCGGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTg                                                                                                        <  1:290525/147‑1 (MQ=14)
                                                 aaGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCACCCTGGCGCAAGACATGCGCTGTa                                                                                                      >  1:2661556/1‑151 (MQ=255)
                                                   gACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAACACAACCGGGCGCAAGACATGCGCAGTaaa                                                                                                    >  1:5163045/1‑151 (MQ=14)
                                                       ggcgcGGAAAGTGTAGGGCCGGGTTCACCCGGAGGGGGGGCAAGGCGGGGGATGGCACCGGCTCC‑GCCGGTGTTCGGGGGGGAACCCGCCCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTa                                                                                                <  1:1540321/150‑1 (MQ=9)
                                                          gcTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCTGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGGGCCAGGAAACCAACCTGCCCCAAGACCTGCGGCGCAAATGTTgggg                                                                                            >  2:671538/1‑149 (MQ=11)
                                                              tAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCCAGAAACCCACCGGGCCCACCACATGCGCCCTAAAGGTCGGCGCgg                                                                                         >  2:124583/1‑151 (MQ=14)
                                                                    ggggAGGGGGTTCGCCCGCGAGGGAGGCAAGGGGGGGAATGGCACCGGCTGG‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTAc                                                                                   <  1:1706369/149‑1 (MQ=14)
                                                                          cGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTACGAGCGGGTTTACCCGCGaa                                                                            >  1:4785418/1‑151 (MQ=255)
                                                                           gggggtACCCGCGAGGCAGGCAGGGCGGGGGAGGGCCCGGGCCAC‑GCCGGTGTTCGCGGGGGACCCCGCTCCTA‑GGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGaa                                                                            <  2:1498842/146‑1 (MQ=17)
                                                                           gggTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGGGGGGTTACCCGCGaa                                                                            >  1:4254037/1‑151 (MQ=255)
                                                                           gggTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAAACATGCGCGGTAAATGTAGAAGGGGGTTTAACCGCGaa                                                                            >  1:183940/1‑151 (MQ=17)
                                                                             gTTTACCGGTGAAGCAGGCAGCGCGGGGGAGGGCCGCGCCTAC‑GCGGGTGTGCGGGGGCGACCCCGCTCCTA‑GGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGc                                                                          <  2:4087533/151‑1 (MQ=17)
                                                                             gTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCAAGAAACCAACCTGCCCCAAGACACGCGCTGTACAGGTCGAAGCGGGTTCACCCGCGACGc                                                                          >  2:3753951/1‑151 (MQ=17)
                                                                                tACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGGAAATGTGGGAGCGGGTTTACCCGCGAAGCAgg                                                                       >  1:604810/1‑151 (MQ=255)
                                                                                    cgcgAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCCACCTGGCGCAAGACATGCGCTGTAAATGTAGGCGCGGGGTTGCCCGGGAAGCAGGccaa                                                                   >  1:2383917/1‑148 (MQ=17)
                                                                                        aaGCGGGCAAGGCGGGGCA‑GGCGCCGGGTCCGGCCGGTGTTGGGGGGTGACCCCCCCCTTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCgg                                                               <  1:982553/151‑1 (MQ=14)
                                                                                            aGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCGAGAAACCAACCTGGCGCAACACGCGCCCCGTAAATGTAGGAGCGGGTTTACCCCCCGAGCACGCCACGCcgggac                                                           >  2:842553/1‑149 (MQ=17)
                                                                                                aaCGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGGATGGc                                                       >  1:1845308/1‑151 (MQ=255)
                                                                                                aaCGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCCCAAGACACGCGCTGTAAATGTCGGAGCGGCTTTACCCCCGCAACGGGCAACGCCGGGCATgcc                                                       >  2:2433072/1‑149 (MQ=17)
                                                                                                     ggggCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCAAGAAACCAACCTGGCCCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACCCGGGGCCTGGCCCCgg                                                  >  1:942555/1‑151 (MQ=14)
                                                                                                     ggggCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGc                                                 >  1:1016040/1‑151 (MQ=255)
                                                                                                      gggCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGGGGGGCATGGCACCGGCt                                                >  2:1082603/1‑151 (MQ=255)
                                                                                                          aTGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACCTGCGCTGTAACTGTAGGAGCGGGTTGACCCCCCACGCAAGCACCCCGGGGCAGGGCCCCgggccgc                                             >  2:4123564/1‑146 (MQ=14)
                                                                                                          aTGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTCCCCCCCGAAGCAGGCAACGCGGGGCATGGCACCGGCtttcc                                            >  2:3817349/1‑149 (MQ=14)
                                                                                                            ggCACCGGCTGC‑GCCAGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGC‑CGCCAAGACACCAAGCTGGCGCAAGACATGCGTTGTAAAAGGAGGAGGGGGTTTCCCCGCGAAGGCGGCAAGGGGGGGGGCGGCGCCGGcgtcggcg                                          >  2:2144452/1‑144 (MQ=2)
                                                                                                             gCACCGGCTAC‑GCGGGTGTGCGCGGGCGACCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCg                                          <  1:4930877/151‑1 (MQ=17)
                                                                                                              cACGGGCTAC‑GCCGTGGTTAGCGGGTGAACCCGCACCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGgt                                        <  1:2144452/151‑1 (MQ=14)
                                                                                                                  ggCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCCCGAAGCAGGCAACGCGGGGGAGGGCCCCGGGTTTGCCgggggt                                     >  1:2490670/1‑147 (MQ=17)
                                                                                                                    cTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACTCGGGGCATGGCACCGGCTTTGCCGGTGTTcg                                   >  1:116305/1‑151 (MQ=255)
                                                                                                                    cTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGCAAATGTAGGAGCGGGTTTCCCCGCGAAGCCGGCAACGCGGGGCATGGCCACGGCCTTGCCggccgtct                                   >  2:3297319/1‑145 (MQ=17)
                                                                                                                     tAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCGGGCGCAAGACATGCGCTGAAAATGTAGGAGCGGGTTCACCCGCGAAACCGGCCACGCGGGGCCCGGCACCGGCTTTGCCCGGGTTcgc                                  >  2:877739/1‑151 (MQ=17)
                                                                                                                         gCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCCCTAAATGTAGGAGCGGGGTTCCCCCCGACGCACGCAACGCGGGGCTGCGCCCCGGCCTTGCGGCGGTTCCCggg                               >  2:1518187/1‑151 (MQ=14)
                                                                                                                            ggTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGGGa                            >  1:2177491/1‑151 (MQ=255)
                                                                                                                                 tCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGGTGACCCGCGAAGCAGGCCACGCGGGGGATGGCACCGGCGTTGCCGGTGGTCGCGGGGGGAcccc                       >  2:780332/1‑150 (MQ=14)
                                                                                                                                  ngcgGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCCCCCGCTTTGCCGGTGTTCGCGGGTTAACCCGc                      >  1:422363/2‑151 (MQ=255)
                                                                                                                                   gcgGGTGAACCGGCCCCTGC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCt                     <  1:137421/151‑1 (MQ=255)
                                                                                                                                    cgGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGccc                    >  1:838667/1‑149 (MQ=255)
                                                                                                                                    cgGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGGGTTCGCGGGGGAACCCGCTc                    >  1:4412421/1‑151 (MQ=255)
                                                                                                                                    cgGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCAGGGCCCCGGCTTTGCCGGGGGTCGCGGGTGAACCCGCTc                    >  1:127912/1‑151 (MQ=14)
                                                                                                                                      ggTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGTAGAGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGGGAACTCGCTCCt                  >  1:668006/1‑151 (MQ=14)
                                                                                                                                       gTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTAc                >  1:892311/1‑151 (MQ=255)
                                                                                                                                          aaCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACCCGGGGCATGGCCCCGGCTCTGCCGGGGTTCGCGGGGGAACCCGCTCCTACgg              >  2:1518582/1‑151 (MQ=14)
                                                                                                                                           gcccGCGCCTTC‑GGGGGGGGGGGCGGGGCAAGTAACCAACCGGGCGCAAGATATGCGAGGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACggg             <  1:4580561/150‑1 (MQ=14)
                                                                                                                                             ccGCTCTTAC‑GGGGGGGGGGGCGCGCAAAGAAACCAACCTGGAGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACggggg           <  1:4685173/151‑1 (MQ=255)
                                                                                                                                              cGCTCCTAC‑GGGGGGGGGGGCGCCCCAAAAAACCAACCTTGCCCAAGACATCCGCCTTAATTGTCCGAGCGGGTCTCCCCCCGAAGCCCGCACCCCCGGGCATGGCCCCGGCGTTCCCGGTCTGCCCGGGTGGACCGCCCCGTCCgggggg          >  2:534295/1‑151 (MQ=255)
                                                                                                                                               gCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGGGTTCGCGGGTGAACCCGCCCCTACggggggg         >  1:1859594/1‑151 (MQ=255)
                                                                                                                                                       ggggggcgggggCGCGCCATGACTTCTGCGTGGAGCAAGACGGGCGCGGTAGATGTAGGCGCTGGTTCACCCGCGGAGCAGGCAGCGCGGGGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGgcgcgc  <  2:4572768/151‑1 (MQ=14)
                                                                                                                                                       |                                                                                                                                                      
GGGTGAACCCGCTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTAC‑GCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGCGCGC  >  CP009974/5704191‑5704491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGCGCGCCAA  >  CP009974/5704202‑5704494
                                                                                                                                                      |                                                                                                                                              
ctcctgcgggggggcccgaaaaaacaccaaccgggctcaagacgggggctgtaagggtgGGAGCGGGTTTCCCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGC                                                                                                                                                <  SRR10918231.3823026/92‑1 (MQ=30)
   cgacgggggagagcgcacaaaagccacgcgggctgaggaacggcgctgaacgtggcggcggGGGTTTCCCCGCGAGGCGGGAAACGGGGGGCATGGCACCGGCTACCCCGGTGTTCGCGGGTGAACCCCCTCCTAC‑GGGGGGGGGGGCGCG                                                                                                                                            <  SRR10918231.3494723/90‑1 (MQ=22)
         GGGGGGGGGCCAAGAAACCACCCGGGCACAGGACGGGCGCTGTAAGTGTAGGGGGGGGTTTACCCGCGAAGCAGGCACCGCGGGGCAGGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGA                                                                                                                                      <  SRR10918231.3483576/151‑1 (MQ=19)
         GGGGGCGCGCCAAGAACCCAACCTGGCTCAAGGCAGGCGCTGTAAAGGAAGGAGCGGGTTTCCCCGCGAAGCAGGCAACGCGGGGCATGGCCCCGGCTGCGCCGGTGTCCGCGGGTGACCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGA                                                                                                                                      <  SRR10918231.3607180/151‑1 (MQ=20)
          GGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGGGGCG‑CAAGAA                                                                                                                                     >  SRR10918231.4455718/1‑151 (MQ=36)
               gcgcaaagaaaccacctggggtgaagaaagggggtgaaagTGTAGGGGCGGGTTGCCCCGCGAAGCGGGCAAGGCGGGGCAGGGCCCGGGCTGCGCCGGTGTTCGGGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAA                                                                                                                                <  SRR10918231.4753800/111‑1 (MQ=6)
                   CAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCCCCGGGCACGCCGGTGTTCGCGGGAAAACCCGCCCCTAC‑‑‑GGGGGGGGGCGCGCCAAGAAAACAACCcagc                                                                                                                          >  SRR10918231.3026175/1‑147 (MQ=7)
                     aaaaaccgcccggcgcggggatgcggcgggcaatTGTAGGGGCGGGTTTACCCGCGAAGCAGGAAAGGGGGGGCATGGCACCGGCTACGCCGGGGTTCGCGGGTGACCCCGCTCTTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCG                                                                                                                         <  SRR10918231.4707862/117‑1 (MQ=19)
                      aaacgcgaccgcgaagaaggtgggcgctgtaagggtgggaggGGGTTTACCCGCGAAGCAGGCAAGGGGGGGAATGGCACCGGCTAGGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCG                                                                                                                         <  SRR10918231.3450456/109‑1 (MQ=18)
                       aaaggacccggcgacaaggccgtgggggacaagggcaggcgcgggtacaccccggaggcaggcaaggcggggctgcccgcgggacgcgcgggggtccgggggtgaccccgctcaACGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCG                                                                                                                         <  SRR10918231.870913/37‑1 (MQ=33)
                         ACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAAAAACAAACATGGCGCAAA                                                                                                                     >  SRR10918231.4324003/1‑151 (MQ=34)
                          ccaccttggcaaaagcCATGCGCTGTAAATGAAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAG                                                                                                                     <  SRR10918231.1557022/135‑1 (MQ=37)
                          CCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCAAGAAACCAACCgggcgacag                                                                                                                     >  SRR10918231.4045063/1‑142 (MQ=37)
                              cccgcatgagGACATGCGCTGTAAATGTAGGAGCGGGTTTCCCCGCAAAGAAGGCAACGCGGGGCATGGCACCGCCTACGCCGGTGTTCGGGGGTGACCCCGCCCCGACGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACA                                                                                                                  <  SRR10918231.3933025/141‑1 (MQ=47)
                              CCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAgatcggaagagcacac                                                                                                                 >  SRR10918231.1059219/1‑135 (MQ=39)
                               cgggtTCAAGACATGCGCTGGAAATGTAGGAGCGGGTTTCCCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTA‑GGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATG                                                                                                                <  SRR10918231.3405072/146‑1 (MQ=47)
                                ggggtgaggacaggggcgGTAAATGGAGGGGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGC                                                                                                               <  SRR10918231.1392335/133‑1 (MQ=36)
                                gggcccaagaaacgcgctggaaaggtaggggcgggtttcccggggaggcaggaacgGCGGGGCAGGGCACCGGCGGCGCCGGTGTTCGCGGGGGACCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGC                                                                                                               <  SRR10918231.3499743/95‑1 (MQ=22)
                                ggaccccaggccagcgcggcagatgacGGAGCGGGTTTCCCCGCGAAGGGGGAAAGGCGGGGCATGGGACCGTCTACTCCAGGGTTCGCGGGTGAACCCGCTCAAAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGC                                                                                                               <  SRR10918231.4130453/124‑1 (MQ=7)
                                    gcaggacatgcgcggaaaaggtgGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCGGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTG                                                                                                           <  SRR10918231.291561/128‑1 (MQ=40)
                                    TCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAAACAACCGGGCGCAAGACATGGGCTGG                                                                                                          >  SRR10918231.3711907/1‑151 (MQ=30)
                                     acaagcccggtgcttcgatacgtggggggtgttcccccgggagccgggccaggcgggcgcggggccgggcgcccctgtgtttccggggtgccccccccctcgcGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTG                                                                                                           <  SRR10918231.4741258/48‑1 (MQ=48)
                                      AAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCACCCTGGCGCAAGACATGCGCTGTA                                                                                                         >  SRR10918231.2673809/1‑151 (MQ=42)
                                        GACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAACACAACCGGGCGCAAGACATGCGCAGTAAA                                                                                                       >  SRR10918231.5188446/1‑151 (MQ=35)
                                            ggcgcggaaagtgtagggccgggttcacccggagggggGGCAAGGCGGGGGATGGCACCGGCTCCGCCGGTGTTCGGGGGGGAACCCGCCCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTA                                                                                                   <  SRR10918231.1547422/113‑1 (MQ=27)
                                                 tgtttgggtgggtggcgggttgacctccgagaaggggaacgcgggatgtggaggggtggtggtgggtttttttgggcgatcccgcCCTACCGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAG                                                                                               <  SRR10918231.2955763/66‑1 (MQ=60)
                                                    aaatgtaggacagggtttaaccaggaggcagaaagagggggccacgccacaggcaaggcggttgtggggGGGTGAACCCGCTCCAAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGG                                                                                           <  SRR10918231.224283/82‑1 (MQ=29)
                                                    AAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGGCAAGAAACCAACCTGGCGCAAcacacgcccctgaaaaggtagaggggg                                                                                           >  SRR10918231.1192809/1‑124 (MQ=36)
                                                       ttggggggggggggtacccgggagggAGGCAACGCGGGGCATGGCACCGGCTAGGCGGGGGTGCGGGGGTGACCCCGCTCCTA‑GGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTT                                                                                        <  SRR10918231.5082344/125‑1 (MQ=60)
                                                        aggcagaggctggttttcccgagaggcagtcagcggggggcagggccgcggcagcgacggggttggagtgtgaccccgcacctacGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTT                                                                                        <  SRR10918231.936908/66‑1 (MQ=25)
                                                         tttttttttgttttcctGCGAAGCAGGCATCGCGGGGCTTGGCTCGGGCTAGCCCGGTTTTCGTGGGTGAGCCCTCTCCTTC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTAC                                                                                      <  SRR10918231.4526616/134‑1 (MQ=18)
                                                         ggggagggggttcgcccgcgagggaggcaagggggggaATGGCACCGGCTGGGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTAC                                                                                      <  SRR10918231.1714203/113‑1 (MQ=25)
                                                           ctccgccagtcccccccccagaaagaaaaaacggggaggccccggggccagctgggcgcgtccgctcacccccgccccctcGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACC                                                                                     <  SRR10918231.1190029/70‑1 (MQ=27)
                                                            ggtcgggtcacggggcgagagagccgccgcgggccctgcccccggctcggccggggttggggggggagccggcccctcgGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCC                                                                                    <  SRR10918231.53414/72‑1 (MQ=30)
                                                            gagctggtgcagtggctggggtgcagacgtgtgccctgccgctggcccggccggggtaggggggggagccggcccctcCGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCC                                                                                    <  SRR10918231.46340/73‑1 (MQ=32)
                                                             gagggggtttctccgcagaaccggcaaggcggggaagggcaccggctagctcggtgtttgcgggtgactgcgcccctacGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCG                                                                                   <  SRR10918231.85440/72‑1 (MQ=27)
                                                               CGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTACGAGCGGGTTTACCCGCGAA                                                                               >  SRR10918231.4808373/1‑151 (MQ=57)
                                                                GGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGGGGGGTTACCCGCGAA                                                                               >  SRR10918231.4274168/1‑151 (MQ=60)
                                                                GGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAAACATGCGCGGTAAATGTAGAAGGGGGTTTAACCGCGAA                                                                               >  SRR10918231.184604/1‑151 (MQ=50)
                                                                     TACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGGAAATGTGGGAGCGGGTTTACCCGCGAAGCAGG                                                                          >  SRR10918231.607046/1‑151 (MQ=60)
                                                                      tggcggcgcaggcggagagcggggccaggaccggggtagcggcggggttggcgggtgcccgcgatgtgaGGGGGGGGGGGGCGCGCAAAGAACCCAACCTGGCGCAAGACATGCGATGTAAATGTAGGGGCGGGTTTACCCGCGAAGCAGG                                                                          <  SRR10918231.201307/82‑1 (MQ=23)
                                                                       ccctgcgagcagagaagggtggggccggcaggcccccgggcgggtcgtggggggcgacccaaccccaacGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGC                                                                         <  SRR10918231.2590156/82‑1 (MQ=25)
                                                                       ccccgcaacgcaagcaccgtgtggtttttcacgcaccaggcggcatacgagagctgcaccagtgcctgGGGGGGGGGCGTCGCCCCAACAAACCAACCTCGCGCGAGACGCGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGC                                                                         <  SRR10918231.2590497/83‑1 (MQ=3)
                                                                         ggaggagaagacaaagcgcggcgtgacacgcgccaggcgggtggttgcgcggaacCCCGCTCCTAC‑GGAGGGGGGGGAGCGGCAAGAAACCAACCTGGCGCAAGACATGCGCGGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAAC                                                                      <  SRR10918231.1729594/96‑1 (MQ=11)
                                                                         CGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCCACCTGGCGCAAGACATGCGCTGTAAATGTAGGCGCGGGGTTGCCCGGGAAGCAGGCcaa                                                                      >  SRR10918231.2394789/1‑148 (MQ=60)
                                                                              aagcgggcaaggcggggcaggcgccgggtccggccggtgttggggggtgaccccccccttacGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGG                                                                  <  SRR10918231.986892/89‑1 (MQ=24)
                                                                                gcaggcaccggggggcgggccagagtttttgttgtgatgcaggggtcggcctacgagctcGGGGGGAGGGACGCGCCAACAAACCAACCTGGCGCAAGACTTGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGG                                                                <  SRR10918231.555300/91‑1 (MQ=10)
                                                                                aaggaaacgcagggcggggccccggctggcgggggggtgGCGGGTGGACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGC                                                               <  SRR10918231.4975443/112‑1 (MQ=27)
                                                                                 ggcgaggggggggcggggcacggcccgggcgggggtagggggggagccggaccccgcgGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGC                                                               <  SRR10918231.3419883/93‑1 (MQ=26)
                                                                                  GGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGAattgcgctctaaagggtgagcggggtgctgccgggcagtagggaacggggggaatg                                                            >  SRR10918231.3235390/1‑95 (MQ=36)
                                                                                     aaacgcgagcaagggccccggcaacgacggtgagagggggtaaacccgcccaaacGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGG                                                           <  SRR10918231.410810/96‑1 (MQ=23)
                                                                                     AACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGGATGGC                                                          >  SRR10918231.1853631/1‑151 (MQ=50)
                                                                                          GGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCCCCAAGAAACCAACCTGGCCCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACCCGGGGCCTGGCCCCGG                                                     >  SRR10918231.946626/1‑151 (MQ=31)
                                                                                          GGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGC                                                    >  SRR10918231.1020594/1‑151 (MQ=39)
                                                                                                gagtctgcaatcgacgccgcgcctcgcggcgacgccccaccctcgcgggagaaGGCGCGCCAAGAAACCAACCTGCCCCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTG                                                <  SRR10918231.4098489/98‑1 (MQ=4)
                                                                                                agggccgccgctgggcgggtgtttgtgggggagtgcgctcgcaGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTG                                                <  SRR10918231.3836443/108‑1 (MQ=32)
                                                                                                  GCACCGGCTACGCGGGTGTGCGCGGGCGACCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCG                                             <  SRR10918231.4954673/151‑1 (MQ=52)
                                                                                                   aacgggcgacggagttgtgggggggtaaaacagagtcaaaGGGGGGGGGAGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCG                                             <  SRR10918231.489382/111‑1 (MQ=38)
                                                                                                   CACGGGCTACGCCGTGGTTAGCGGGTGAACCCGCACCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGT                                           <  SRR10918231.2154066/151‑1 (MQ=32)
                                                                                                       GGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCCCGAAGCAGGCAACGCGGGGgagggccccgggtttgccgggggt                                        >  SRR10918231.2502085/1‑127 (MQ=47)
                                                                                                        GCTGCGCCGGTGTTCGGGGGTGACCCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTC                                       <  SRR10918231.4871360/151‑1 (MQ=44)
                                                                                                         CTACGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACTCGGGGCATGGCACCGGCTTTGCCGGTGTTCG                                      >  SRR10918231.116732/1‑151 (MQ=60)
                                                                                                            CGCCGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGCCGCGCCAAGAAACCAACCTGGCGCAAGACATGGGCTGTAAATGTAGGAGCGGGTTTACCCGGGAAGCAGCCAACGCGGGGCCCGGCACCGGGTTTGCCGGGGTTCGGGG                                   >  SRR10918231.2328279/1‑151 (MQ=60)
                                                                                                               CGGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTG                                >  SRR10918231.3718240/1‑151 (MQ=49)
                                                                                                                GGTGTTCGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGGGA                               >  SRR10918231.2187268/1‑151 (MQ=47)
                                                                                                                     TCGCGAGTGAACCCGCTCCTAC‑‑‑‑‑GGGGGGTCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGTGGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCC                      <  SRR10918231.4688081/151‑1 (MQ=4)
                                                                                                                      NGCGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCCCCCGCTTTGCCGGTGTTCGCGGGTTAACCCGC                         >  SRR10918231.423903/1‑151 (MQ=42)
                                                                                                                       ggcgggataccgcgcccgcaccGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGC                         <  SRR10918231.1608821/129‑1 (MQ=18)
                                                                                                                       ggcgcgtgaccccgcgcctacGGGGGGGGGGGCGCGCCAATACACCAGCCTGGCGCAAGACGTGCGCTATAACTGTACGAGCGGGTTTACCCGCGAAGCGGGCAACGCGGGGCTTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGC                         <  SRR10918231.845089/130‑1 (MQ=33)
                                                                                                                       GCGGGTGAACCGGCCCCTGC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCT                        <  SRR10918231.137923/151‑1 (MQ=36)
                                                                                                                        CGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCCC                       >  SRR10918231.842148/1‑151 (MQ=45)
                                                                                                                        CGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGGGTTCGCGGGGGAACCCGCTC                       >  SRR10918231.4433346/1‑151 (MQ=43)
                                                                                                                        CGGGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCAGGGCCCCGGCTTTGCCGGGGGTCGCGGGTGAACCCGCTC                       >  SRR10918231.128376/1‑151 (MQ=41)
                                                                                                                          gggcgaccgcccgcatacGGGGGGGGGGGCGCGCCAAGAAGCCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCC                      <  SRR10918231.348326/133‑1 (MQ=42)
                                                                                                                          GGTGAACCCGCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGTAGAGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGGGAACTCGCTCCT                     >  SRR10918231.670566/1‑151 (MQ=41)
                                                                                                                           GTGAACCCGCTCCTAC‑‑GGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTAC                   >  SRR10918231.896083/1‑151 (MQ=39)
                                                                                                                            TGAACCCGCTCCTAC‑‑‑‑‑GGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTGACCCGCGAAACAGGCAACGCGGGGCATGGCCCCGGCATTGCCGGTGTTCGgcgcaggacccgactccacgggg               >  SRR10918231.1583147/1‑128 (MQ=2)
                                                                                                                             gtgaaagagcccctccGGGGGGGGGACGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTAC                   <  SRR10918231.406970/135‑1 (MQ=39)
                                                                                                                             GAACCCGCTCCTACGGGGGGGGGGGGCGCGCCAAGAAACCAACCGGGGGCAAGAAATGCGCTGAAAATGTAGGAGCGGGTTTACCCCCGAACCAGGCAACGCGGGGCAGGGCCCCGGCTTGCCCGGGGTTCGCGGGGGAACCCGCcccccc                   >  SRR10918231.3616454/1‑145 (MQ=40)
                                                                                                                                gcccgcgccttcGGGGGGGGGGGCGGGGCAAGTAACCAACCGGGCGCAAGATATGCGAGGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGG                <  SRR10918231.4602414/139‑1 (MQ=34)
                                                                                                                                 CCGCTCTTAC‑GGGGGGGGGGGCGCGCAAAGAAACCAACCTGGAGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGG              <  SRR10918231.4707572/151‑1 (MQ=36)
                                                                                                                                   GCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGGGTTCGCGGGTGAACCCGCCCCTACGGGGGGG            >  SRR10918231.1867975/1‑151 (MQ=39)
                                                                                                                                   GCTCCTAC‑GGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAACACATGCGCTGTAAATGTAGGAGCGGGGTTGCCCGCGAAGCcggtaaagaggggcaggggcccggttcggcaggcgtttgcgggtgagccggtcttcagtgggggg            >  SRR10918231.3829848/1‑86 (MQ=22)
                                                                                                                                              GGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGCAAATGTAGGAGCGGGTTTcccagggacacaggcgacggggggcatggcaccgggtatgtcggcggttgggtggggaccctcacccgcgggggggggcgcgccaa  >  SRR10918231.4182804/1‑65 (MQ=58)
                                                                                                                                                      |                                                                                                                                              
CTCCTACGGGGGGGCGCGCCAAGAAACCAACCTGGCTCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGGGGCGCGCCAAGAAACCAACCTGGCGCAAGACATGCGCTGTAAATGTAGGAGCGGGTTTACCCGCGAAGCAGGCAACGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGGGGGGCGCGCCAA  >  CP009974/5704202‑5704494

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