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,025,930 G→A S42F (TCC→TTC)  b1956 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,025,9300GA96.0% 423.3 / 2.1 124S42F (TCC→TTC) b1956hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (58/61);  total (60/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 129 total aligned reads displayed.

GAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  >  NC_000913/2025839‑2026019
                                                                                           |                                                                                         
gAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGa                                                                                          <  2:3108074/93‑1 (MQ=255)
 aGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaa                                                                                         >  2:4339448/1‑93 (MQ=255)
  gcgcGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaaa                                                                                        >  2:2224674/1‑93 (MQ=255)
  gcgcGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaaa                                                                                        >  1:2752161/1‑93 (MQ=255)
  gcgcGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaaa                                                                                        <  2:1377684/93‑1 (MQ=255)
   cgcgTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAaaaa                                                                                       <  1:2872081/93‑1 (MQ=255)
      gTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAg                                                                                    <  1:3247685/93‑1 (MQ=255)
      gTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAg                                                                                    >  2:3187912/1‑93 (MQ=255)
       ttGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGc                                                                                   >  2:441521/1‑93 (MQ=255)
        tGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCa                                                                                  <  2:378032/93‑1 (MQ=255)
        tGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCa                                                                                  >  1:1693464/1‑93 (MQ=255)
        tGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCa                                                                                  >  1:266637/1‑93 (MQ=255)
          gCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAAt                                                                                <  1:1491794/93‑1 (MQ=255)
           cAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATa                                                                               <  2:3680147/93‑1 (MQ=255)
           cAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATa                                                                               <  1:510987/93‑1 (MQ=255)
           cAACGTTGTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATa                                                                               <  1:130797/93‑1 (MQ=255)
            aaCGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGCGAGCAAGGTGAAAAAAAGCAATAc                                                                              >  1:3265796/1‑93 (MQ=255)
             aCGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATAcc                                                                             <  1:2700816/93‑1 (MQ=255)
               gTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCAc                                                                           >  2:4077860/1‑93 (MQ=255)
                ttttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACg                                                                          >  1:550444/1‑93 (MQ=255)
                 tttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGa                                                                         <  2:2462406/93‑1 (MQ=255)
                 tttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGa                                                                         >  2:1422045/1‑93 (MQ=255)
                 tttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGa                                                                         <  2:375198/93‑1 (MQ=255)
                     cAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATaaa                                                                     <  1:4072125/93‑1 (MQ=255)
                     cAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATaaa                                                                     <  1:2179950/93‑1 (MQ=255)
                       gatgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGc                                                                   >  1:248377/1‑93 (MQ=255)
                       gatgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGc                                                                   <  2:408919/93‑1 (MQ=255)
                        atgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCa                                                                  <  2:451345/93‑1 (MQ=255)
                         tgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCag                                                                 >  1:462031/1‑93 (MQ=255)
                         tgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCag                                                                 >  2:501168/1‑93 (MQ=255)
                          gatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCaga                                                                >  1:4534981/1‑93 (MQ=255)
                            ttACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAgaga                                                              <  1:2408933/93‑1 (MQ=255)
                            ttACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGAagagt                                                              >  1:825017/1‑92 (MQ=255)
                             tACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAt                                                             <  1:2667808/93‑1 (MQ=255)
                              aCGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAtt                                                            <  1:1625087/93‑1 (MQ=255)
                              aCGCTGGCTGGAGATATAGGCATCTACCAGCACCACCACTACACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAtt                                                            <  2:2598484/93‑1 (MQ=255)
                                gctggctgGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTaa                                                          <  1:780984/93‑1 (MQ=255)
                                gctggctgGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTaa                                                          >  2:2292128/1‑93 (MQ=255)
                                gctggctgGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTaa                                                          <  2:330014/93‑1 (MQ=255)
                                  tggctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACg                                                        >  2:1281030/1‑93 (MQ=255)
                                   ggctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGa                                                       >  1:2331319/1‑93 (MQ=255)
                                   ggctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGa                                                       <  1:4255370/93‑1 (MQ=255)
                                     ctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACa                                                     >  2:1493970/1‑93 (MQ=255)
                                      tggAGATATAGGAATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACAt                                                    <  1:1161790/93‑1 (MQ=255)
                                        gagaTATAGGCATCTTCCAGCACCACCACTTCGCGCCAGGTGAGCAAGGTGAAAACAAGCAATACCACGATCAAGAGGAGGTTAACGACATGa                                                  >  1:3503264/1‑93 (MQ=255)
                                        gagaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATgc                                                  >  1:3503263/1‑92 (MQ=255)
                                         agaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAc                                                 >  1:3810970/1‑93 (MQ=255)
                                          gaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAcc                                                >  2:2035108/1‑93 (MQ=255)
                                          gaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAcc                                                <  2:2672602/93‑1 (MQ=255)
                                          gaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAcc                                                >  1:1116279/1‑93 (MQ=255)
                                            tataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              <  1:457773/93‑1 (MQ=255)
                                            tataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              >  1:103488/1‑93 (MQ=255)
                                            tataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              <  2:2029995/93‑1 (MQ=255)
                                            tataGGCATCATCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              <  2:2679864/93‑1 (MQ=255)
                                             ataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAgg                                             >  1:4624128/1‑93 (MQ=255)
                                             ataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAgg                                             <  2:735153/93‑1 (MQ=255)
                                                 gCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGcccc                                         <  1:2006781/93‑1 (MQ=255)
                                                  cATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCa                                        >  1:4462362/1‑93 (MQ=255)
                                                  cATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCa                                        <  1:1626487/93‑1 (MQ=255)
                                                  cATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCa                                        <  2:1601466/93‑1 (MQ=255)
                                                   aTCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAg                                       <  2:219334/93‑1 (MQ=255)
                                                   aTCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAg                                       >  2:797474/1‑93 (MQ=255)
                                                    tCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAgg                                      <  2:2257475/93‑1 (MQ=255)
                                                    tCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAgg                                      >  1:1476681/1‑93 (MQ=255)
                                                     cTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggc                                     <  1:3620971/93‑1 (MQ=255)
                                                     cTTCCAGCACCACCACTGCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggc                                     <  2:166895/93‑1 (MQ=255)
                                                      ttCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggcg                                    >  2:3055662/1‑93 (MQ=255)
                                                       tCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggcgg                                   <  2:3585341/93‑1 (MQ=255)
                                                         cAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGgcggcg                                 >  2:444587/1‑93 (MQ=255)
                                                            caccaccacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGc                              <  1:3304649/93‑1 (MQ=255)
                                                               aaccacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCTGAGATTAACGACATGACAAGGCCCCAGGCGGCGTGCGAg                           >  1:313219/2‑93 (MQ=255)
                                                                accacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCATGCGGCGTGCGAGt                          >  2:3092579/1‑93 (MQ=255)
                                                                accacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGt                          >  2:288233/1‑93 (MQ=255)
                                                                  cacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttt                        >  1:4115626/1‑93 (MQ=255)
                                                                  cacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttt                        <  1:776945/93‑1 (MQ=255)
                                                                   acTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttt                       >  2:3637219/1‑93 (MQ=255)
                                                                   acTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttt                       >  1:3857200/1‑93 (MQ=255)
                                                                   acTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttt                       <  1:346209/93‑1 (MQ=255)
                                                                    cTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttttt                      >  1:963119/1‑93 (MQ=255)
                                                                    cTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttttt                      <  2:4413586/93‑1 (MQ=255)
                                                                    cTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttttt                      >  2:123125/1‑93 (MQ=255)
                                                                         cGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAAc                 <  2:4681120/93‑1 (MQ=255)
                                                                           ccAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCa               >  1:4306029/1‑93 (MQ=255)
                                                                           ccAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCa               <  2:1026780/93‑1 (MQ=255)
                                                                             aGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGc             <  2:391759/93‑1 (MQ=255)
                                                                              ggTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGct            >  1:3872789/1‑93 (MQ=255)
                                                                              ggTGAGCAAGGTGAAAAAAAGCAAAACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGct            <  2:3376635/93‑1 (MQ=255)
                                                                               gTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctc           <  2:3934862/93‑1 (MQ=255)
                                                                               gTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctc           >  1:489389/1‑93 (MQ=255)
                                                                                tGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctct          <  1:2877142/93‑1 (MQ=255)
                                                                                tGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctct          >  2:779013/1‑93 (MQ=255)
                                                                                 gAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTg         <  2:4443086/93‑1 (MQ=255)
                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       >  1:4077359/1‑93 (MQ=255)
                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       <  1:926096/93‑1 (MQ=255)
                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       <  2:1009718/93‑1 (MQ=255)
                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       >  2:2380721/1‑93 (MQ=255)
                                                                                    cAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGtt      <  1:1441806/93‑1 (MQ=255)
                                                                                    cAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGtt      >  2:3280056/1‑93 (MQ=255)
                                                                                     aaGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGttt     >  2:1884309/1‑93 (MQ=255)
                                                                                        gTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTcc  >  2:2057618/1‑93 (MQ=255)
                                                                                           |                                                                                         
GAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  >  NC_000913/2025839‑2026019

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
TCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  >  NC_000913/2025837‑2026019
                                                                                             |                                                                                         
TCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGA                                                                                           <  SRR4302544.3206830/94‑1 (MQ=60)
 CGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAA                                                                                          >  SRR4302544.1099576/1‑94 (MQ=60)
   AGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAA                                                                                        >  SRR4302544.2778662/1‑94 (MQ=60)
     CGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAA                                                                                      <  SRR4302544.2899802/94‑1 (MQ=60)
       CGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGTCATCTTCCATCACCACCACTTCCCGCCAGGTGAGCAAGGTGAAAAAAAG                                                                                    >  SRR4302544.1451269/1‑94 (MQ=60)
        GTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGC                                                                                   <  SRR4302544.3278770/94‑1 (MQ=60)
         TTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCA                                                                                  >  SRR4302544.1710174/1‑94 (MQ=60)
         TTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCA                                                                                  >  SRR4302544.269451/1‑94 (MQ=60)
            GCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATA                                                                               <  SRR4302544.1506546/94‑1 (MQ=60)
             CAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATAC                                                                              <  SRR4302544.516028/94‑1 (MQ=60)
             CAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGCGAGCAAGGTGAAAAAAAGCAATAC                                                                              >  SRR4302544.3296972/1‑94 (MQ=60)
             CAACGTTGTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATAC                                                                              <  SRR4302544.131999/94‑1 (MQ=60)
               ACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCA                                                                            <  SRR4302544.2726781/94‑1 (MQ=60)
                 GTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACG                                                                          >  SRR4302544.555837/1‑94 (MQ=60)
                   TTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGAT                                                                        >  SRR4302544.3466344/1‑94 (MQ=60)
                       CAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAG                                                                    <  SRR4302544.2200591/94‑1 (MQ=60)
                       CAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAG                                                                    <  SRR4302544.4112649/94‑1 (MQ=60)
                       CAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAC                                                                    >  SRR4302544.310578/1‑94 (MQ=60)
                        AGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGC                                                                   >  SRR4302544.250981/1‑94 (MQ=60)
                          ATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAG                                                                 >  SRR4302544.466570/1‑94 (MQ=60)
                          ATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAG                                                                 <  SRR4302544.1537334/94‑1 (MQ=60)
                           TGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGA                                                                >  SRR4302544.4580330/1‑94 (MQ=60)
                             ATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGAAGAGT                                                              >  SRR4302544.833102/1‑94 (MQ=60)
                              TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAT                                                             <  SRR4302544.2431819/94‑1 (MQ=60)
                              TTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAT                                                             <  SRR4302544.4702454/94‑1 (MQ=60)
                               TACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATT                                                            <  SRR4302544.2693382/94‑1 (MQ=60)
                                ACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTA                                                           <  SRR4302544.1641099/94‑1 (MQ=60)
                                  GCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAAC                                                         <  SRR4302544.788556/94‑1 (MQ=60)
                                   GTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACG                                                        <  SRR4302544.2442240/94‑1 (MQ=60)
                                    TGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGA                                                       >  SRR4302544.2353426/1‑94 (MQ=60)
                                     GGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGAC                                                      <  SRR4302544.4298751/94‑1 (MQ=60)
                                     GGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGAC                                                      <  SRR4302544.4168774/94‑1 (MQ=60)
                                       CTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACAT                                                    <  SRR4302544.582622/94‑1 (MQ=60)
                                        TGGAGATATAGGAATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATG                                                   <  SRR4302544.1173366/94‑1 (MQ=60)
                                         GGAGATATAGGCATCTTCCAGCACCACCACTTCGCGCCAGGTGAGCAAGGTGAAAACAAGCAATACCACGATCAAGAGGAGGTTAACGACATGA                                                  >  SRR4302544.3536683/1‑94 (MQ=60)
                                         GGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGC                                                  >  SRR4302544.3536682/1‑94 (MQ=60)
                                          GAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAC                                                 >  SRR4302544.3848737/1‑94 (MQ=60)
                                           AGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACC                                                >  SRR4302544.1127422/1‑94 (MQ=60)
                                            GATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCA                                               >  SRR4302544.3568605/1‑94 (MQ=60)
                                             ATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAG                                              >  SRR4302544.104439/1‑94 (MQ=60)
                                              TATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGG                                             >  SRR4302544.4669997/1‑94 (MQ=60)
                                              TATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGG                                             <  SRR4302544.462280/94‑1 (MQ=60)
                                                   GCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCA                                        >  SRR4302544.4507273/1‑94 (MQ=60)
                                                   GCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCA                                        <  SRR4302544.2026279/94‑1 (MQ=60)
                                                    CATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAG                                       <  SRR4302544.4076276/94‑1 (MQ=60)
                                                    CATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAG                                       <  SRR4302544.1642511/94‑1 (MQ=60)
                                                     ATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGG                                      >  SRR4302544.1491288/1‑94 (MQ=60)
                                                       CTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCG                                    <  SRR4302544.3657048/94‑1 (MQ=60)
                                                              CACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCG                             <  SRR4302544.3336128/94‑1 (MQ=60)
                                                                CAACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCTGAGATTAACGACATGACAAGGCCCCAGGCGGCGTGCGAG                           >  SRR4302544.316471/1‑94 (MQ=60)
                                                                   CCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTT                        >  SRR4302544.4156508/1‑94 (MQ=60)
                                                                    CACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTT                       >  SRR4302544.3895490/1‑94 (MQ=60)
                                                                    CACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTT                       <  SRR4302544.784483/94‑1 (MQ=60)
                                                                     ACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTT                      <  SRR4302544.349771/94‑1 (MQ=60)
                                                                     ACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTT                      >  SRR4302544.972577/1‑94 (MQ=60)
                                                                            GCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCA               >  SRR4302544.4349937/1‑94 (MQ=60)
                                                                              CAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGC             >  SRR4302544.2073982/1‑94 (MQ=60)
                                                                               AGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCT            >  SRR4302544.3911250/1‑94 (MQ=60)
                                                                                GGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTC           >  SRR4302544.494190/1‑94 (MQ=60)
                                                                                  TGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTG         <  SRR4302544.2904912/94‑1 (MQ=60)
                                                                                    AGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGT       >  SRR4302544.4117935/1‑94 (MQ=60)
                                                                                     GCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTT      <  SRR4302544.935220/94‑1 (MQ=60)
                                                                                     GCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTT      <  SRR4302544.935221/94‑1 (MQ=60)
                                                                                      CAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTT     <  SRR4302544.1456150/94‑1 (MQ=60)
                                                                                      CAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTT     >  SRR4302544.2761330/1‑94 (MQ=60)
                                                                                         GGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  <  SRR4302544.4534192/94‑1 (MQ=60)
                                                                                             |                                                                                         
TCGAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  >  NC_000913/2025837‑2026019

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