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
A3 F181 I0 R1 65 126.1 6126155 96.5% 5911739 99.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,598:1 +T 100% pseudogene (57/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,5971.T96.9% 453.8 / 6.1 129pseudogene (56/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (56/69);  total (57/72)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 137 total aligned reads displayed.

TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691
                                                                                                    |                                                                                              
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                 >  2:399256/1‑100 (MQ=255)
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                 <  2:2974929/100‑1 (MQ=255)
  tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑tt                                                                                              >  2:1573261/1‑99 (MQ=255)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttc                                                                                             >  1:1457413/1‑98 (MQ=255)
     aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttccc                                                                                           >  1:1112845/1‑96 (MQ=255)
     aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttccc                                                                                           >  1:2887061/1‑96 (MQ=255)
     aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttccc                                                                                           >  1:1196889/1‑96 (MQ=255)
     aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑Ttccc                                                                                           <  1:2483408/100‑2 (MQ=255)
      tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttcccg                                                                                          >  1:2200239/1‑95 (MQ=255)
      tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑Ttcccg                                                                                          <  2:713641/100‑3 (MQ=255)
       gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCgg                                                                                          <  2:2578213/100‑1 (MQ=255)
       gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCgg                                                                                          >  2:2557599/1‑100 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGt                                                                                         >  2:2414790/1‑100 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAACAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGt                                                                                         >  2:1540423/1‑100 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTg                                                                                        <  1:2771581/100‑1 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGAc                                                                                      >  1:2303230/1‑100 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGAc                                                                                      >  2:2053593/1‑100 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATTAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGCCGCt                                                                                   <  1:380592/100‑1 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACg                                                                                     >  1:621783/1‑98 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCt                                                                                   >  2:211212/1‑100 (MQ=255)
               aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTa                                                                                  <  2:1112851/100‑1 (MQ=255)
               aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTa                                                                                  >  1:684185/1‑100 (MQ=255)
                 cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATc                                                                                <  2:1977180/100‑1 (MQ=255)
                 cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATc                                                                                <  1:413766/100‑1 (MQ=255)
                     aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcg                                                                            >  1:1156378/1‑100 (MQ=255)
                       tACATCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcgcg                                                                          <  2:1782465/100‑1 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACgcgc                                                                         <  2:474478/100‑1 (MQ=255)
                         cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGcc                                                                        >  1:97660/1‑100 (MQ=255)
                         cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGcc                                                                        >  1:97511/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgt                                                                      >  1:1648506/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgt                                                                      >  1:861188/1‑100 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgtt                                                                     >  1:1763699/1‑100 (MQ=255)
                             agcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttg                                                                    >  2:2911388/1‑100 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgt                                                                   <  2:2007677/100‑1 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgt                                                                   <  2:1763143/100‑1 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgt                                                                   <  2:489426/100‑1 (MQ=255)
                                cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttg                                                                 >  1:1889224/1‑100 (MQ=255)
                                cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttg                                                                 <  1:257469/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                >  2:2404485/1‑100 (MQ=255)
                                  ttAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgtt                                                               <  2:2942350/100‑1 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                                                             >  1:2195303/1‑100 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                                                              <  1:648828/99‑1 (MQ=255)
                                     aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                            <  2:968670/100‑1 (MQ=255)
                                      aTCAATGACTGGCAGCAACACCGCGCGCCGCTGCGTGCTGAACATTCCTGGCGTCCCGCCCATTCCCGGTGACGCTACCTACCCGCCgtcgtt                                                                  >  1:2644319/1‑93 (MQ=255)
                                       tCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaaa                                                           <  1:1239932/99‑1 (MQ=255)
                                        cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                                                         >  2:2128575/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                        >  2:2661387/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                        >  2:1368438/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                          <  1:670125/98‑1 (MQ=255)
                                                ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg                                                 >  2:2200609/1‑100 (MQ=255)
                                                ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg                                                 >  2:2285133/1‑100 (MQ=255)
                                                   agcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc                                              >  2:2010624/1‑100 (MQ=255)
                                                   agcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc                                              <  1:1567266/100‑1 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  1:71251/100‑1 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  2:2702739/100‑1 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  1:72678/100‑1 (MQ=255)
                                                         acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACACCTGTCGGATCGTAAAc                                        >  2:2989491/1‑100 (MQ=255)
                                                         acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc                                        <  1:511521/100‑1 (MQ=255)
                                                         acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc                                        >  2:362779/1‑100 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       <  1:58927/100‑1 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       >  1:2996918/1‑100 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       <  1:155917/100‑1 (MQ=255)
                                                           tGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt                                      <  2:900899/100‑1 (MQ=255)
                                                           tGCGCGCAGCGGCGTGATGAACATTCCTGGCGCTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       <  1:302728/99‑1 (MQ=255)
                                                              gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg                                   <  1:2911370/100‑1 (MQ=255)
                                                              gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg                                   >  1:968087/1‑100 (MQ=255)
                                                               cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                  <  2:937301/100‑1 (MQ=255)
                                                                   ncTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACACTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTg                              >  1:1097440/2‑100 (MQ=255)
                                                                   gCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTg                              <  1:388761/100‑1 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                            <  2:1481640/100‑1 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACACTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                            <  2:1097446/100‑1 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgt                           >  1:2122991/1‑100 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                          <  2:2391658/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         <  2:1713944/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         <  2:926111/100‑1 (MQ=255)
                                                                         tgatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcgt                        <  2:651271/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                       <  2:352617/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                       <  1:2758573/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                       >  1:1642832/1‑100 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                      <  1:1978891/100‑1 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                      <  1:2301174/100‑1 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                      <  2:1901530/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                    <  2:1864821/100‑1 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                  <  1:1887151/100‑1 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAg                <  1:687626/100‑1 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAg                <  1:2592884/100‑1 (MQ=255)
                                                                                   tCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAg                <  1:1005641/98‑1 (MQ=255)
                                                                                   tCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa               <  1:1962789/99‑1 (MQ=255)
                                                                                   tCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAt              >  2:1154192/1‑100 (MQ=255)
                                                                                       ggCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg          <  2:968092/100‑1 (MQ=255)
                                                                                        gCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTggg         <  2:673624/100‑1 (MQ=255)
                                                                                        gCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTggg         <  2:1867496/100‑1 (MQ=255)
                                                                                         cGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgggg        >  1:3012918/1‑100 (MQ=255)
                                                                                         cGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgggg        <  2:2557571/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       >  1:2878019/1‑100 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:2882221/100‑1 (MQ=255)
                                                                                           ttACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      <  1:56937/100‑1 (MQ=255)
                                                                                            tACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcag     >  2:1319945/1‑100 (MQ=255)
                                                                                            tACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcag     >  2:2312727/1‑100 (MQ=255)
                                                                                               gACCA‑TCCCGGTGACGCTATCTACGAGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAc  <  2:2681172/99‑1 (MQ=255)
                                                                                                    |                                                                                              
TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691

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

GATK/CNVnator alignment

N/A