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
A1 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,391,348 A→G A573A (GCT→GCC nuoL ← NADH dehydrogenase I chain L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,3480AG97.6% 380.6 / ‑1.0 125A573A (GCT→GCCnuoLNADH dehydrogenase I chain L
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (67/55);  total (68/57)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 129 total aligned reads displayed.

CGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAG  >  NC_000913/2391257‑2391437
                                                                                           |                                                                                         
cGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAg                                                                                          <  2:2154844/93‑1 (MQ=255)
 gACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCATACCTTTACCTGCAAAGCGGGAAAGGACGGc                                                                                         >  1:3568082/1‑93 (MQ=255)
 gACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGc                                                                                         >  2:633142/1‑93 (MQ=255)
 gACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGc                                                                                         <  2:4326121/93‑1 (MQ=255)
 gACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGc                                                                                         <  2:106220/93‑1 (MQ=255)
  aCCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGcc                                                                                        >  1:4134191/1‑93 (MQ=255)
   ccGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCg                                                                                       <  1:953348/93‑1 (MQ=255)
    cGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCgg                                                                                      >  2:412552/1‑93 (MQ=255)
     gCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCggg                                                                                     >  1:4442828/1‑93 (MQ=255)
      gtccGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCTCTTAACAGCAGAGCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGa                                                                                    <  2:3183837/91‑1 (MQ=255)
       aCCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGAt                                                                                   <  1:3110187/93‑1 (MQ=255)
       aCCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGAt                                                                                   <  1:2166450/93‑1 (MQ=255)
        ccGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATg                                                                                  <  2:3884681/93‑1 (MQ=255)
         cGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGt                                                                                 >  2:3015529/1‑93 (MQ=255)
         cGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGt                                                                                 >  2:3015530/1‑93 (MQ=255)
           aTGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGCtc                                                                               >  2:1991043/1‑93 (MQ=255)
            tGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTtca                                                                              <  2:3098554/93‑1 (MQ=255)
            ggCTCATAGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTtca                                                                              >  2:2882108/2‑93 (MQ=255)
              cTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTtcatc                                                                            <  2:4494437/93‑1 (MQ=255)
                cATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTcatcat                                                                          >  2:3476029/1‑93 (MQ=255)
                    gATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAg                                                                      <  2:2415033/93‑1 (MQ=255)
                         cacaTACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAg                                                                 >  1:1215632/1‑93 (MQ=255)
                         cacaTACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAg                                                                 <  2:4044056/93‑1 (MQ=255)
                          acaTACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGc                                                                >  1:1266867/1‑93 (MQ=255)
                             tACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGa                                                             >  2:117429/1‑93 (MQ=255)
                               ccAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATc                                                           <  1:2517690/93‑1 (MQ=255)
                                cAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATcg                                                          >  2:280145/1‑93 (MQ=255)
                                 aGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATcgc                                                         <  1:3613689/93‑1 (MQ=255)
                                    gcAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCACAGGATCGCAtt                                                      >  1:4108202/1‑93 (MQ=255)
                                      aGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTc                                                    >  1:650955/1‑93 (MQ=255)
                                      aGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTc                                                    >  2:4608426/1‑93 (MQ=255)
                                       gATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCa                                                   >  1:1476938/1‑93 (MQ=255)
                                       gATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCa                                                   >  1:3925000/1‑93 (MQ=255)
                                       gATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCa                                                   <  1:3948636/93‑1 (MQ=255)
                                       gATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCa                                                   <  2:771072/93‑1 (MQ=255)
                                         tAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGc                                                 >  1:1719780/1‑93 (MQ=255)
                                         aaGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGc                                                 <  1:2739361/92‑1 (MQ=255)
                                            ccGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAAc                                              <  2:802205/93‑1 (MQ=255)
                                               ttCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAg                                           >  2:3856901/1‑93 (MQ=255)
                                                tctcACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAgg                                          >  1:1471852/1‑93 (MQ=255)
                                                 ctcACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGc                                         >  2:636704/1‑93 (MQ=255)
                                                   cACTTAACAGCAGGCCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCaa                                       <  2:3848583/93‑1 (MQ=255)
                                                   cACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCaa                                       <  1:4161408/93‑1 (MQ=255)
                                                   cACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCaa                                       <  1:1029426/93‑1 (MQ=255)
                                                    aCTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAAt                                      <  1:3296403/93‑1 (MQ=255)
                                                    aCTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAAt                                      <  2:1121573/93‑1 (MQ=255)
                                                     cTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATa                                     <  1:2744964/93‑1 (MQ=255)
                                                     cTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATa                                     <  2:1266867/93‑1 (MQ=255)
                                                     cTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATa                                     <  2:882160/93‑1 (MQ=255)
                                                      ttAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATAc                                    >  1:2634549/1‑93 (MQ=255)
                                                      ttAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATAc                                    >  2:1422856/1‑93 (MQ=255)
                                                      ttAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGAGGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACACGCCAATAc                                    >  1:2634547/1‑93 (MQ=255)
                                                        aaCAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATAccc                                  <  1:295764/93‑1 (MQ=255)
                                                        aaCAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATAccc                                  <  2:2200369/93‑1 (MQ=255)
                                                         aCAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCa                                 <  2:330473/93‑1 (MQ=255)
                                                         aCAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCa                                 >  2:1403669/1‑93 (MQ=255)
                                                          cagcagACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAg                                >  1:1149864/1‑93 (MQ=255)
                                                          cagcagACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAg                                <  1:1925026/93‑1 (MQ=255)
                                                           ataagACCTTTACCAGCCAAACGAGATAGGACAGCCGGGATGTTCATCATTGAGCTCAGCGGATCGCGTTCAACCAACAAGGCAATACCCAgg                               >  2:3331541/4‑93 (MQ=255)
                                                            gcagACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGa                              >  1:597291/1‑93 (MQ=255)
                                                             cagACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGaa                             <  1:4358080/93‑1 (MQ=255)
                                                               gACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACg                           >  1:1788638/1‑93 (MQ=255)
                                                               gACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACg                           >  1:2060372/1‑93 (MQ=255)
                                                               gACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACg                           >  1:3539775/1‑93 (MQ=255)
                                                               gACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATATTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACg                           >  2:4149846/1‑93 (MQ=255)
                                                                aCCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACgg                          <  1:3848063/93‑1 (MQ=255)
                                                                aCCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACgg                          >  1:4632728/1‑93 (MQ=255)
                                                                aCCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACgg                          >  1:2887155/1‑93 (MQ=255)
                                                                aCCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACgg                          <  2:2203913/93‑1 (MQ=255)
                                                                aCCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACgg                          <  2:889802/93‑1 (MQ=255)
                                                                 ccTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGc                         <  1:2622161/93‑1 (MQ=255)
                                                                  cTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATGCCCAGGAACGGCt                        <  2:2224038/93‑1 (MQ=255)
                                                                  cTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCt                        <  2:3956234/93‑1 (MQ=255)
                                                                   tttACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCtt                       >  2:2115229/1‑93 (MQ=255)
                                                                   tttACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCtt                       >  1:855845/1‑93 (MQ=255)
                                                                   tttACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCtt                       >  2:2300295/1‑93 (MQ=255)
                                                                   tttACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCtt                       >  2:1639521/1‑93 (MQ=255)
                                                                    ttACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTg                      <  1:1714005/93‑1 (MQ=255)
                                                                      aCCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGAc                    >  2:3310072/1‑93 (MQ=255)
                                                                      aCCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGAc                    >  1:3298912/1‑93 (MQ=255)
                                                                       ccTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACg                   >  1:4240877/1‑93 (MQ=255)
                                                                       ccTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTAAGCAACCAGGCAATACCCACGAACGGCTTGACg                   >  1:2739059/1‑93 (MQ=255)
                                                                        cTGCAAAGCGGGAAAGGACGGCCGGGATTTTCATCATTGAGTTCAGCGGATCGAGTTTCAGCAACCAGGCAATACCCAGCAACGGCTTgcaca                  >  1:4404681/1‑89 (MQ=255)
                                                                        cTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGa                  >  2:2654762/1‑93 (MQ=255)
                                                                         tGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGaa                 <  2:3432762/93‑1 (MQ=255)
                                                                            aaaGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAacac              <  1:70325/93‑1 (MQ=255)
                                                                              aGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACtt            >  2:531565/1‑93 (MQ=255)
                                                                              aGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCGTGACGAACACtt            >  1:3351348/1‑93 (MQ=255)
                                                                                 gggAAAGTACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGt         <  2:3374841/93‑1 (MQ=255)
                                                                                 gggAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGt         <  2:358295/93‑1 (MQ=255)
                                                                                  ggAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTc        >  2:60670/1‑93 (MQ=255)
                                                                                   gAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCa       >  1:4234783/1‑93 (MQ=255)
                                                                                   gAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCa       <  2:2177346/93‑1 (MQ=255)
                                                                                    aaaGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAt      >  1:3175245/1‑93 (MQ=255)
                                                                                      aGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATAc    >  2:4366644/1‑93 (MQ=255)
                                                                                      aGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATAc    >  2:672777/1‑93 (MQ=255)
                                                                                      aGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATAc    <  2:2009624/93‑1 (MQ=255)
                                                                                       ggACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACa   <  1:2349107/93‑1 (MQ=255)
                                                                                        gACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAg  >  1:3572651/1‑93 (MQ=255)
                                                                                        gACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAg  <  1:1206083/93‑1 (MQ=255)
                                                                                           |                                                                                         
CGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAG  >  NC_000913/2391257‑2391437

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGCCA  >  NC_000913/2391256‑2391440
                                                                                            |                                                                                            
ACGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGG                                                                                             <  SRR4302544.1590593/94‑1 (MQ=60)
 CGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCATACCTTTACCTGCAAAGCGGGAAAGGACGGC                                                                                            >  SRR4302544.3602294/1‑94 (MQ=60)
  GACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCC                                                                                           >  SRR4302544.4175151/1‑94 (MQ=60)
    CCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGG                                                                                         <  SRR4302544.962720/94‑1 (MQ=60)
     CGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGG                                                                                        >  SRR4302544.4487612/1‑94 (MQ=60)
        ACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATG                                                                                     <  SRR4302544.3139782/94‑1 (MQ=60)
        ACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATG                                                                                     <  SRR4302544.2187012/94‑1 (MQ=60)
                      ATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTT                                                                       >  SRR4302544.4220057/1‑94 (MQ=60)
                         CCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAG                                                                    >  SRR4302544.1227665/1‑94 (MQ=60)
                          CACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGC                                                                   >  SRR4302544.1279322/1‑94 (MQ=60)
                                CCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCG                                                             <  SRR4302544.2541868/94‑1 (MQ=60)
                                  AGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCG                                                           <  SRR4302544.3649180/94‑1 (MQ=60)
                                    CGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCACAGGATCGCATT                                                         >  SRR4302544.4149046/1‑94 (MQ=60)
                                      CAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTC                                                       >  SRR4302544.657396/1‑94 (MQ=60)
                                      CAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTC                                                       >  SRR4302544.3787014/1‑94 (MQ=60)
                                       AGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCA                                                      >  SRR4302544.1491545/1‑94 (MQ=60)
                                       AGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCA                                                      >  SRR4302544.3963988/1‑94 (MQ=60)
                                        GATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAG                                                     <  SRR4302544.3987878/94‑1 (MQ=60)
                                         ATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGC                                                    >  SRR4302544.1736679/1‑94 (MQ=60)
                                          AAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCA                                                   <  SRR4302544.2765770/94‑1 (MQ=60)
                                              CGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCA                                               >  SRR4302544.266381/1‑94 (MQ=60)
                                               ctgcgtcctgacccacagcccttgacgtgcaacgcgtGAAAGGACGGCCGGGATGTTCATCATTGAGGTCAGCGGATCGCGTTTCAGCAACCAG                                              <  SRR4302544.1999155/57‑1 (MQ=60)
                                                TTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGG                                             >  SRR4302544.1486427/1‑94 (MQ=60)
                                                  CTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCA                                           >  SRR4302544.3544501/1‑94 (MQ=60)
                                                   TCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTGCAGCGGATCGCGTTTCAGCAACCAGGCAA                                          >  SRR4302544.647323/1‑94 (MQ=60)
                                                    CACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAAT                                         <  SRR4302544.4202630/94‑1 (MQ=60)
                                                    CACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAAT                                         <  SRR4302544.1039686/94‑1 (MQ=60)
                                                     ACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATA                                        <  SRR4302544.3327797/94‑1 (MQ=60)
                                                      CTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATAC                                       >  SRR4302544.2659799/1‑94 (MQ=60)
                                                      CTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATAC                                       <  SRR4302544.2771413/94‑1 (MQ=60)
                                                      CTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGAGGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACacgccaatac                                       >  SRR4302544.2659797/1‑84 (MQ=60)
                                                         AACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCA                                    <  SRR4302544.298822/94‑1 (MQ=60)
                                                          ACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAG                                   >  SRR4302544.1161302/1‑94 (MQ=60)
                                                           CAGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGG                                  <  SRR4302544.1943650/94‑1 (MQ=60)
                                                            AGCAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGA                                 >  SRR4302544.603145/1‑94 (MQ=60)
                                                              CAGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAAC                               <  SRR4302544.4402317/94‑1 (MQ=60)
                                                               AGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACG                              >  SRR4302544.3573645/1‑94 (MQ=60)
                                                               AGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACG                              >  SRR4302544.1806121/1‑94 (MQ=60)
                                                               AGACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACG                              >  SRR4302544.2080187/1‑94 (MQ=60)
                                                                GACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGG                             >  SRR4302544.2914998/1‑94 (MQ=60)
                                                                GACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGG                             >  SRR4302544.4678646/1‑94 (MQ=60)
                                                                GACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGG                             >  SRR4302544.2914997/1‑94 (MQ=60)
                                                                 ACCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGC                            <  SRR4302544.3886251/94‑1 (MQ=60)
                                                                  CCTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCT                           <  SRR4302544.2647308/94‑1 (MQ=60)
                                                                   CTTTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTT                          >  SRR4302544.864312/1‑94 (MQ=60)
                                                                     TTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGA                        <  SRR4302544.3227647/94‑1 (MQ=60)
                                                                     TTACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGA                        <  SRR4302544.1730864/94‑1 (MQ=60)
                                                                      TACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGAC                       <  SRR4302544.1951834/94‑1 (MQ=60)
                                                                      TACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGAC                       <  SRR4302544.4451335/94‑1 (MQ=60)
                                                                      TACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGAC                       >  SRR4302544.3330335/1‑94 (MQ=60)
                                                                       ACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACG                      >  SRR4302544.4284099/1‑94 (MQ=60)
                                                                       ACCTGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTAAGCAACCAGGCAATACCCACGAACGGCTTGACG                      >  SRR4302544.2765464/1‑94 (MQ=60)
                                                                        CCTGCAAAGCGGGAAAGGACGGCCGGGATTTTCATCATTGAGTTCAGCGGATCGAGTTTCAGCAACCAGGCAATACCCAGCAACGGCTTGcaca                     >  SRR4302544.4449171/1‑90 (MQ=60)
                                                                          TGCAAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAAC                   <  SRR4302544.1254050/94‑1 (MQ=60)
                                                                             AAAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACT                <  SRR4302544.71011/94‑1 (MQ=60)
                                                                              AAGCGGGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCGTGACGAACACTT               >  SRR4302544.3383328/1‑94 (MQ=60)
                                                                                   GGAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCA          >  SRR4302544.4277947/1‑94 (MQ=60)
                                                                                    GAAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAT         >  SRR4302544.3205460/1‑94 (MQ=60)
                                                                                     AAAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATA        >  SRR4302544.2449459/1‑94 (MQ=60)
                                                                                      AAGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATAC       >  SRR4302544.1551110/1‑94 (MQ=60)
                                                                                       AGGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACA      <  SRR4302544.4029446/94‑1 (MQ=60)
                                                                                        GGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAG     >  SRR4302544.3606909/1‑94 (MQ=60)
                                                                                        GGACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAG     <  SRR4302544.2371359/94‑1 (MQ=60)
                                                                                         GACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGC    <  SRR4302544.1218028/94‑1 (MQ=60)
                                                                                          ACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGCC   >  SRR4302544.3516537/1‑94 (MQ=60)
                                                                                          ACGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGCC   <  SRR4302544.1784553/94‑1 (MQ=60)
                                                                                           CGGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGCTTCAGCAACCAGGCAATACCCAGGAACGGCTTTACGAACACCTTGTCATACAGCCT  >  SRR4302544.3393305/1‑94 (MQ=60)
                                                                                            |                                                                                            
ACGACCGCACCGATGCTCATGGATGCCACATACCAGCGCAGATAGCCGTTCTCACTTAACAGCAGACCTTTACCTGCAAAGCGGGAAAGGACAGCCGGGATGTTCATCATTGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGCCA  >  NC_000913/2391256‑2391440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: