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 F28 I1 R1 120 93.5 3664254 88.3% 3235536 136.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,819,492 G→A G840G (GGC→GGT)  alaS ← alanine‑‑tRNA ligase/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,819,4920GA98.0% 323.7 / ‑2.2 99G840G (GGC→GGT) alaSalanine‑‑tRNA ligase/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (41/56);  total (42/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

Only 100 of 102 total aligned reads displayed.

CGGCATTGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAG  >  NC_000913/2819360‑2819623
                                                                                                                                    |                                                                                                                                   
cGGCATTGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGccacccac                                                                                                                                <  1:1182265/138‑1 (MQ=255)
 ggCATTGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCcacccacc                                                                                                                               <  2:1760626/138‑1 (MQ=255)
   cATTGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACctg                                                                                                                             <  1:1494176/138‑1 (MQ=255)
    aTTGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACctgc                                                                                                                            <  2:277223/138‑1 (MQ=255)
      tGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCATCCACctgctg                                                                                                                          <  2:1830058/138‑1 (MQ=255)
         cTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGc                                                                                                                       <  2:1117282/138‑1 (MQ=255)
         cTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGc                                                                                                                       <  1:1441850/138‑1 (MQ=255)
         cTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGc                                                                                                                       <  1:1421286/138‑1 (MQ=255)
              gCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCa                                                                                                                  <  2:1214850/138‑1 (MQ=255)
              gCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCa                                                                                                                  <  1:755786/138‑1 (MQ=255)
              gCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCa                                                                                                                  <  1:861315/138‑1 (MQ=255)
                   tattatTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCa                                                                                                             <  2:1667839/138‑1 (MQ=255)
                   tattatTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCa                                                                                                             >  1:301579/1‑138 (MQ=255)
                      tattGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATc                                                                                                          <  2:1182906/138‑1 (MQ=255)
                        ttGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAg                                                                                                        >  2:1733180/1‑138 (MQ=255)
                         tGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGt                                                                                                       >  2:461201/1‑138 (MQ=255)
                          gCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGtt                                                                                                      >  1:1670673/1‑138 (MQ=255)
                          gCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGtt                                                                                                      >  2:735/1‑138 (MQ=255)
                          gCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGtt                                                                                                      >  2:1212329/1‑138 (MQ=255)
                           cAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTc                                                                                                     <  2:1192193/138‑1 (MQ=255)
                             aTTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTccc                                                                                                   <  1:458540/138‑1 (MQ=255)
                              tttCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTcccc                                                                                                  >  2:709122/1‑138 (MQ=255)
                               ttCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCt                                                                                                 >  2:276542/1‑138 (MQ=255)
                               ttCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCt                                                                                                 >  1:1307137/1‑138 (MQ=255)
                               ttCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCt                                                                                                 >  2:752977/1‑138 (MQ=255)
                                tCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTg                                                                                                >  2:769557/1‑138 (MQ=255)
                                  gcgcTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCt                                                                                              <  1:276542/138‑1 (MQ=255)
                                  gcgcTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCt                                                                                              <  1:1167323/138‑1 (MQ=255)
                                     cTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTc                                                                                           <  1:709122/138‑1 (MQ=255)
                                     cTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTc                                                                                           <  2:820224/138‑1 (MQ=255)
                                     cTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTc                                                                                           >  2:1609017/1‑138 (MQ=255)
                                     cTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTc                                                                                           >  1:1344746/1‑138 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          <  2:1307137/138‑1 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          >  2:1323591/1‑138 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          <  2:1670673/138‑1 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          <  2:301579/138‑1 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          <  1:1212329/138‑1 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          <  1:675976/138‑1 (MQ=255)
                                      tGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTca                                                                                          <  1:607829/138‑1 (MQ=255)
                                          ccAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACg                                                                                      >  1:109240/1‑138 (MQ=255)
                                          ccAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACg                                                                                      <  1:1733180/138‑1 (MQ=255)
                                          ccAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACg                                                                                      <  1:461201/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  1:1607655/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  1:1353842/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  2:109240/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  1:1323591/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  1:769557/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  1:730833/138‑1 (MQ=255)
                                                tttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgt                                                                                <  1:702905/138‑1 (MQ=255)
                                                 ttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgtg                                                                               <  1:1419003/138‑1 (MQ=255)
                                                 ttCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgtg                                                                               <  2:343419/138‑1 (MQ=255)
                                                  tCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGa                                                                              <  2:1344746/138‑1 (MQ=255)
                                                    acacTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACg                                                                            >  1:1181324/1‑138 (MQ=255)
                                                          gCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTa                                                                      >  2:675750/1‑138 (MQ=255)
                                                          gCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGACATGTCAGGACGTCCACCACCCTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtg                                                                                 >  2:1360039/1‑127 (MQ=255)
                                                          gCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGACATGTCAGGACGTCCACCACCCTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtg                                                                                 <  1:1360039/127‑1 (MQ=255)
                                                           cTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAg                                                                     <  1:390845/138‑1 (MQ=255)
                                                           cTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAg                                                                     >  2:200985/1‑138 (MQ=255)
                                                            tAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGa                                                                    >  1:47196/1‑138 (MQ=255)
                                                             aaCGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgtg                                                                               >  2:1558613/1‑126 (MQ=255)
                                                             aaCGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCtgtg                                                                               <  1:1558613/126‑1 (MQ=255)
                                                             aaCGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGAt                                                                   >  1:1208530/1‑138 (MQ=255)
                                                               cGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATAc                                                                 >  1:1776722/1‑138 (MQ=255)
                                                                gCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACg                                                                <  1:675750/138‑1 (MQ=255)
                                                                    cAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTg                                                            <  2:1814824/138‑1 (MQ=255)
                                                                      ggTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCa                                                          <  1:215697/138‑1 (MQ=255)
                                                                      ggTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCa                                                          >  2:1805564/1‑138 (MQ=255)
                                                                      ggTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCa                                                          >  1:1372675/1‑138 (MQ=255)
                                                                      ggTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCa                                                          >  1:24806/1‑138 (MQ=255)
                                                                      ggTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCa                                                          >  1:1018174/1‑138 (MQ=255)
                                                                           ggCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCag                                                     >  2:1280850/1‑138 (MQ=255)
                                                                            gCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCaga                                                    <  2:47196/138‑1 (MQ=255)
                                                                            gCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCaga                                                    >  1:875082/1‑138 (MQ=255)
                                                                            gCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTACACCTCCCATGCCACCCACCTGCTGAGCGACCATGCCAATCAGTTCCCCTGCTTTCACGCGATCTGTGACGTCCTTAGATGCGCCTGCAATCaga                                                    <  2:291617/138‑1 (MQ=255)
                                                                              cGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAgaga                                                  >  1:660098/1‑138 (MQ=255)
                                                                                cagcaTCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGaaa                                                <  1:1609017/138‑1 (MQ=255)
                                                                                cagcaTCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTT‑CCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAac                                               >  2:1240148/1‑138 (MQ=255)
                                                                                  gcaTCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAacc                                              >  1:623993/1‑138 (MQ=255)
                                                                                   caTCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAacct                                             <  1:200985/138‑1 (MQ=255)
                                                                                    aTCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAacctt                                            <  1:488468/138‑1 (MQ=255)
                                                                                     tCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAacctta                                           <  2:1018174/138‑1 (MQ=255)
                                                                                          accaccGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTc                                      <  2:623993/138‑1 (MQ=255)
                                                                                          accaccGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTc                                      <  1:735/138‑1 (MQ=255)
                                                                                                 gCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTAcc                               >  2:1250972/1‑138 (MQ=255)
                                                                                                   ttGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGTGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGt                             >  1:1528893/1‑138 (MQ=255)
                                                                                                   ttGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGt                             >  1:89150/1‑138 (MQ=255)
                                                                                                    tGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGtt                            <  2:875082/138‑1 (MQ=255)
                                                                                                    tGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGtt                            <  1:1805564/138‑1 (MQ=255)
                                                                                                    tGCGCCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGtt                            <  1:1250972/138‑1 (MQ=255)
                                                                                                     gcgcCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTg                           <  1:1412049/138‑1 (MQ=255)
                                                                                                     gcgcCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTg                           >  1:1219111/1‑138 (MQ=255)
                                                                                                       gcCATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGcc                         <  2:1372675/138‑1 (MQ=255)
                                                                                                         cATGTCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAg                       <  2:1208530/138‑1 (MQ=255)
                                                                                                             tCAGGACGTCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACg                   >  2:403122/1‑138 (MQ=255)
                                                                                                                  aCGTCCACCACCCTTGCCACCCACCTGCTGTGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAAt              <  2:1528893/138‑1 (MQ=255)
                                                                                                                     tCCACCACCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGa                  <  1:999559/131‑1 (MQ=255)
                                                                                                                       caccacCCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCg         >  1:1727696/1‑138 (MQ=255)
                                                                                                                         ccaccCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGAc       >  1:631390/1‑138 (MQ=255)
                                                                                                                         ccaccCTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGAc       >  1:584747/1‑138 (MQ=255)
                                                                                                                              cTTGCCACCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAg  >  2:1668687/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                   
CGGCATTGCCTGACGCTTATATTATTGCAATTTCGCGCTGACCCAGCCTTTCACACTGGCTAACGCTGCAGGTAAGGCCGCAGCATCCGTACCACCGGCTTGCGCCATGTCAGGACGTCCACCACCCTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAG  >  NC_000913/2819360‑2819623

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

GATK/CNVnator alignment

N/A