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
A4 F244 I0 R1 270 114.4 5840133 93.3% 5448844 99.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,559:1 +A 100% pseudogene (18/663 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,5591.A86.7% 332.7 / 46.3 113pseudogene (18/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (5/10);  new base A (49/49);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGC  >  NC_000913/3951460‑3951657
                                                                                                      |                                                                                                   
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑c                                                                                                      >  2:2218963/1‑100 (MQ=255)
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑c                                                                                                      <  1:2096921/100‑1 (MQ=255)
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑c                                                                                                      >  2:2456137/1‑100 (MQ=255)
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑ca                                                                                                     >  2:568491/1‑99 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑cagcg                                                                                                  <  2:44251/100‑5 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCGGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑cagcg                                                                                                  >  2:131635/1‑96 (MQ=255)
     gCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑g‑cagcgc                                                                                                 <  2:49174/100‑6 (MQ=255)
     gCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGACGATAAATCAATGACTGGCAGCAACACT‑g‑cagcgc                                                                                                 >  1:1281028/1‑95 (MQ=255)
           tCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGc                                                                                           >  2:1730986/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGTTGCGtg                                                                                        <  1:451149/100‑1 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGtg                                                                                        <  2:1817212/100‑1 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGtg                                                                                        >  2:1124995/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGtg                                                                                        >  2:2394052/1‑100 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGtgat                                                                                      <  2:665813/100‑1 (MQ=255)
                 acttgttGCATTACAATGTGATTTAAATGCTCTGTTACCAGCAGTACAACAGCCTTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGtgat                                                                                      <  1:422258/96‑1 (MQ=255)
                  caTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGAAACACT‑G‑CAGCGCAGCTGCGtgatga                                                                                    <  1:970198/100‑1 (MQ=255)
                     gtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACa                                                                                 >  1:1712939/1‑100 (MQ=255)
                     gtgGCATTACAAGGTCATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACa                                                                                 <  2:2876367/100‑1 (MQ=255)
                      tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCATCTGCGTGATGAACAt                                                                                >  2:2870053/1‑100 (MQ=255)
                       ggCATTACAAGGAGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACAtt                                                                               >  2:1912962/1‑100 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTcc                                                                             >  1:121036/1‑100 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTcc                                                                             <  1:133918/100‑1 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCt                                                                            >  2:1464792/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCCATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTAATGAACATTCCTgg                                                                          >  1:1293230/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTgg                                                                          <  2:1202652/100‑1 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGc                                                                         >  2:906222/1‑100 (MQ=255)
                              cgaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        <  1:1602830/98‑1 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        >  2:498466/1‑100 (MQ=255)
                               aaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGtt                                                                      >  2:1738439/1‑100 (MQ=255)
                               aaGGTGATTTAAATGCTCTGTTACCAGCACCACAGCAGCCGTTAAATCAATGACTGGCAGCAACACC‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGt                                                                       <  1:1756450/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                     <  1:2864690/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                     <  1:607713/100‑1 (MQ=255)
                                  gTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                    <  1:1227089/100‑1 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                 <  1:782892/100‑1 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCTCCGCTGCATGATGAACATTCCTGGCGTTACGAcc                                                                >  1:2753928/1‑100 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                                >  1:1503115/1‑100 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                                >  1:2577062/1‑100 (MQ=255)
                                        tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                              >  1:1728280/1‑100 (MQ=255)
                                        tAAATGCTCTGTTACCAGCATTACAGAAGCCGTTAAATCAATCACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                              <  2:2775083/100‑1 (MQ=255)
                                           aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCTGCTGCGTGATGAACATTCCTGGCGTTACGACCATccc                                                           <  1:1556300/100‑1 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                         <  2:2448306/100‑1 (MQ=255)
                                             gCTCTGTTACCACCATTACAGCAGCCGTTCAATCAATCCCTCGCACCACCACT‑C‑CAGCCCAGCAGCGTCACGAACCTTCCTGGCGTCACGACCGTCCCg                                                          >  1:66875/1‑99 (MQ=255)
                                              ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        >  2:1232813/1‑100 (MQ=255)
                                                 tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     >  1:678311/1‑100 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGCCCATCCCGGTGACGc                                                   >  2:2231818/1‑100 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                   <  2:426516/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                                >  2:2203692/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                               <  2:1446745/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                               >  2:1099545/1‑100 (MQ=255)
                                                       cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  2:2428478/1‑100 (MQ=255)
                                                        aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  1:1483625/1‑100 (MQ=255)
                                                        aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  2:1449282/1‑100 (MQ=255)
                                                          ccatacacaagcCGTAATATCACTCACTGGCAGCAACACC‑GCC‑GCGCATATGCGTGATGAACATTCCACGCGTAACGCCCATCCCGCTGACGCTATc                                               <  1:1308044/88‑1 (MQ=255)
                                                           aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                           <  2:566024/100‑1 (MQ=255)
                                                            ttACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgc                                          >  1:2542800/1‑100 (MQ=255)
                                                             tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgcg                                         <  2:2592638/100‑1 (MQ=255)
                                                               cagcagCCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                       >  2:77090/1‑100 (MQ=255)
                                                                 gcagcCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  1:1029789/100‑1 (MQ=255)
                                                                   agcCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                   <  2:1863621/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  <  1:1349405/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  <  1:1349350/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACT‑C‑CAGCGCAGCTGCGTGATGAACATCCCTGGCGTTACGACCATCCCGGTGCCGCTATCTACGCGCCgttg                                   >  1:2562345/1‑99 (MQ=255)
                                                                     ccGTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 >  2:491419/1‑100 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGG‑A‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 <  2:1094469/98‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  2:815789/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  2:2003859/100‑1 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              <  1:1730978/100‑1 (MQ=255)
                                                                         tAAATCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             >  2:2544607/1‑100 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                            <  2:2629428/100‑1 (MQ=255)
                                                                          aaaTCAATGACTAGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             >  2:841132/1‑99 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                          >  2:1521799/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                          <  2:29262/100‑1 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                          <  2:206297/100‑1 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                          >  2:1276516/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             <  1:424260/98‑1 (MQ=255)
                                                                              cAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aaca                       >  1:602327/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aac                        <  1:1739581/99‑1 (MQ=255)
                                                                               aaTGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCGCGGTGACGCTATCTGCGCGCCGTTGTTGTTAA‑aacaa                      >  1:1963253/1‑100 (MQ=255)
                                                                                aTGACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaac                     <  1:691761/100‑1 (MQ=255)
                                                                                  gACTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTT‑TTaacaacaacgg                   <  1:988805/100‑3 (MQ=255)
                                                                                  gACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGt                  <  1:2231804/100‑1 (MQ=255)
                                                                                   aCTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTg                   <  1:2655793/99‑1 (MQ=255)
                                                                                   aCTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTc                 >  1:2219056/1‑100 (MQ=255)
                                                                                    cTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTT‑TTaacaacaacggcg                 <  2:1404191/100‑5 (MQ=255)
                                                                                    cTGGCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTT‑TTaacaacaacggcg                 <  1:1512116/100‑5 (MQ=255)
                                                                                      ggCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTTT‑aacaacaacggc               <  1:896793/100‑7 (MQ=255)
                                                                                      ggCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCgg               >  1:2524719/1‑100 (MQ=255)
                                                                                      ggCAGCAACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCgg               >  2:2356298/1‑100 (MQ=255)
                                                                                       gcagcaACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGa              >  1:1601644/1‑100 (MQ=255)
                                                                                        cagcaACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGAt             <  2:1642286/100‑1 (MQ=255)
                                                                                         agcaACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATc            >  2:173193/1‑100 (MQ=255)
                                                                                          gcagcaCT‑G‑CAGCGCAGATACGCGATGCACCCTCCTAAGGTTACGACAATCCCGGTGACGCTATCTACGCGACGTTGTTGTTAA‑AACAACTGTGGGATCg           <  1:966797/100‑1 (MQ=255)
                                                                                           caACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGt          >  1:2802037/1‑100 (MQ=255)
                                                                                           caACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGt          <  2:1269992/100‑1 (MQ=255)
                                                                                           caACACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGt          <  2:891249/100‑1 (MQ=255)
                                                                                            aaCACT‑G‑CAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCg           >  2:2617478/1‑98 (MQ=255)
                                                                                            aaCACT‑G‑CAGCGCAGCTGCGTGATCAACATTCCTGGCTTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACACCTGTCGGATCGTa         <  1:1799145/100‑1 (MQ=255)
                                                                                             acacT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTaaa       <  2:2540573/100‑1 (MQ=255)
                                                                                               acT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAAcc     <  1:1664745/100‑1 (MQ=255)
                                                                                                   g‑c‑gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGc  >  1:1564373/1‑100 (MQ=255)
                                                                                                      |                                                                                                   
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACT‑G‑C‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGC  >  NC_000913/3951460‑3951657

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

GATK/CNVnator alignment

N/A