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 F1 I1 R1 343 268.7 10873516 97.0% 10547310 150.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,837,678 C→T intergenic (+229/‑322) RPPX_17095 → / → RPPX_17100 proteophosphoglycan precursor/cupin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,837,6780CT85.4% 150.8 / 21.0 89intergenic (+229/‑322)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base C (9/4);  new base T (31/45);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 7.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTA  >  CP009974/3837529‑3837816
                                                                                                                                                     |                                                                                                                                          
ttCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCACTTTTCCCCCGCCGCTCACCGTTGCTGCCGCCAccc                                                                                                                                        >  2:2202478/1‑150 (MQ=2)
 tCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAc                                                                                                                                          >  1:2443675/1‑151 (MQ=255)
 ccACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAc                                                                                                                                          <  2:709736/150‑1 (MQ=255)
  cACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAcc                                                                                                                                         <  1:992327/151‑1 (MQ=255)
  cACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCAGCTGGCCAATTTCCACCCGCCCTTCACCCTTGCCGCCGCTAcc                                                                                                                                         >  2:633630/1‑151 (MQ=255)
   aCCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGGTCACCCGCGCAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAccc                                                                                                                                        >  1:217464/1‑150 (MQ=255)
    ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGt                                                                                                                                       <  1:70630/151‑1 (MQ=255)
    ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTCGCCAATTACCACCCGCCGTTCACCGTTGCCGCCGCTACCGt                                                                                                                                       >  1:802435/1‑151 (MQ=255)
    ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTGGCAGGAGCGGGGTCACCCGCGAGAGAGCCGGGCCAGCTGGCCAACTTCCCCCCGCCCTTCACCGTCGCGGCCGCCACTGt                                                                                                                                       >  2:2229877/1‑151 (MQ=14)
     cATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTCCCGtt                                                                                                                                      >  1:2421116/1‑150 (MQ=255)
     cATTGCCGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTa                                                                                                                                      <  1:1134483/151‑1 (MQ=255)
      aTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGCGGGCGGCACAGGCCGCTGTAGGCGCGGGGTCACCCGCGCAAGGGCCGGGCCTGCTGGCCAATTTCCAACCGCCGTTCCCCGTTGCCGCCGCTACCGTAt                                                                                                                                     >  2:4766797/1‑151 (MQ=11)
       ttGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAtt                                                                                                                                    <  1:2013621/151‑1 (MQ=38)
       ttGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAtt                                                                                                                                    <  1:2013048/151‑1 (MQ=38)
        gggcgccgcTACCGTAGGCGCGGGTAGCCCCGCTCCGGCACACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTc                                                                                                                                   <  2:2407442/148‑1 (MQ=11)
         gccgccGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTcg                                                                                                                                  <  1:799036/151‑1 (MQ=38)
           cgccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTcgcg                                                                                                                                >  1:2949723/1‑151 (MQ=37)
            gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                                               <  1:5104766/151‑1 (MQ=38)
            gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                                               >  2:4294160/1‑151 (MQ=38)
              cgcgagcGTATTCGCGGGTAAACCCGCTCCTACACACGCAGCGAGGGCGGCAGAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:5269376/145‑1 (MQ=11)
              cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:4069824/151‑1 (MQ=37)
              cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:573860/151‑1 (MQ=38)
              cgcTACCGTATGCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:4071443/151‑1 (MQ=255)
               gctagggTATTCGCGGGTAAACCCGCTCCTACGCACGCGGCGAGGGCGGCACAGGCCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTa                                                                                                                            <  2:2284682/145‑1 (MQ=11)
                cTACGGGATTCGCTGGTAAACCCGCTGCTACACGCGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaa                                                                                                                           <  2:1353523/151‑1 (MQ=11)
                cTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaa                                                                                                                           <  1:3941279/151‑1 (MQ=37)
                 tACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaaa                                                                                                                          >  1:2836770/1‑151 (MQ=21)
                 tACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCCGTGCTGCCCCTTTCTTATTCGCGGGTaag                                                                                                                          >  2:4340562/1‑150 (MQ=11)
                  agggTAGTCGCGGGTAGACCGGGTCCAACACCGGCGTCGAGGGCGGCACAGGCCGTCGTATGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAc                                                                                                                         <  2:3287090/148‑1 (MQ=11)
                  aCCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAc                                                                                                                         >  2:943127/1‑151 (MQ=25)
                   gcggaTTCGCGGGTAAACCCGCCGCTACACACGCATCGAGGGCGGCCCAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAcc                                                                                                                        <  2:397137/147‑1 (MQ=11)
                   ccGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAcc                                                                                                                        >  2:3492075/1‑151 (MQ=18)
                     gTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                                                                                      >  1:1120971/1‑151 (MQ=25)
                      tATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                                                                                     <  1:785044/151‑1 (MQ=25)
                      tATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                                                                                     <  1:786730/151‑1 (MQ=25)
                       aTTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGCAAACCCGcc                                                                                                                    >  2:2885058/1‑150 (MQ=25)
                        ttCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTc                                                                                                                   >  1:132710/1‑151 (MQ=21)
                         tCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                                                  >  1:899208/1‑151 (MQ=21)
                          tgcgGGTGAACCGGGTCCTACACATGCAGCGAGGGCGGCACAGGCCGTTGTAGGGGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 <  2:2703079/150‑1 (MQ=11)
                          cgcgGGTAAGACCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 <  1:1055907/151‑1 (MQ=255)
                          cgcgGGTAAACCCGCTCCTACCCACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 <  1:2410448/151‑1 (MQ=255)
                          cgcgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 >  1:3332573/1‑151 (MQ=25)
                           gcgGGTAAACCCGCTCCTGCACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtg                                                                                                                <  2:1552571/151‑2 (MQ=19)
                            cgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTACACCCGCTCCtgc                                                                                                               >  2:2985084/1‑149 (MQ=12)
                            cgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtgc                                                                                                               >  1:493947/1‑149 (MQ=18)
                            cgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtgc                                                                                                               >  1:5261179/1‑149 (MQ=18)
                             gggTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgca                                                                                                              <  1:1680035/151‑1 (MQ=18)
                              ggTAACCCGGCTCCGACACACGCATCGAGGGTGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcac                                                                                                             <  2:807970/151‑1 (MQ=11)
                              ggTAAACCCGCTCCTACACACGCATCGGGGGCGGCACAGGCCGTTTGAGGAGCGGGTTCACCCGCGAAAGGGCCGTGCCTGCTGGCCACTTTCTACCCGCAGTCCACCGTTGCGCCGGCTACTCTATTCGCTGGTACAACGGCTCCTtcac                                                                                                             >  2:3456076/1‑151 (MQ=11)
                              gctagaCCCGCCCCTACGCAGGCATCGAGGGCTGCAGAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcac                                                                                                             <  2:4290510/146‑1 (MQ=11)
                               gTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcccgcacc                                                                                                            >  1:1958521/1‑145 (MQ=16)
                                tAACCCCGCTCCTACACACGCCTCGTGGGCGGCGCTGCCCCTTGTACGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcacac                                                                                                           <  2:4356744/151‑1 (MQ=11)
                                 aaaCCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACg                                                                                                          >  2:3752476/1‑151 (MQ=18)
                                 aaaCCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACg                                                                                                          >  2:731777/1‑151 (MQ=18)
                                     ccGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCCCGCACACGCCTc                                                                                                      >  2:2268122/1‑151 (MQ=255)
                                       gCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCgg                                                                                                    >  2:5288804/1‑150 (MQ=2)
                                        cTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCggt                                                                                                   <  1:5078801/151‑3 (MQ=9)
                                        cTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCggt                                                                                                   >  1:162919/1‑149 (MQ=9)
                                                                                        cagccGCGAAAGGGGCGGGCCTGCTGGTCAATTTTCACCGGCCGATCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGGGTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCct                                                <  1:4181969/148‑1 (MQ=2)
                                                                                         cccagcgATAGGGCCGGGTCTGCTGGCCAATTCTCACCTGCCGTTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACTC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                                               <  2:1787585/147‑1 (MQ=9)
                                                                                         aCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                                               <  1:4284316/151‑1 (MQ=2)
                                                                                         aCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                                               <  1:5286651/151‑1 (MQ=2)
                                                                                           ccGCGAAGGGGCCGGTCCTGCTGGCCAATTTTCACCAGCCGTTCTCGGTTGCTGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  2:2345951/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCGTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCGGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  1:3502551/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCGAATTTTCACCAGCCGTTCAACGTTGCTGCCGCTACCGTATTCGCGGGTAACCCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  1:4434930/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTCCCGTATTCGCGGGTAACCCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  2:31401/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             >  2:235082/1‑151 (MQ=11)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  1:170097/151‑1 (MQ=2)
                                                                                            cgcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCGGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                                            <  2:150742/151‑1 (MQ=11)
                                                                                            cgcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                                            <  2:3227558/151‑1 (MQ=2)
                                                                                             ccgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTgg                                           <  2:2073831/150‑1 (MQ=2)
                                                                                              cgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGACAGGGCCGGGCCTGCTGGc                                          >  2:4041131/1‑151 (MQ=2)
                                                                                               gAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCCACCGAATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGcc                                         <  2:5069675/151‑1 (MQ=9)
                                                                                                aaaGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCa                                        <  1:3455756/151‑1 (MQ=2)
                                                                                                 aaGGGCCGGGCCTGCTGGCCAGTTTTCACCAGCCGTTCACCGTTGCTGCCGCTCCCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCaa                                       <  2:1451941/151‑1 (MQ=2)
                                                                                                 aaGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCaa                                       <  2:3097971/151‑1 (MQ=2)
                                                                                                  aGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                      >  1:5408570/1‑151 (MQ=2)
                                                                                                  aGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                      <  1:618014/151‑1 (MQ=2)
                                                                                                  aGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                      <  1:39988/151‑1 (MQ=2)
                                                                                                       cgggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAttttca                                 <  1:4166300/151‑6 (MQ=2)
                                                                                                       cgggccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAttttcc                                 >  1:3461420/1‑146 (MQ=2)
                                                                                                         ggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTACACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAcc                               <  1:615587/151‑1 (MQ=2)
                                                                                                           ccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccag                             >  1:3869106/1‑149 (MQ=2)
                                                                                                            ctgctgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagc                            >  2:613653/1‑148 (MQ=2)
                                                                                                              gctgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccg                          >  1:2480953/1‑146 (MQ=2)
                                                                                                                tgGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgcc                        >  2:2099875/1‑144 (MQ=2)
                                                                                                                tgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgct                        >  2:862381/1‑144 (MQ=2)
                                                                                                                tgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgcc                        >  2:567042/1‑144 (MQ=2)
                                                                                                                                      cACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACAC‑TG‑ATC‑TGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCCGGCCCATCTTCACCACCCGCTCACCCCTGCCCCCGCCACCGTa  >  2:1106872/1‑151 (MQ=2)
                                                                                                                                                     |                                                                                                                                          
TTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTA  >  CP009974/3837529‑3837816

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTC  >  CP009974/3837529‑3837835
                                                                                                                                                     |                                                                                                                                                             
TTCACCATTGCCGCCGCTACCGTATGCGCGGGTACACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTA                                                                                                                                                              <  SRR10918230.3273024/151‑1 (MQ=46)
 TCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAC                                                                                                                                                             >  SRR10918230.4543343/1‑151 (MQ=47)
 TCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAC                                                                                                                                                             >  SRR10918230.2454162/1‑151 (MQ=47)
  CACCATGGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACC                                                                                                                                                            <  SRR10918230.996365/151‑1 (MQ=37)
    CCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGT                                                                                                                                                          <  SRR10918230.70858/151‑1 (MQ=38)
    CCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTCGCCAATTACCACCCGCCGTTCACCGTTGCCGCCGCTACCGT                                                                                                                                                          >  SRR10918230.805510/1‑151 (MQ=36)
     CATTGCCGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTA                                                                                                                                                         <  SRR10918230.1139333/151‑1 (MQ=35)
       TTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATT                                                                                                                                                       <  SRR10918230.2021549/151‑1 (MQ=31)
       TTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATT                                                                                                                                                       <  SRR10918230.2022123/151‑1 (MQ=31)
       TTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGGTCCCCGGCGGAAGGGCCGTGCCTGCTGGCCAATTTCCACCCGCCCTTCACCGTGGCCGCAGCTACCGTATT                                                                                                                                                       >  SRR10918230.5283687/1‑151 (MQ=25)
         GCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCG                                                                                                                                                     <  SRR10918230.802099/151‑1 (MQ=31)
           CGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCG                                                                                                                                                   >  SRR10918230.2962547/1‑151 (MQ=31)
            GCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGG                                                                                                                                                  <  SRR10918230.5128189/151‑1 (MQ=31)
              cgctacggggtgtgctggtagACCCGCCCCTACACCTGCATCGTGGGCGGGCCAGGCCGTCGTATGAGCGGGTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                                                                                                                                                <  SRR10918230.4941086/130‑1 (MQ=0)
              CGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                                                                                                                                                <  SRR10918230.4087963/151‑1 (MQ=29)
              CGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                                                                                                                                                <  SRR10918230.575904/151‑1 (MQ=29)
              CGCTACCGTATGCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                                                                                                                                                <  SRR10918230.4089591/151‑1 (MQ=29)
              CGCGAGCGTATTCGCGGGTAAACCCGCTCCTACACACGCAGCGAGGGCGGCAGAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                                                                                                                                                <  SRR10918230.5293812/151‑1 (MQ=19)
                CTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAA                                                                                                                                              <  SRR10918230.3958572/151‑1 (MQ=25)
                 TACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAA                                                                                                                                             >  SRR10918230.2849152/1‑151 (MQ=23)
                   CCGTATTCGCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGCAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGAGTAAACCCGCTCCTACA                                                                                                                                 >  SRR10918230.1820179/1‑151 (MQ=0)
                     GTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCG                                                                                                                                         >  SRR10918230.1125753/1‑151 (MQ=21)
                     GTATTCGCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACA                                                                                                                               >  SRR10918230.3420515/1‑151 (MQ=0)
                      TATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGC                                                                                                                                        <  SRR10918230.788041/151‑1 (MQ=21)
                      TATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGC                                                                                                                                        <  SRR10918230.789733/151‑1 (MQ=21)
                        TTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTC                                                                                                                                      >  SRR10918230.133170/1‑151 (MQ=21)
                         TCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCC                                                                                                                                     >  SRR10918230.902765/1‑151 (MQ=21)
                          CGCGGGTAAGACCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCT                                                                                                                                    <  SRR10918230.1060332/151‑1 (MQ=19)
                          CGCGGGTAAACCCGCTCCTACCCACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCT                                                                                                                                    <  SRR10918230.2420778/151‑1 (MQ=20)
                          CGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCGGTTGTAGTAGCGGGTTCCCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCT                                                                                                                                    <  SRR10918230.3838877/151‑1 (MQ=19)
                          CGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCT                                                                                                                                    >  SRR10918230.3347033/1‑151 (MQ=21)
                           GCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACatgca                                                                                                                         >  SRR10918230.1010947/1‑146 (MQ=0)
                            CGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGC                                                                                                                                  >  SRR10918230.495664/1‑151 (MQ=17)
                            CGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGC                                                                                                                                  >  SRR10918230.5285576/1‑151 (MQ=17)
                             GGGTAACCCCGCTCTTCCACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCA                                                                                                                                 <  SRR10918230.1287265/151‑1 (MQ=10)
                             GGGTAACCCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCA                                                                                                                                 <  SRR10918230.1687388/151‑1 (MQ=11)
                               GTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCcgcacc                                                                                                                               >  SRR10918230.1966808/1‑145 (MQ=21)
                                        CTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgt                                                                                                                      <  SRR10918230.5102089/151‑3 (MQ=3)
                                        CTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgt                                                                                                                      >  SRR10918230.163475/1‑149 (MQ=3)
                                           CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCctcggtcgt                                                                                                                   >  SRR10918230.4495702/1‑142 (MQ=5)
                                           CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtggt                                                                                                                   <  SRR10918230.1935105/151‑6 (MQ=0)
                                                                                                                                                  CGCCACCGTATTCGCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTC  >  SRR10918230.2847175/1‑151 (MQ=0)
                                                                                                                                                     |                                                                                                                                                             
TTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTC  >  CP009974/3837529‑3837835

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