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 I2 R1 884 257.2 12958910 96.9% 12557183 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,421,120 Δ1 bp noncoding (60/121 nt) rrfF ← 5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,421,1190G.96.0% 183.5 / 3.9 76noncoding (61/121 nt)rrfF5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (28/45);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 6.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207
                                                                                            |                                                                                        
cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          <  2:2967793/93‑1 (MQ=40)
cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:5356440/1‑93 (MQ=255)
cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:2051032/1‑93 (MQ=255)
cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:2703116/1‑93 (MQ=255)
 aGTCTTCCGACTGAGCCTTTTTTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  2:5203559/1‑92 (MQ=255)
 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  2:1741685/1‑92 (MQ=255)
 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         <  1:5608654/93‑2 (MQ=39)
 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         <  2:3467923/93‑2 (MQ=39)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        <  1:1426510/93‑3 (MQ=38)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        <  2:6217306/93‑3 (MQ=38)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  1:3916907/1‑91 (MQ=255)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        <  1:3028555/93‑3 (MQ=38)
   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  2:90743/1‑93 (MQ=255)
   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  1:4573655/1‑93 (MQ=255)
    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      >  2:1728091/1‑92 (MQ=255)
     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                                                    >  2:763175/1‑93 (MQ=19)
      tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                                                                   >  1:2905930/1‑93 (MQ=19)
      tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                                                                   <  1:3040515/93‑1 (MQ=16)
       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                                                  <  1:916787/93‑1 (MQ=12)
       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                                                  <  1:2406960/93‑1 (MQ=12)
       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                                                  <  1:5907208/93‑1 (MQ=12)
        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                                                                 <  1:145145/93‑1 (MQ=9)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCCGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                >  1:766777/1‑93 (MQ=9)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                >  2:3073363/1‑93 (MQ=12)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                <  2:4552273/93‑1 (MQ=9)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                <  1:4633335/93‑1 (MQ=9)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               <  2:2184662/93‑1 (MQ=9)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               >  2:2069517/1‑93 (MQ=16)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               >  2:3260236/1‑93 (MQ=16)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               >  2:2901095/1‑93 (MQ=12)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              <  1:2255576/93‑1 (MQ=9)
            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                                             >  1:943905/1‑93 (MQ=16)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            >  2:5756221/1‑93 (MQ=12)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            <  2:167601/93‑1 (MQ=9)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            >  1:717023/1‑93 (MQ=16)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          >  1:2696897/1‑93 (MQ=16)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          >  1:3832670/1‑93 (MQ=16)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          <  2:3308853/93‑1 (MQ=9)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          <  2:3925324/93‑1 (MQ=9)
                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                                                         <  1:6067680/93‑1 (MQ=9)
                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                                                         <  1:4462578/93‑1 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        <  1:6402854/93‑1 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        <  1:6402855/93‑1 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        <  1:5023873/93‑1 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        >  1:4616374/1‑93 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        >  1:3180334/1‑93 (MQ=12)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCCGAGTTCg                                                                        <  1:3073363/93‑1 (MQ=9)
                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                                                       >  1:5178838/1‑93 (MQ=12)
                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                                                      <  2:3403101/93‑1 (MQ=9)
                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                                                      >  2:3326084/1‑93 (MQ=9)
                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                                                                     <  2:3109369/93‑1 (MQ=0)
                                                           ggggAGGCCCCCCACTTGCAACGGGGCTTCGGGGTTTCCCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACa                              <  2:5707737/93‑1 (MQ=11)
                                                                   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  1:697570/93‑1 (MQ=0)
                                                                   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:288361/93‑1 (MQ=9)
                                                                   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACAACCGCGCTACTGCCGCCAGACAAATTCttt                      >  1:5644535/1‑93 (MQ=0)
                                                                    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCCCGCGCTACTGCCGCCAGACAAATTCtttt                     >  2:4214202/1‑93 (MQ=9)
                                                                    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     >  2:5749724/1‑93 (MQ=9)
                                                                    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     <  2:2649501/93‑1 (MQ=12)
                                                                    ccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     <  2:720864/93‑1 (MQ=12)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttttt                    >  1:3598535/1‑92 (MQ=0)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    >  2:3861748/1‑93 (MQ=12)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:1390203/93‑1 (MQ=12)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:3484623/93‑1 (MQ=12)
                                                                        ccTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  1:6043251/92‑1 (MQ=255)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTat                 >  1:3752861/1‑92 (MQ=16)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  1:1144225/93‑1 (MQ=255)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:6246045/93‑1 (MQ=255)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACAACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  1:2084426/93‑1 (MQ=255)
                                                                         cTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                >  2:4587246/1‑93 (MQ=255)
                                                                          tACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               >  1:5800261/1‑93 (MQ=255)
                                                                          tACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               >  2:6219271/1‑93 (MQ=255)
                                                                          tACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  2:5140924/93‑1 (MQ=255)
                                                                          tACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCAGCCAGACAAATTCTTTTCTAATc               >  2:3203981/1‑93 (MQ=255)
                                                                           aCCATCGGCGCTGCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt              <  2:3388403/93‑1 (MQ=255)
                                                                              aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           <  2:1284962/93‑1 (MQ=255)
                                                                              aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           <  2:5352052/93‑1 (MQ=255)
                                                                              aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           <  2:132769/93‑1 (MQ=255)
                                                                              aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTTTAATATGcc           >  1:438384/1‑93 (MQ=11)
                                                                               tCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg          >  1:4127249/1‑93 (MQ=255)
                                                                                cGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGa         <  1:1534163/93‑1 (MQ=255)
                                                                                  gcgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  2:6278446/93‑1 (MQ=255)
                                                                                  gcgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  1:1470590/93‑1 (MQ=255)
                                                                                  gcgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  2:4664072/93‑1 (MQ=255)
                                                                                   cgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt      <  1:2901095/93‑1 (MQ=255)
                                                                                   cgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt      <  2:5342489/93‑1 (MQ=255)
                                                                                    gcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  1:3361270/93‑1 (MQ=255)
                                                                                    gcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  1:5466217/93‑1 (MQ=255)
                                                                                    gcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  2:949714/93‑1 (MQ=255)
                                                                                      tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   >  1:393716/1‑93 (MQ=255)
                                                                                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTTCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  <  1:3868762/93‑1 (MQ=255)
                                                                                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  <  1:1584086/93‑1 (MQ=255)
                                                                                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  2:3427887/1‑93 (MQ=255)
                                                                                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  1:6058470/1‑93 (MQ=255)
                                                                                            |                                                                                        
CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAA  >  NC_000913/3421028‑3421213
                                                                                            |                                                                                             
AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTT                                                                                              >  SRR4302534.4021906/1‑94 (MQ=34)
AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTT                                                                                              <  SRR4302534.5760376/94‑1 (MQ=34)
 GTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTT                                                                                             >  SRR4302534.4697686/1‑94 (MQ=31)
 GTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTT                                                                                             <  SRR4302534.1464025/94‑1 (MQ=31)
 GTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTT                                                                                             <  SRR4302534.3107953/94‑1 (MQ=31)
   atgcgagcggatcttttagttgtagttggtgctgggcagttgtctaaTCTCGCATGAGCAGACCCCACACTACAATCGGCGCTACGGCGGTTTT                                                                                           <  SRR4302534.2371167/47‑1 (MQ=54)
    TTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAC                                                                                         >  SRR4302534.2982219/1‑94 (MQ=14)
     TCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACT                                                                                        <  SRR4302534.3120259/94‑1 (MQ=9)
      CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                       <  SRR4302534.940111/94‑1 (MQ=7)
      CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                       <  SRR4302534.6067607/94‑1 (MQ=7)
      CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                       <  SRR4302534.2470926/94‑1 (MQ=7)
       CGACTGAGCCTTTTGTTTTATTTGATGTCTGGCCGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTC                                                                                      >  SRR4302534.785835/1‑94 (MQ=4)
       CGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTC                                                                                      <  SRR4302534.148934/94‑1 (MQ=4)
        GACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCT                                                                                     <  SRR4302534.4759061/94‑1 (MQ=4)
          CTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGA                                                                                   >  SRR4302534.967921/1‑94 (MQ=4)
          CTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGA                                                                                   <  SRR4302534.2315365/94‑1 (MQ=4)
           TGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAG                                                                                  >  SRR4302534.734696/1‑94 (MQ=4)
             AGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTT                                                                                >  SRR4302534.2768148/1‑94 (MQ=4)
             AGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTT                                                                                >  SRR4302534.3935598/1‑94 (MQ=4)
               CCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCG                                                                              <  SRR4302534.6232644/94‑1 (MQ=4)
               CCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCG                                                                              >  SRR4302534.3264344/1‑94 (MQ=4)
               CCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCG                                                                              >  SRR4302534.4741644/1‑94 (MQ=4)
               CCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCG                                                                              <  SRR4302534.4583233/94‑1 (MQ=4)
                CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                             <  SRR4302534.6575717/94‑1 (MQ=4)
                CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                             <  SRR4302534.6575718/94‑1 (MQ=4)
                CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                             >  SRR4302534.5318365/1‑94 (MQ=4)
                CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                             <  SRR4302534.5159579/94‑1 (MQ=4)
                                                                  CCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTT                           <  SRR4302534.714650/94‑1 (MQ=0)
                                                                   CCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTT                          >  SRR4302534.3695310/1‑94 (MQ=0)
                                                                    CACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCT                         <  SRR4302534.1426791/94‑1 (MQ=1)
                                                                    CACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCT                         <  SRR4302534.3578295/94‑1 (MQ=1)
                                                                      CACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAT                       >  SRR4302534.3853644/1‑94 (MQ=4)
                                                                       CCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAT                      <  SRR4302534.6207531/94‑1 (MQ=9)
                                                                       ACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAT                      <  SRR4302534.1174179/94‑1 (MQ=9)
                                                                       ACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAT                      >  SRR4302534.2186477/1‑94 (MQ=9)
                                                                       ACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACAACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAT                      <  SRR4302534.2139040/94‑1 (MQ=8)
                                                                        CTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATC                     >  SRR4302534.5957662/1‑94 (MQ=11)
                                                                            CATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCG                 >  SRR4302534.4941962/1‑94 (MQ=24)
                                                                             ATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCG                >  SRR4302534.4238385/1‑94 (MQ=29)
                                                                              TCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAGATTCTTTTCTAATCTGGCGA               >  SRR4302534.5660674/1‑94 (MQ=18)
                                                                              TCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACACATTCTTTTCTAATCTGCCGA               >  SRR4302534.4842616/1‑94 (MQ=28)
                                                                               CGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAA              <  SRR4302534.1574400/94‑1 (MQ=34)
                                                                                 GCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACT            <  SRR4302534.1509275/94‑1 (MQ=39)
                                                                                  CGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTT           <  SRR4302534.2977226/94‑1 (MQ=41)
                                                                                   GCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTT          <  SRR4302534.3451303/94‑1 (MQ=44)
                                                                                   GCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTT          <  SRR4302534.5613558/94‑1 (MQ=44)
                                                                                    CTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTA         >  SRR4302534.403783/1‑94 (MQ=46)
                                                                                     TACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA        >  SRR4302534.6223208/1‑94 (MQ=52)
                                                                                     TACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA        >  SRR4302534.3848541/1‑94 (MQ=52)
                                                                                      ACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTTCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAC       <  SRR4302534.3972638/94‑1 (MQ=55)
                                                                                      ACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAC       >  SRR4302534.1970725/1‑94 (MQ=57)
                                                                                      ACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAC       <  SRR4302534.1625521/94‑1 (MQ=57)
                                                                                      ACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAC       >  SRR4302534.4782893/1‑94 (MQ=57)
                                                                                       CGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATATGCCGcatcttaacc      >  SRR4302534.734291/1‑84 (MQ=25)
                                                                                        GGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCT     >  SRR4302534.371411/1‑94 (MQ=60)
                                                                                          gcGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA    >  SRR4302534.5152364/3‑94 (MQ=60)
                                                                                           CGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAA   <  SRR4302534.91618/94‑1 (MQ=60)
                                                                                           CGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTGTTCTAtcctgccgcactttcgcctta   >  SRR4302534.4385924/1‑73 (MQ=0)
                                                                                            GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTCATCTGCCGAACTTTAACCTAAA  <  SRR4302534.4202295/94‑1 (MQ=60)
                                                                                            GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAA  >  SRR4302534.5152218/1‑94 (MQ=60)
                                                                                            |                                                                                             
AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAA  >  NC_000913/3421028‑3421213

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