Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F3 I1 R1 15 28.0 1143766 97.4% 1114028 135.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,964,929 C→T E234E (GAG→GAA rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,964,9290CT100.0% 237.8 / NA 68E234E (GAG→GAArhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/31);  total (37/31)

CCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCG  >  NC_000913/3964792‑3965071
                                                                                                                                         |                                                                                                                                              
ccATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAAt                                                                                                                                     <  5:103947/149‑1 (MQ=255)
 cATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATa                                                                                                                                    <  5:121833/149‑1 (MQ=255)
   tCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACg                                                                                                                                  >  5:56803/1‑149 (MQ=255)
       ctgcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTgg                                                                                                                              >  5:32392/1‑149 (MQ=255)
         gcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGcc                                                                                                                            >  5:106341/1‑149 (MQ=255)
              gTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGttt                                                                                                                       <  6:32392/149‑1 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  6:72497/1‑149 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  3:28984/1‑149 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  6:75800/1‑149 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  3:74206/1‑149 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  5:29523/1‑149 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  1:130408/1‑149 (MQ=255)
                   cccAGATCTCTTCACAACGGTGTTAGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                  >  6:128029/1‑149 (MQ=255)
                    ccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGtt                                                                                                                 >  3:38133/1‑149 (MQ=255)
                             ttCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCgg                                                                                                             <  6:29523/144‑1 (MQ=255)
                              tCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAc                                                                                                       >  7:66605/1‑149 (MQ=255)
                               cacaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACt                                                                                                      >  5:76964/1‑149 (MQ=255)
                                acaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACtt                                                                                                     <  5:72497/149‑1 (MQ=255)
                                  aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTc                                                                                                    >  3:32591/1‑148 (MQ=255)
                                    cGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                 <  4:74206/149‑1 (MQ=255)
                                                ttGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGc                                                                                        >  7:25598/1‑146 (MQ=255)
                                                ttGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGc                                                                                        <  8:25598/146‑1 (MQ=255)
                                                   gCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGt                                                                                                                         <  3:39811/110‑1 (MQ=255)
                                                   gCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGt                                                                                                                         >  4:39811/1‑110 (MQ=255)
                                                      aaaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtctt                                                                                                                                            <  7:95491/88‑1 (MQ=255)
                                                      aaaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtctt                                                                                                                                            >  8:95491/1‑88 (MQ=255)
                                                       aaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATc                                                                              <  6:56803/149‑1 (MQ=255)
                                                          taatCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGc                                                                           <  4:32591/149‑1 (MQ=255)
                                                           aatCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCt                                                                          >  8:64631/1‑149 (MQ=255)
                                                            atCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTc                                                                         <  4:28984/149‑1 (MQ=255)
                                                                        ggCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGtt                                                              >  5:19926/1‑148 (MQ=255)
                                                                           cACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACg                                                          >  5:115969/1‑149 (MQ=255)
                                                                            aCTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAt                                                                                      <  8:97504/120‑1 (MQ=255)
                                                                                ttcttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCTGCATTGTTCATCTGCTCGAACTCCAGTTCACGTAccc                                                     >  8:26235/1‑149 (MQ=255)
                                                                                ttcttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAccc                                                     >  3:101742/1‑149 (MQ=255)
                                                                                ttcttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAc                                                       <  5:75800/147‑1 (MQ=255)
                                                                                  cttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGTGTAGAAAAGTTATTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCgg                                                   >  3:92406/1‑149 (MQ=255)
                                                                                  cttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAccc                                                     <  7:64631/147‑1 (MQ=255)
                                                                                   ttcGATCAGCGGTTGCAGCAAACGCATTTTTTCTTAGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGTTTCCACTTAAATATATTCGGCATTGATCATCTGCTCGAACGCCAGTTCACGTACCCgg                                                   >  3:41755/1‑148 (MQ=255)
                                                                                          aGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGaaa                                            <  4:38133/148‑1 (MQ=255)
                                                                                          aGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGaaa                                            <  6:76964/148‑1 (MQ=255)
                                                                                           gCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTgg                                                                                                                              <  1:60241/65‑1 (MQ=255)
                                                                                           gCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTgg                                                                                                                              >  2:60241/1‑65 (MQ=255)
                                                                                                tGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTgg                                      <  6:106341/148‑1 (MQ=255)
                                                                                                 gcagcaAACGCATTTTTTCTTCGTTAGTAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGGTCGACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGc                                     <  5:128029/148‑1 (MQ=255)
                                                                                                 gcagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGc                                     <  4:101742/148‑1 (MQ=255)
                                                                                                   agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCg                                    >  2:133605/1‑147 (MQ=255)
                                                                                                       aaCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTACACGTACCCGGTACGAAAGCGTGGCGGAGaa                               >  1:138432/1‑148 (MQ=255)
                                                                                                        aCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGaa                               >  1:137765/1‑147 (MQ=255)
                                                                                                         cGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAg                            <  6:19926/149‑1 (MQ=255)
                                                                                                           cATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCa                          <  6:115969/149‑1 (MQ=255)
                                                                                                                 ttcttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGtt                      >  6:102046/1‑147 (MQ=255)
                                                                                                                    ttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGtt                      <  1:100621/144‑1 (MQ=255)
                                                                                                                    ttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGtt                      >  2:100621/1‑144 (MQ=255)
                                                                                                                    ttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGgc                 >  2:97810/1‑149 (MQ=255)
                                                                                                                      cGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACt                                                                                                      >  3:18911/1‑62 (MQ=255)
                                                                                                                      cGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACt                                                                                                      <  4:18911/62‑1 (MQ=255)
                                                                                                                      cGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGgc                 <  5:102046/147‑1 (MQ=255)
                                                                                                                        ttAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                    >  7:102753/1‑142 (MQ=255)
                                                                                                                        ttAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                    <  8:102753/142‑1 (MQ=255)
                                                                                                                        ttAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTg             <  8:66605/149‑1 (MQ=255)
                                                                                                                              gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTg        >  5:91493/1‑148 (MQ=255)
                                                                                                                              gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTg        >  3:47447/1‑148 (MQ=255)
                                                                                                                              gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTg        >  3:49550/1‑148 (MQ=255)
                                                                                                                              gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGc       <  1:97810/149‑1 (MQ=255)
                                                                                                                                  aGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCa      >  2:102440/1‑146 (MQ=255)
                                                                                                                                  aGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCa      <  1:102440/146‑1 (MQ=255)
                                                                                                                                   gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcg  <  1:133605/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCG  >  NC_000913/3964792‑3965071

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG  >  NC_000913/3964791‑3965072
                                                                                                                                          |                                                                                                                                               
ACCATATGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGGGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATA                                                                                                                                     <  NB501016:228:HFK2FAFXY:3:21502:14477:18532/151‑1 (MQ=60)
 CCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATAC                                                                                                                                    <  NB501016:228:HFK2FAFXY:3:21512:15738:3921/151‑1 (MQ=60)
   ATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGG                                                                                                                                  >  NB501016:228:HFK2FAFXY:3:11603:15913:17050/1‑151 (MQ=60)
       GCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGC                                                                                                                              >  NB501016:228:HFK2FAFXY:3:11504:14513:10652/1‑151 (MQ=60)
         TGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCC                                                                                                                            >  NB501016:228:HFK2FAFXY:3:21504:22630:3817/1‑151 (MQ=60)
                   CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTT                                                                                                                  >  NB501016:228:HFK2FAFXY:1:21209:12916:5166/1‑151 (MQ=60)
                   CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTT                                                                                                                  >  NB501016:228:HFK2FAFXY:3:11503:17331:3856/1‑151 (MQ=60)
                   CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTT                                                                                                                  >  NB501016:228:HFK2FAFXY:2:21105:22410:10243/1‑151 (MQ=60)
                   CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTT                                                                                                                  >  NB501016:228:HFK2FAFXY:2:11205:19519:12178/1‑151 (MQ=60)
                    CCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTC                                                                                                                 >  NB501016:228:HFK2FAFXY:2:11210:5943:15169/1‑151 (MQ=60)
                              TTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACT                                                                                                       >  NB501016:228:HFK2FAFXY:4:21402:19646:3720/1‑151 (MQ=60)
                               TCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTT                                                                                                      >  NB501016:228:HFK2FAFXY:3:21401:4480:6188/1‑151 (MQ=60)
                                CACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTC                                                                                                     <  NB501016:228:HFK2FAFXY:3:11611:7242:3705/151‑1 (MQ=60)
                                  CAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCA                                                                                                    >  NB501016:228:HFK2FAFXY:2:11207:4887:13365/1‑150 (MQ=60)
                                                GTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCA                                                                                        >  NB501016:228:HFK2FAFXY:4:11502:11544:4736/1‑148 (MQ=60)
                                                   GGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTT                                                                                                                         <  NB501016:228:HFK2FAFXY:2:11211:11570:13126/112‑1 (MQ=60)
                                                      GAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTC                                                                                                                                            <  NB501016:228:HFK2FAFXY:4:21505:21032:19071/90‑1 (MQ=60)
                                                                   CTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTC                                                                                                                                               >  NB501016:228:HFK2FAFXY:3:11605:12584:13793/1‑74 (MQ=60)
                                                                        TGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTC                                                              >  NB501016:228:HFK2FAFXY:3:11410:8436:12091/1‑150 (MQ=60)
                                                                           CCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGT                                                          >  NB501016:228:HFK2FAFXY:3:21509:22805:3785/1‑151 (MQ=60)
                                                                            CACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTagatcggaagagcacacgtctgactcc                                                           >  NB501016:228:HFK2FAFXY:4:21507:7375:5293/1‑122 (MQ=60)
                                                                                CTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCG                                                     >  NB501016:228:HFK2FAFXY:2:21209:21373:14168/1‑151 (MQ=60)
                                                                                CTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACC                                                       <  NB501016:228:HFK2FAFXY:3:11612:22412:13917/149‑1 (MQ=60)
                                                                                  TCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGTGTAGAAAAGTTATTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGT                                                   >  NB501016:228:HFK2FAFXY:2:21204:25812:6874/1‑151 (MQ=60)
                                                                                  TCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGACCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCG                                                     <  NB501016:228:HFK2FAFXY:4:21401:4111:1552/149‑1 (MQ=60)
                                                                                   CTTCGATCAGCGGTTGCAGCAAACGCATTTTTTCTTAGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGTTTCCACTTAAATATATTCGGCATTGATCATCTGCTCGAACGCCAGTTCACGTACCCGGT                                                   >  NB501016:228:HFK2FAFXY:2:11212:26104:14317/1‑150 (MQ=60)
                                                                                           AGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGC                                                                                                                              <  NB501016:228:HFK2FAFXY:1:11302:21751:13504/67‑1 (MQ=60)
                                                                                                 TGCAGCAAACGCATTTTTTCTTCGTTAGTAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGGTCGACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCG                                     <  NB501016:228:HFK2FAFXY:3:21603:26436:3947/150‑1 (MQ=60)
                                                                                                       AAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTACACGTACCCGGTACGAAAGCGTGGCGGAGAAC                               >  NB501016:228:HFK2FAFXY:1:21212:4753:15852/1‑150 (MQ=60)
                                                                                                        AACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAAC                               >  NB501016:228:HFK2FAFXY:1:21212:24436:10132/1‑149 (MQ=60)
                                                                                                                    CTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTG                      <  NB501016:228:HFK2FAFXY:1:21108:11346:1184/146‑1 (MQ=60)
                                                                                                                      TCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCG                 <  NB501016:228:HFK2FAFXY:3:21501:1483:19358/149‑1 (MQ=60)
                                                                                                                      TCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGGTTTCTGTTCCGGTTCCACTT                                                                                                      >  NB501016:228:HFK2FAFXY:2:11111:13241:17206/1‑64 (MQ=60)
                                                                                                                        GTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAG                    >  NB501016:228:HFK2FAFXY:4:21510:14527:9787/1‑144 (MQ=60)
                                                                                                                              AGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCA       <  NB501016:228:HFK2FAFXY:1:21106:20650:16209/151‑1 (MQ=60)
                                                                                                                              AGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGC        >  NB501016:228:HFK2FAFXY:3:21408:21816:12879/1‑150 (MQ=60)
                                                                                                                              AGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGC        >  NB501016:228:HFK2FAFXY:2:11304:26732:15097/1‑150 (MQ=60)
                                                                                                                              AGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGC        >  NB501016:228:HFK2FAFXY:2:11303:1285:15103/1‑150 (MQ=60)
                                                                                                                                  TAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAG      <  NB501016:228:HFK2FAFXY:1:21108:18461:16328/148‑1 (MQ=60)
                                                                                                                                   AGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG  <  NB501016:228:HFK2FAFXY:1:21210:19588:12901/151‑1 (MQ=60)
                                                                                                                                         GTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGG     <  NB501016:228:HFK2FAFXY:3:11608:22544:4054/142‑1 (MQ=60)
                                                                                                                                         GTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGG     <  NB501016:228:HFK2FAFXY:2:21203:12219:9914/142‑1 (MQ=60)
                                                                                                                                          |                                                                                                                                               
ACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG  >  NC_000913/3964791‑3965072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: