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
A21 F39 I0 R1 217 77.2 2771496 98.2% 2721609 143.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,151,873 T→A 6.2% intergenic (+22/‑66) acpP → / → fabF acyl carrier protein (ACP)/3‑oxoacyl‑[acyl‑carrier‑protein] synthase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,8730TA 6.2% 250.1 / 5.4 76intergenic (+22/‑66)acpP/fabFacyl carrier protein (ACP)/3‑oxoacyl‑[acyl‑carrier‑protein] synthase II
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base A (3/2);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAG  >  NC_000913/1151727‑1152017
                                                                                                                                                  |                                                                                                                                                
ttGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                >  1:132703/1‑149 (MQ=255)
ttGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                >  1:564498/1‑149 (MQ=255)
       cGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGcc                                                                                                                                         >  1:143551/1‑149 (MQ=255)
        gTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCt                                                                                                                                        <  3:429842/149‑1 (MQ=255)
        gTTGAGCTGGGAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCt                                                                                                                                        <  4:435966/149‑1 (MQ=255)
          tGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGa                                                                                                                                      <  3:584372/149‑1 (MQ=255)
               tGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGtttt                                                                                                                                 <  2:651358/149‑1 (MQ=255)
                  tAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATc                                                                                                                              >  1:614983/1‑149 (MQ=255)
                   aaTGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCt                                                                                                                             >  4:397281/1‑149 (MQ=255)
                     tGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCttt                                                                                                                           >  3:346170/1‑149 (MQ=255)
                     tGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTCTTATCttt                                                                                                                           <  4:60814/149‑1 (MQ=255)
                         tctGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGt                                                                                                                       >  4:440129/1‑149 (MQ=255)
                           tGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTcc                                                                                                                     >  2:357747/1‑149 (MQ=255)
                           tGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTcc                                                                                                                     >  2:275314/1‑149 (MQ=255)
                                        tGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCAACCAGGCGTAAGTGAACATCTCCAGGCGGTCGCTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAtt                                                                                                        >  4:443520/1‑149 (MQ=255)
                                              tGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtcc                                                                                                  >  1:387970/1‑149 (MQ=255)
                                                    tCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTccctccct                                                                                            <  4:664219/149‑1 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                         <  2:523908/149‑1 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                         >  2:390290/1‑149 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                         <  2:176784/149‑1 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                         <  4:524885/149‑1 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                         >  4:510912/1‑149 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGAGTACATCAACGGCCACCAGGCGTAAGTTAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                         >  4:248658/1‑149 (MQ=255)
                                                         aCGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACTGCCTGAGTTTTATCtttt                                                                                                                          >  2:624859/1‑114 (MQ=255)
                                                         aCGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACTGCCTGAGTTTTATCtttt                                                                                                                          <  1:624859/114‑1 (MQ=255)
                                                           gaagaaGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCtttt                                                                                                                          <  1:510452/112‑1 (MQ=255)
                                                           gaagaaGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCtttt                                                                                                                          >  2:510452/1‑112 (MQ=255)
                                                              gaagCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACaaa                                                                                  <  2:387970/149‑1 (MQ=255)
                                                               aagCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAAc                                                                                 <  2:331342/149‑1 (MQ=255)
                                                                 gCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgt                                                                               >  3:310414/1‑149 (MQ=255)
                                                                       aaaaTCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCt                                                                                                                                        >  4:492774/1‑86 (MQ=255)
                                                                       aaaaTCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCt                                                                                                                                        <  3:492774/86‑1 (MQ=255)
                                                                         aaTCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGc                                                                       >  1:151248/1‑149 (MQ=255)
                                                                          ctCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcg                                                                      >  3:443475/2‑149 (MQ=255)
                                                                          aTCACCACCGTTCAGGCTGCCATTGATTACATCAACGGGCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcg                                                                      >  4:282340/1‑149 (MQ=255)
                                                                           tCACCACCGTTCAGGCTGCCATGGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCTCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcgt                                                                     >  2:53812/1‑149 (MQ=255)
                                                                            caccacCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcgtc                                                                    <  1:275314/149‑1 (MQ=255)
                                                                              ccaccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcgtcgt                                                                  >  1:3619/1‑149 (MQ=255)
                                                                              ccaccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcgtcgt                                                                  <  4:471948/149‑1 (MQ=255)
                                                                              ccaccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcgtcgt                                                                  >  1:10824/1‑149 (MQ=255)
                                                                              ccaccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGcgtcgt                                                                  <  2:614983/149‑1 (MQ=255)
                                                                                   gTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGt                                                             >  1:270390/1‑149 (MQ=255)
                                                                                      cAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTtgt                                                          <  1:390290/149‑1 (MQ=255)
                                                                                      cAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTtgt                                                          <  3:530904/149‑1 (MQ=255)
                                                                                           tGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAg                                                                        >  1:222413/1‑130 (MQ=255)
                                                                                           tGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAg                                                                        <  2:222413/130‑1 (MQ=255)
                                                                                             ccATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGa                                                   >  1:370699/1‑149 (MQ=255)
                                                                                                   aTTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGc                                             >  2:71439/1‑149 (MQ=255)
                                                                                                        aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgtt                                        >  3:32113/1‑149 (MQ=255)
                                                                                                        aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgtt                                        >  4:639433/1‑149 (MQ=255)
                                                                                                         tCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgttg                                       >  4:514016/1‑149 (MQ=255)
                                                                                                          cAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgttgt                                      <  2:370699/149‑1 (MQ=255)
                                                                                                           aaCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtc                                     <  2:447000/149‑1 (MQ=255)
                                                                                                            aCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtct                                    >  3:349291/1‑149 (MQ=255)
                                                                                                            aCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtct                                    <  1:357747/149‑1 (MQ=255)
                                                                                                            aCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtct                                    >  1:663791/1‑149 (MQ=255)
                                                                                                                  accaGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTc                              >  3:483504/1‑149 (MQ=255)
                                                                                                                    caGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCgg                            >  4:318039/1‑149 (MQ=255)
                                                                                                                     aGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGc                           >  4:273127/1‑149 (MQ=255)
                                                                                                                       gCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCaa                         <  3:510912/149‑1 (MQ=255)
                                                                                                                       gCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCaa                         <  1:367917/149‑1 (MQ=255)
                                                                                                                             gTGAACATCTTCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCGCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGt                   <  1:92022/149‑1 (MQ=255)
                                                                                                                             gTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGt                   >  1:304967/1‑149 (MQ=255)
                                                                                                                             gTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGt                   >  3:182915/1‑149 (MQ=255)
                                                                                                                              tGAACATCTCCAGGCGGTCGTTCGACCGCCTTAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGGCAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGACGGCAATACCGTa                  >  4:398520/1‑149 (MQ=255)
                                                                                                                              tGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTa                  <  4:349291/149‑1 (MQ=255)
                                                                                                                               gAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTag                 >  4:436717/1‑149 (MQ=255)
                                                                                                                                aaCATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTaga                <  3:440129/149‑1 (MQ=255)
                                                                                                                                  cATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACCGGGCATGTTGTCTCCTGTCGGCAATACTGTAGAGt              <  2:143551/149‑1 (MQ=255)
                                                                                                                                      gccAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTAc          <  2:132703/148‑1 (MQ=255)
                                                                                                                                      ctcAGGCGGTCGAACGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTAc          >  1:557794/3‑149 (MQ=255)
                                                                                                                                       tcAGGCGGTCGAACGACCGCCTGGGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTAcc         >  2:358647/2‑149 (MQ=255)
                                                                                                                                       ccAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTAcc         <  1:53812/149‑1 (MQ=255)
                                                                                                                                        taagCGGTGGTTCGACCGCCTGAGTTATATCTTTTTGTCCCACCAGAATCATTTTTTCCCTCCCTGGAGGACATACGTGTCTAAGCGTCGTGTAGTTGTGACCGCACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCt        <  4:218210/146‑1 (MQ=255)
                                                                                                                                         aGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTg       <  4:167458/149‑1 (MQ=255)
                                                                                                                                              gTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAg  >  1:608760/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAG  >  NC_000913/1151727‑1152017

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A