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 1,788,045 A→G S408S (AGT→AGC b1706 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,0450AG96.4% 337.4 / 1.3 112S408S (AGT→AGCb1706hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (60/48);  total (61/51)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 119 total aligned reads displayed.

ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGC  >  NC_000913/1787955‑1788135
                                                                                          |                                                                                          
aCCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGAcgctcgct                                                                                          >  2:239865/1‑93 (MQ=255)
 cccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTa                                                                                         >  1:2465511/1‑93 (MQ=255)
 cccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTa                                                                                         <  1:4109750/93‑1 (MQ=255)
  ccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTAt                                                                                        <  2:50648/93‑1 (MQ=255)
  ccTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTAt                                                                                        <  1:2065187/93‑1 (MQ=255)
    ttttCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCg                                                                                      >  2:3067647/1‑93 (MQ=255)
     tttCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGc                                                                                     >  1:346341/1‑93 (MQ=255)
       tCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCtcgctatcgcta                                                                                   <  2:1100768/93‑1 (MQ=255)
        cTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATGCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTaa                                                                                  <  1:1989368/93‑1 (MQ=255)
         tGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAc                                                                                 >  1:3163011/1‑93 (MQ=255)
          gCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAcc                                                                                <  2:3340778/93‑1 (MQ=255)
          gCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAcc                                                                                <  2:3478287/93‑1 (MQ=255)
              cTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGt                                                                            >  1:3157780/1‑93 (MQ=255)
              cTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGt                                                                            <  2:3044657/93‑1 (MQ=255)
               ttCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTc                                                                           >  2:66019/1‑93 (MQ=255)
               ttCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTc                                                                           >  1:3806556/1‑93 (MQ=255)
                tCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCt                                                                          <  2:587095/93‑1 (MQ=255)
                 cAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCtt                                                                         <  1:2613252/93‑1 (MQ=255)
                 cAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCtt                                                                         >  1:82417/1‑93 (MQ=255)
                 cAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCCGACGCTCGCTATCGCTAACCTCGTCtt                                                                         >  1:2401013/1‑93 (MQ=255)
                  aaTCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttg                                                                        <  1:1208104/93‑1 (MQ=255)
                   aTCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgt                                                                       >  1:713143/1‑93 (MQ=255)
                    tCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgtt                                                                      >  1:2281715/1‑93 (MQ=255)
                    tCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgtt                                                                      <  2:39201/93‑1 (MQ=255)
                     cGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCttgttg                                                                     <  1:4652604/93‑1 (MQ=255)
                      gCCCGTTGCGCCAACCAGGTGCGCAATACCAGAGCGGGGCTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGc                                                                    >  2:3648182/1‑93 (MQ=255)
                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  <  2:4638748/93‑1 (MQ=255)
                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  >  1:1650108/1‑93 (MQ=255)
                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  >  1:3945428/1‑93 (MQ=255)
                        ccGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaa                                                                  >  2:2102628/1‑93 (MQ=255)
                         cGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaaa                                                                 <  1:4686356/93‑1 (MQ=255)
                         cGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCaaa                                                                 <  1:4686357/93‑1 (MQ=255)
                         cGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCaaa                                                                 >  2:3987089/1‑93 (MQ=255)
                          gtagcgcCAACCAGTTGCGCAATTCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCTTCTTGTTGCAAAc                                                                <  2:4327266/90‑1 (MQ=255)
                          gTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAAc                                                                <  1:1326850/93‑1 (MQ=255)
                             gcgcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcc                                                             <  1:4035146/93‑1 (MQ=255)
                              cgcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccc                                                            <  1:3008649/93‑1 (MQ=255)
                              cgcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCACTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccc                                                            >  2:1124868/1‑93 (MQ=255)
                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           <  2:2178404/93‑1 (MQ=255)
                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           <  1:2519809/93‑1 (MQ=255)
                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           >  2:2141349/1‑93 (MQ=255)
                               gcCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGcccc                                                           >  1:2957916/1‑93 (MQ=255)
                                ccAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccccc                                                          >  1:3539701/1‑93 (MQ=255)
                                ccAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGccccc                                                          <  1:2877459/93‑1 (MQ=255)
                                  aaCCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGa                                                        >  1:2309654/1‑93 (MQ=255)
                                  aaCCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGa                                                        >  2:2656894/1‑93 (MQ=255)
                                    ccAGTTGCGCAATCCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTCTCGCTAACCTCGTCTTGTTGCAAACGCCCCCCATa                                                      >  1:810546/1‑93 (MQ=255)
                                    ccAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATa                                                      >  1:810547/1‑93 (MQ=255)
                                    ccAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATa                                                      >  2:1859728/1‑93 (MQ=255)
                                    ccAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATa                                                      >  2:1859725/1‑93 (MQ=255)
                                       gTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACg                                                   <  1:2737013/93‑1 (MQ=255)
                                            gcAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCaa                                              >  1:3928424/1‑93 (MQ=255)
                                            gcAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCaa                                              >  2:3048490/1‑93 (MQ=255)
                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  2:2357736/1‑93 (MQ=255)
                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  1:4569313/1‑93 (MQ=255)
                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  1:1390827/1‑93 (MQ=255)
                                              aaTACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAc                                            >  1:1768913/1‑93 (MQ=255)
                                                tACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCa                                          >  2:1104708/1‑93 (MQ=255)
                                                tACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCa                                          >  2:3854426/1‑93 (MQ=255)
                                                 aCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAg                                         >  2:690945/1‑93 (MQ=255)
                                                   cAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtc                                       >  1:1958130/1‑93 (MQ=255)
                                                    agagCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtca                                      <  2:1293119/93‑1 (MQ=255)
                                                     gagCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtcat                                     >  1:2821823/1‑93 (MQ=255)
                                                     gagCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtcat                                     <  2:654343/93‑1 (MQ=255)
                                                       gCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGtcatca                                   >  2:4609622/1‑93 (MQ=255)
                                                         tgggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCaaa                                 >  2:3617225/2‑93 (MQ=255)
                                                         ggggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCaaa                                 <  1:3897064/93‑1 (MQ=255)
                                                          gtgTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAAt                                <  2:4285778/91‑1 (MQ=255)
                                                          gggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAAt                                <  1:1697328/93‑1 (MQ=255)
                                                           ggTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATg                               <  2:1364536/93‑1 (MQ=255)
                                                            gTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGc                              <  1:4605302/93‑1 (MQ=255)
                                                             ttAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGcc                             >  2:4735560/1‑93 (MQ=255)
                                                             ttAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGcc                             >  1:1931638/1‑93 (MQ=255)
                                                                 cGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAc                         <  2:526407/93‑1 (MQ=255)
                                                                  gCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACg                        >  2:2346575/1‑93 (MQ=255)
                                                                   cTTTGCATCAGTTGCTGGCGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGa                       <  2:220286/93‑1 (MQ=255)
                                                                   cTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGa                       <  2:1344109/93‑1 (MQ=255)
                                                                   cTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGa                       <  1:2973917/93‑1 (MQ=255)
                                                                     ttGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATc                     <  2:2306188/93‑1 (MQ=255)
                                                                      tGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCa                    >  1:1679229/1‑93 (MQ=255)
                                                                      tGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCa                    >  1:2709964/1‑93 (MQ=255)
                                                                        cATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAAt                  <  2:1239960/93‑1 (MQ=255)
                                                                         aTCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATa                 <  1:2455028/93‑1 (MQ=255)
                                                                          tCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATaa                >  1:2990263/1‑93 (MQ=255)
                                                                          tCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATaa                >  2:2118235/1‑93 (MQ=255)
                                                                           cAGTTGCAGACGCTCGCTATCGCTAACCTCGTCGTGTTGCAAACGCCTCCGATAACGGGCAACCCAGTCATCAAATGCCGCACGATCAATaaa               <  2:1371/93‑1 (MQ=255)
                                                                            aGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAAc              >  1:3525419/1‑93 (MQ=255)
                                                                            aGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAAc              <  2:673471/93‑1 (MQ=255)
                                                                              ttGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtc            <  1:1026765/93‑1 (MQ=255)
                                                                              ttGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtc            >  2:4471426/1‑93 (MQ=255)
                                                                              ttGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtc            >  2:2917223/1‑93 (MQ=255)
                                                                                gCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCAGCACGATCAATAAACtcat          >  1:1436927/1‑93 (MQ=255)
                                                                                  tGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATCAACtcataa        >  2:1858437/1‑91 (MQ=255)
                                                                                  tGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACtcatca        <  1:1307553/93‑1 (MQ=255)
                                                                                   gACGCTCGCTCTCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCAc       >  2:4347789/1‑93 (MQ=255)
                                                                                     cgctcgctATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGt     >  1:208304/1‑93 (MQ=255)
                                                                                     cgctcgctATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGt     <  1:2793508/93‑1 (MQ=255)
                                                                                      gCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTa    <  2:3896453/93‑1 (MQ=255)
                                                                                      ccTCGCTATCGCTAACCTCGTCTTTTTGCAAACGCCCCCGATAACGGGCCAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTa    <  1:1095027/92‑1 (MQ=255)
                                                                                        tCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGc  <  1:649534/93‑1 (MQ=255)
                                                                                          |                                                                                          
ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGC  >  NC_000913/1787955‑1788135

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
ATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGCGGT  >  NC_000913/1787952‑1788138
                                                                                             |                                                                                             
ATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCG                                                                                               >  SRR4302544.3642396/1‑94 (MQ=60)
  CACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCT                                                                                             <  SRR4302544.3409808/94‑1 (MQ=60)
  CACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCT                                                                                             <  SRR4302544.3915644/94‑1 (MQ=60)
   ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTA                                                                                            >  SRR4302544.2489058/1‑94 (MQ=60)
    CCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTAT                                                                                           >  SRR4302544.510850/1‑94 (MQ=60)
    CCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTAT                                                                                           <  SRR4302544.4150601/94‑1 (MQ=60)
     CCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATC                                                                                          <  SRR4302544.2085033/94‑1 (MQ=60)
     CCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATC                                                                                          <  SRR4302544.3692649/94‑1 (MQ=60)
       TTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGC                                                                                        >  SRR4302544.349905/1‑94 (MQ=60)
           CTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATGCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAC                                                                                    <  SRR4302544.2008681/94‑1 (MQ=60)
           CTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAAC                                                                                    >  SRR4302544.3193098/1‑94 (MQ=60)
                GCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGT                                                                               >  SRR4302544.3187816/1‑94 (MQ=60)
                 CTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTC                                                                              >  SRR4302544.3844277/1‑94 (MQ=60)
                   TCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTT                                                                            >  SRR4302544.83209/1‑94 (MQ=60)
                   TCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCCGACGCTCGCTATCGCTAACCTCGTCTT                                                                            >  SRR4302544.2423794/1‑94 (MQ=60)
                    CAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTG                                                                           <  SRR4302544.2638343/94‑1 (MQ=60)
                     AATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGT                                                                          <  SRR4302544.1220070/94‑1 (MQ=60)
                     AATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGT                                                                          >  SRR4302544.720206/1‑94 (MQ=60)
                      ATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTT                                                                         >  SRR4302544.2303213/1‑94 (MQ=60)
                        CGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGC                                                                       <  SRR4302544.4698649/94‑1 (MQ=60)
                          CCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAA                                                                     >  SRR4302544.1666386/1‑94 (MQ=60)
                          CCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAA                                                                     >  SRR4302544.3984631/1‑94 (MQ=60)
                          CCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAA                                                                     <  SRR4302544.3172810/94‑1 (MQ=60)
                            CGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAAC                                                                   <  SRR4302544.4732679/94‑1 (MQ=60)
                            CGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAAC                                                                   <  SRR4302544.4732678/94‑1 (MQ=60)
                            CGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAAC                                                                   <  SRR4302544.2151079/94‑1 (MQ=60)
                             GTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACG                                                                  <  SRR4302544.1339897/94‑1 (MQ=60)
                                GCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCC                                                               <  SRR4302544.4075337/94‑1 (MQ=60)
                                 CGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCC                                                              >  SRR4302544.2883322/1‑94 (MQ=60)
                                 CGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCC                                                              >  SRR4302544.2986316/1‑94 (MQ=60)
                                 CGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCC                                                              <  SRR4302544.3037486/94‑1 (MQ=60)
                                  GCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCC                                                             <  SRR4302544.2544017/94‑1 (MQ=60)
                                  GCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCC                                                             >  SRR4302544.3573571/1‑94 (MQ=60)
                                  GCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCC                                                             <  SRR4302544.1846871/94‑1 (MQ=60)
                                   CCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCG                                                            <  SRR4302544.2905230/94‑1 (MQ=60)
                                    CAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGA                                                           >  SRR4302544.2331489/1‑94 (MQ=60)
                                      ACCAGTTGCGCAATCCCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTCTCGCTAACCTCGTCTTGTTGCAAACGCCCCCCATA                                                         >  SRR4302544.818449/1‑94 (MQ=60)
                                      ACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATA                                                         >  SRR4302544.818450/1‑94 (MQ=60)
                                         AGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACG                                                      <  SRR4302544.256391/94‑1 (MQ=60)
                                          GTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGG                                                     <  SRR4302544.2763398/94‑1 (MQ=60)
                                              CGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAA                                                 >  SRR4302544.3967438/1‑94 (MQ=60)
                                                CAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGCTAACGGGCAAAC                                               >  SRR4302544.1404710/1‑94 (MQ=60)
                                                CAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAC                                               >  SRR4302544.1786228/1‑94 (MQ=60)
                                                CAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAC                                               >  SRR4302544.4614861/1‑94 (MQ=60)
                                                CAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAAC                                               >  SRR4302544.1404709/1‑94 (MQ=60)
                                                    ACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGT                                           <  SRR4302544.856132/94‑1 (MQ=60)
                                                     CCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTC                                          >  SRR4302544.1977125/1‑94 (MQ=60)
                                                       AGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCAT                                        >  SRR4302544.2849082/1‑94 (MQ=60)
                                                            GGGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAAT                                   <  SRR4302544.3935768/94‑1 (MQ=60)
                                                             GGGTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATG                                  <  SRR4302544.1714076/94‑1 (MQ=60)
                                                               GTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCC                                >  SRR4302544.1950330/1‑94 (MQ=60)
                                                               GTTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCC                                <  SRR4302544.4651062/94‑1 (MQ=60)
                                                                TTAACGCTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCG                               <  SRR4302544.4493201/94‑1 (MQ=60)
                                                                      CTTTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGAT                         <  SRR4302544.3002443/94‑1 (MQ=60)
                                                                        TTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCA                       >  SRR4302544.1695765/1‑94 (MQ=60)
                                                                        TTGCATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCA                       >  SRR4302544.2736027/1‑94 (MQ=60)
                                                                            ATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAA                   >  SRR4302544.3018940/1‑94 (MQ=60)
                                                                            ATCAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAA                   <  SRR4302544.2478448/94‑1 (MQ=60)
                                                                              CAGTTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAAC                 >  SRR4302544.3559107/1‑94 (MQ=60)
                                                                                 TTGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCA              <  SRR4302544.1037007/94‑1 (MQ=60)
                                                                                  TGCTGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCAGCACGATCAATAAACTCAT             >  SRR4302544.1451236/1‑94 (MQ=60)
                                                                                     TGACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCAC          <  SRR4302544.1320411/94‑1 (MQ=60)
                                                                                       ACGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGT        >  SRR4302544.210442/1‑94 (MQ=60)
                                                                                        CGCTCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTA       <  SRR4302544.2820407/94‑1 (MQ=60)
                                                                                           TCGCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGCG    <  SRR4302544.3057243/94‑1 (MQ=60)
                                                                                           TCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGCG    <  SRR4302544.655955/94‑1 (MQ=60)
                                                                                             GCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGCGGT  >  SRR4302544.3625273/1‑94 (MQ=60)
                                                                                             |                                                                                             
ATCACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCACGATCAATAAACTCATCACGTAGCGGT  >  NC_000913/1787952‑1788138

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