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 F13 I1 R1 108 324.7 15372274 99.7% 15326157 98.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,890,657 G→A G414G (GGC→GGT)  cysJ ← sulfite reductase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,6570GA99.7% 1247.1 / NA 306G414G (GGC→GGT) cysJsulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (184/121);  minor base C (0/1);  total (184/122)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

AGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCG  >  NC_000913/2890559‑2890754
                                                                                                  |                                                                                                 
aGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGAcc                                                                                                 >  1:4705686/1‑101 (MQ=255)
   cGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTc                                                                                              >  2:4705690/1‑101 (MQ=255)
   cGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTc                                                                                              <  2:4705337/101‑1 (MQ=255)
    ggAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTcc                                                                                             >  2:4705693/1‑101 (MQ=255)
     gAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCa                                                                                            <  2:4705370/101‑1 (MQ=255)
      aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAc                                                                                           >  1:4705698/1‑101 (MQ=255)
      aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAc                                                                                           <  1:4705145/101‑1 (MQ=255)
      aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAc                                                                                           >  2:4705700/1‑101 (MQ=255)
      aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAc                                                                                           <  1:4705304/101‑1 (MQ=255)
       aGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcg                                                                                          >  1:4705706/1‑101 (MQ=255)
        gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgt                                                                                         <  1:4705469/101‑1 (MQ=255)
        gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgt                                                                                         <  1:4704988/101‑1 (MQ=255)
         ttATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgtc                                                                                        <  2:4705362/101‑1 (MQ=255)
         ttATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgtc                                                                                        >  2:4705708/1‑101 (MQ=255)
          tATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgtcg                                                                                       <  1:4705201/101‑1 (MQ=255)
           aTCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgtcgt                                                                                      <  1:4705386/101‑1 (MQ=255)
           aTCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgtcgt                                                                                      >  1:4705715/1‑101 (MQ=255)
           aTCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCAcgtcgt                                                                                      >  2:4705714/1‑101 (MQ=255)
            tCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTa                                                                                     >  2:4705718/1‑101 (MQ=255)
              gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTaa                                                                                    <  1:4705069/100‑1 (MQ=255)
              gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAAc                                                                                   >  1:4705721/1‑101 (MQ=255)
               ttATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACg                                                                                  >  1:4705726/1‑101 (MQ=255)
               ttATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACg                                                                                  <  2:4705293/101‑1 (MQ=255)
               ttATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACg                                                                                  >  2:4705727/1‑101 (MQ=255)
                tATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGc                                                                                 >  1:4705731/1‑101 (MQ=255)
                  tGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcac                                                                               >  1:4705737/1‑101 (MQ=255)
                     tCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccac                                                                            >  1:4705741/1‑101 (MQ=255)
                      cGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccaca                                                                           >  2:4705744/1‑101 (MQ=255)
                       gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccacac                                                                          <  2:4705552/101‑1 (MQ=255)
                        aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccacacc                                                                         >  1:4705749/1‑101 (MQ=255)
                         tAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccacacca                                                                        >  2:4705759/1‑101 (MQ=255)
                         tAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccacacca                                                                        >  1:4705757/1‑101 (MQ=255)
                         tAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccacacca                                                                        >  1:4705753/1‑101 (MQ=255)
                           aaTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAAc                                                                      >  1:4705762/1‑101 (MQ=255)
                            aTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGCGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcac                                                                               >  1:4705766/1‑91 (MQ=255)
                             tACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACgg                                                                    >  2:4705768/1‑101 (MQ=255)
                              aCGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGt                                                                   >  1:4705770/1‑101 (MQ=255)
                              aCGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGt                                                                   <  2:4705415/101‑1 (MQ=255)
                               cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTg                                                                  <  1:4705285/101‑1 (MQ=255)
                                gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGa                                                                 <  2:4705352/101‑1 (MQ=255)
                                gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGa                                                                 >  2:4705773/1‑101 (MQ=255)
                                 cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTg                                                                  >  1:4705779/1‑99 (MQ=255)
                                 cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                >  2:4705778/1‑101 (MQ=255)
                                   gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACgt                                                              <  1:4705379/101‑1 (MQ=255)
                                     cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACgtgt                                                            >  2:4705784/1‑101 (MQ=255)
                                      ttCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTa                                                           >  1:4705789/1‑101 (MQ=255)
                                        cGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACt                                                         >  2:4705793/1‑101 (MQ=255)
                                         gCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACtt                                                        >  1:4705794/1‑101 (MQ=255)
                                         gCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACtt                                                        >  2:4705796/1‑101 (MQ=255)
                                          cccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTc                                                       >  2:4705800/1‑101 (MQ=255)
                                           ccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCg                                                      >  2:4705804/1‑101 (MQ=255)
                                            ctctTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGt                                                     <  2:4705493/101‑1 (MQ=255)
                                             tcttctTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTg                                                                  >  1:4705807/1‑87 (MQ=255)
                                              cttcttCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGtt                                                    <  2:4705433/100‑1 (MQ=255)
                                              cttcttCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                   >  2:4705809/1‑101 (MQ=255)
                                               ttcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtct                                                  <  2:4705401/101‑1 (MQ=255)
                                                tcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                 <  2:4705250/101‑1 (MQ=255)
                                                tcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                 <  2:4705216/101‑1 (MQ=255)
                                                tcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                 <  1:4705330/101‑1 (MQ=255)
                                                 cttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTg                                                                  <  2:4704762/83‑1 (MQ=255)
                                                 cttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                <  2:4705240/101‑1 (MQ=255)
                                                 cttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                <  2:4705640/101‑1 (MQ=255)
                                                 cttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                >  1:4705813/1‑101 (MQ=255)
                                                  ttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGa                                               >  2:4705815/1‑101 (MQ=255)
                                                      aCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTc                                           >  2:4705819/1‑101 (MQ=255)
                                                       cgcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGcaccac                                                                            >  1:4705820/1‑67 (MQ=255)
                                                       cgcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                <  2:4705184/95‑1 (MQ=255)
                                                          ggTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGcc                                       >  1:4705824/1‑101 (MQ=255)
                                                            tCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTg                                     <  2:4705316/101‑1 (MQ=255)
                                                            tCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTg                                     >  1:4705827/1‑101 (MQ=255)
                                                             cAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGc                                    >  1:4705830/1‑101 (MQ=255)
                                                              aGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTc                                           >  2:4705836/1‑93 (MQ=255)
                                                              aGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCg                                   >  2:4705833/1‑101 (MQ=255)
                                                               gCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCga                                  <  2:4705147/101‑1 (MQ=255)
                                                                cGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgag                                 >  1:4705840/1‑101 (MQ=255)
                                                                 gAGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgagg                                <  1:4705278/101‑1 (MQ=255)
                                                                  aGGAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgagga                               >  2:4705842/1‑101 (MQ=255)
                                                                   ggAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgaggag                              <  2:4705412/101‑1 (MQ=255)
                                                                    gAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                             <  2:4705084/101‑1 (MQ=255)
                                                                    gAAGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                             <  1:4705423/101‑1 (MQ=255)
                                                                     aaGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGc                            >  2:4705847/1‑101 (MQ=255)
                                                                     aaGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGc                            >  1:4705846/1‑101 (MQ=255)
                                                                      aGCTGGAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCg                           <  1:4705350/101‑1 (MQ=255)
                                                                          ggAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATgg                       >  2:4705852/1‑101 (MQ=255)
                                                                           gAGGCACCACCGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGa                      >  1:4705853/1‑101 (MQ=255)
                                                                                accaccGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                 >  2:4705857/1‑101 (MQ=255)
                                                                                 ccaccGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgagga                               <  2:4705738/86‑1 (MQ=255)
                                                                                  caccGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGa                      >  1:4705863/1‑94 (MQ=255)
                                                                                  caccGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa               >  2:4705860/1‑101 (MQ=255)
                                                                                   accGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAc              <  1:4705215/101‑1 (MQ=255)
                                                                                     cGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc            <  1:4705243/101‑1 (MQ=255)
                                                                                     cGGCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc            >  2:4705867/1‑101 (MQ=255)
                                                                                      ggCACGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgaggag                              >  1:4705868/1‑82 (MQ=255)
                                                                                         aCGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggcg        <  1:4705173/101‑1 (MQ=255)
                                                                                          cGGGCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGt       >  1:4705873/1‑101 (MQ=255)
                                                                                           gggCGCGACCTTCTACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACTGACGCGGCGTc      >  2:4705875/1‑101 (MQ=255)
                                                                                           gggCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTc      <  1:4705490/101‑1 (MQ=255)
                                                                                            ggCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCa     <  2:4705268/101‑1 (MQ=255)
                                                                                            ggCGCGACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCa     <  2:4705303/101‑1 (MQ=255)
                                                                                              cgcgACCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGc   <  1:4705570/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCG  >  NC_000913/2890559‑2890754

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

GATK/CNVnator alignment

N/A