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 I2 R1 593 264.2 13311096 97.8% 13018251 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,621,342 C→T G93G (GGC→GGT)  ppk → polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,621,3420CT87.5% 695.8 / 66.3 239G93G (GGC→GGT) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base T (103/106);  total (121/118)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

Only 100 of 247 total aligned reads displayed.

ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG  >  NC_000913/2621250‑2621434
                                                                                            |                                                                                            
attaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGc                                                                                              >  1:5612695/1‑93 (MQ=255)
  taGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTct                                                                                            <  1:3529133/93‑1 (MQ=255)
   aGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGtctc                                                                                           >  2:2678307/1‑93 (MQ=255)
    gCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAACATTCAGTCCCGGGTGCTGAAAGCCGATCAGGTATTCGACGGTctct                                                                                          >  1:1088832/1‑93 (MQ=255)
    gCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTctct                                                                                          >  2:2820286/1‑93 (MQ=255)
    gCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTctct                                                                                          <  1:4731522/93‑1 (MQ=255)
     cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTa                                                                                         <  1:114724/93‑1 (MQ=255)
     cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTa                                                                                         >  2:1160878/1‑93 (MQ=255)
     cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTa                                                                                         <  1:570987/93‑1 (MQ=255)
      gaagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTAc                                                                                        >  2:2512900/1‑93 (MQ=255)
       aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACa                                                                                       >  2:4684322/1‑93 (MQ=255)
        agaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACaa                                                                                      <  2:3473227/93‑1 (MQ=255)
         gaaCAAGGCTCCAACTCTCATTCCCGCCATGTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAAc                                                                                     <  2:787284/93‑1 (MQ=255)
          aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACg                                                                                    <  2:5064025/93‑1 (MQ=255)
           aCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGa                                                                                   <  1:1968440/93‑1 (MQ=255)
           aCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGa                                                                                   >  1:4488442/1‑93 (MQ=255)
              aGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCt                                                                                >  1:3070862/1‑93 (MQ=255)
              aGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCt                                                                                >  2:2016198/1‑93 (MQ=255)
               ggCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTa                                                                               <  1:4777943/93‑1 (MQ=255)
               agCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTa                                                                               <  2:1427808/92‑1 (MQ=255)
                gCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTAt                                                                              <  1:579255/93‑1 (MQ=255)
                  tCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTg                                                                            >  2:1910257/1‑93 (MQ=255)
                   ccAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGc                                                                           >  1:1119110/1‑93 (MQ=255)
                     aaCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTg                                                                         >  2:3074675/1‑93 (MQ=255)
                      aCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTgg                                                                        >  1:5661545/1‑93 (MQ=255)
                       ctctCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGga                                                                       <  1:556701/93‑1 (MQ=255)
                         ctcATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGgaga                                                                     >  2:2421556/1‑93 (MQ=255)
                          tcATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGAt                                                                    <  1:6245852/93‑1 (MQ=255)
                           cATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATg                                                                   <  2:4048433/93‑1 (MQ=255)
                            aTTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATgg                                                                  <  2:3427137/93‑1 (MQ=255)
                            aTTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATgg                                                                  >  2:45894/1‑93 (MQ=255)
                             ttCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgc                                                                 >  2:5187668/1‑93 (MQ=255)
                              tCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGgcg                                                                >  2:42688/1‑93 (MQ=255)
                              gcccGCCATTGACTGGTCAAAAGTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcg                                                                <  1:186203/92‑1 (MQ=255)
                                ccGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                              >  1:6063185/1‑93 (MQ=255)
                                ccGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                              >  2:2962397/1‑93 (MQ=255)
                                ccGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                              <  1:897583/93‑1 (MQ=255)
                                  gCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCa                                                            >  1:2223731/1‑93 (MQ=255)
                                   ccATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCaa                                                           >  2:3153939/1‑93 (MQ=255)
                                    cATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAc                                                          >  1:5217308/1‑93 (MQ=255)
                                    cATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAc                                                          <  1:3584997/93‑1 (MQ=255)
                                     aTTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAcc                                                         <  1:6447721/93‑1 (MQ=255)
                                     aTTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAcc                                                         >  2:3632241/1‑93 (MQ=255)
                                      tttACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCa                                                        >  1:2360715/1‑93 (MQ=255)
                                      tttACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCa                                                        >  2:2581633/1‑93 (MQ=255)
                                       ttACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAg                                                       <  1:4631330/93‑1 (MQ=255)
                                        tACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGa                                                      >  1:651319/1‑93 (MQ=255)
                                         aCTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGAt                                                     <  1:497153/93‑1 (MQ=255)
                                            gggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCtt                                                  >  2:102017/1‑93 (MQ=255)
                                            gggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCtt                                                  <  2:3205628/93‑1 (MQ=255)
                                             ggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTc                                                 <  1:5870543/93‑1 (MQ=255)
                                             ggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTc                                                 <  2:1658880/93‑1 (MQ=255)
                                               cAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCt                                               <  1:5933950/93‑1 (MQ=255)
                                               cAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCt                                               <  1:134994/93‑1 (MQ=255)
                                                  aaTTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGCCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGAt                                            >  2:5560725/1‑93 (MQ=255)
                                                    ttCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTa                                          >  1:2922992/1‑93 (MQ=255)
                                                     tCAGTCCCGGTTGCTAAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTaa                                         <  2:6101394/93‑1 (MQ=255)
                                                      cAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAAt                                        >  2:5284541/1‑93 (MQ=255)
                                                        gTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGa                                      <  2:4455825/93‑1 (MQ=255)
                                                        gTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGa                                      <  2:6256920/93‑1 (MQ=255)
                                                          cccGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAAc                                    <  2:5375144/93‑1 (MQ=255)
                                                           ccGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACg                                   >  1:1611026/1‑93 (MQ=255)
                                                           ccGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACg                                   >  2:5954324/1‑93 (MQ=255)
                                                            cGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGc                                  >  1:6200863/1‑93 (MQ=255)
                                                            cGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGc                                  <  1:3875261/93‑1 (MQ=255)
                                                             gggTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGcc                                 >  1:3946554/1‑93 (MQ=255)
                                                              ggTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCa                                >  1:538977/1‑93 (MQ=255)
                                                               gTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAg                               <  2:3366874/93‑1 (MQ=255)
                                                                tGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGc                              >  2:728008/1‑93 (MQ=255)
                                                                  cTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctc                            >  2:1335331/1‑93 (MQ=255)
                                                                   tGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctct                           >  2:1949526/1‑93 (MQ=255)
                                                                   tGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctct                           >  1:810787/1‑93 (MQ=255)
                                                                    gAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCtctc                          >  1:3763073/1‑93 (MQ=255)
                                                                    gAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCtctc                          >  1:6356130/1‑93 (MQ=255)
                                                                    gAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCtctc                          <  1:3661365/93‑1 (MQ=255)
                                                                     aaaGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTcc                         <  2:6318735/93‑1 (MQ=255)
                                                                     aaaGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTcc                         <  2:4322485/93‑1 (MQ=255)
                                                                       aGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGt                       >  1:2501334/1‑93 (MQ=255)
                                                                       aGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCCGCTCTCCGt                       >  1:4015648/1‑93 (MQ=255)
                                                                          cGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaa                    <  1:2992860/93‑1 (MQ=255)
                                                                          cGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaa                    <  1:4289726/93‑1 (MQ=255)
                                                                          cGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaa                    <  2:569528/93‑1 (MQ=255)
                                                                            aTCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaatc                  <  2:2307150/93‑1 (MQ=255)
                                                                               aGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaac               >  2:2199570/1‑93 (MQ=255)
                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              >  2:4893309/1‑93 (MQ=255)
                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              <  1:4867024/93‑1 (MQ=255)
                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              >  2:6500772/1‑93 (MQ=255)
                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              <  2:6591122/93‑1 (MQ=255)
                                                                                   aTTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACaaaa           >  1:3500606/1‑93 (MQ=255)
                                                                                   aTTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACaaaa           <  1:1647344/93‑1 (MQ=255)
                                                                                    ttCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAc          <  1:6587756/93‑1 (MQ=255)
                                                                                    ttCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAc          <  2:1603982/93‑1 (MQ=255)
                                                                                     tCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACt         <  1:4565954/93‑1 (MQ=255)
                                                                                     tCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACt         <  1:2612494/93‑1 (MQ=255)
                                                                                      cGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTg        <  2:2102601/93‑1 (MQ=255)
                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      <  1:14702/93‑1 (MQ=255)
                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      <  1:3300903/93‑1 (MQ=255)
                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      >  1:3156469/1‑93 (MQ=255)
                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      <  1:69875/93‑1 (MQ=255)
                                                                                         cGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCt     <  1:4249818/93‑1 (MQ=255)
                                                                                            |                                                                                            
ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG  >  NC_000913/2621250‑2621434

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ 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
TATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGT  >  NC_000913/2621249‑2621435
                                                                                             |                                                                                             
TATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGC                                                                                               >  SRR4302547.5728519/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTC                                                                                              >  SRR4302547.4983042/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTC                                                                                              >  SRR4302547.2364850/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTC                                                                                              >  SRR4302547.2691026/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTC                                                                                              >  SRR4302547.4250255/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTC                                                                                              >  SRR4302547.3341476/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTC                                                                                              >  SRR4302547.3764403/1‑94 (MQ=60)
 ATTAGCGAAGAACAAGGCGCCAACTTTCATTCCGGCCATTTAGTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGCTCAGGAATTCGCCGGTC                                                                                              >  SRR4302547.3341477/1‑94 (MQ=60)
  TTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCT                                                                                             <  SRR4302547.1136745/94‑1 (MQ=60)
  TTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCT                                                                                             >  SRR4302547.2349703/1‑94 (MQ=60)
   TAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTC                                                                                            >  SRR4302547.3405212/1‑94 (MQ=60)
   TAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTC                                                                                            <  SRR4302547.3601183/94‑1 (MQ=60)
    AGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAACATTCAGTCCCGGGTGCTGAAAGCCGATCAGGTATTCGACGGTCTCT                                                                                           >  SRR4302547.1110626/1‑94 (MQ=60)
    AGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCT                                                                                           <  SRR4302547.810870/94‑1 (MQ=60)
     GCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTA                                                                                          <  SRR4302547.4829579/94‑1 (MQ=60)
      CGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTAC                                                                                         <  SRR4302547.581250/94‑1 (MQ=60)
      CGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTAC                                                                                         >  SRR4302547.6607947/1‑94 (MQ=60)
      CGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTAC                                                                                         <  SRR4302547.116811/94‑1 (MQ=60)
        AAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAA                                                                                       >  SRR4302547.1671478/1‑94 (MQ=60)
         AGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAAC                                                                                      >  SRR4302547.3649194/1‑94 (MQ=60)
          GAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACG                                                                                     >  SRR4302547.3251983/1‑94 (MQ=60)
          GAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACG                                                                                     >  SRR4302547.6327298/1‑94 (MQ=60)
          GAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACG                                                                                     >  SRR4302547.3075657/1‑94 (MQ=60)
           AACAAGGCTCCACCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGA                                                                                    <  SRR4302547.3195741/94‑1 (MQ=60)
           AACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGA                                                                                    >  SRR4302547.1635797/1‑94 (MQ=60)
           AACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGA                                                                                    >  SRR4302547.4581161/1‑94 (MQ=60)
            ACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAG                                                                                   <  SRR4302547.2007542/94‑1 (MQ=60)
            ACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAG                                                                                   <  SRR4302547.5561090/94‑1 (MQ=60)
             CAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGC                                                                                  >  SRR4302547.5537089/1‑94 (MQ=60)
             CAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGC                                                                                  >  SRR4302547.4555806/1‑94 (MQ=60)
              AAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCT                                                                                 >  SRR4302547.3131876/1‑94 (MQ=60)
              AAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCT                                                                                 <  SRR4302547.5776884/94‑1 (MQ=60)
               AGGCTCCACCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTA                                                                                <  SRR4302547.1800904/94‑1 (MQ=60)
               AGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTA                                                                                >  SRR4302547.3692472/1‑94 (MQ=60)
                GGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTAT                                                                               <  SRR4302547.4876915/94‑1 (MQ=60)
                 GCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATT                                                                              <  SRR4302547.589643/94‑1 (MQ=60)
                  CTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAAatgtcaagat                                                                             <  SRR4302547.6581516/94‑11 (MQ=60)
                  CTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTG                                                                             >  SRR4302547.2899659/1‑94 (MQ=60)
                  CTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTG                                                                             <  SRR4302547.248048/94‑1 (MQ=60)
                  CTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTG                                                                             >  SRR4302547.3890642/1‑94 (MQ=60)
                   TCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGC                                                                            >  SRR4302547.1141662/1‑94 (MQ=60)
                      AACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGG                                                                         >  SRR4302547.5778412/1‑94 (MQ=60)
                      AACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGG                                                                         >  SRR4302547.2590375/1‑94 (MQ=60)
                       CCTCTAATTCCCGCCAATAACTGGGCAAAATTCAGCCACGGGTGCTGAAAGCCGAAAAGGAATTCGACGGACTCTACAACGAGCTATTGCTAGA                                                                        <  SRR4302547.5864160/94‑1 (MQ=55)
                       ACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGA                                                                        <  SRR4302547.347365/94‑1 (MQ=60)
                        CTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAG                                                                       <  SRR4302547.566726/94‑1 (MQ=60)
                        CTCTCATTCCCGCCATGTACTGGTCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAG                                                                       <  SRR4302547.5906582/94‑1 (MQ=60)
                        CTCTCATGCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAG                                                                       <  SRR4302547.5953336/94‑1 (MQ=60)
                         TCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGA                                                                      <  SRR4302547.823120/94‑1 (MQ=60)
                         TCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGA                                                                      >  SRR4302547.1533646/1‑94 (MQ=60)
                           TCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATG                                                                    <  SRR4302547.6375597/94‑1 (MQ=60)
                            CATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGCGATGG                                                                   >  SRR4302547.3984701/1‑94 (MQ=60)
                             ATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGC                                                                  <  SRR4302547.2420371/94‑1 (MQ=60)
                             ATCCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGC                                                                  <  SRR4302547.5170799/94‑1 (MQ=60)
                              TTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCG                                                                 <  SRR4302547.5718346/94‑1 (MQ=60)
                               TCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGC                                                                <  SRR4302547.964346/94‑1 (MQ=60)
                               TCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGC                                                                <  SRR4302547.5025125/94‑1 (MQ=60)
                               GCCCGCCATTGACTGGTCAAAAGTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGC                                                                <  SRR4302547.189708/94‑1 (MQ=60)
                                CCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCG                                                               >  SRR4302547.6188952/1‑94 (MQ=60)
                                 CCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGC                                                              <  SRR4302547.914839/94‑1 (MQ=60)
                                  CGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCA                                                             >  SRR4302547.2268534/1‑94 (MQ=60)
                                  CGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCA                                                             <  SRR4302547.6356985/94‑1 (MQ=60)
                                   GCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAA                                                            >  SRR4302547.4875392/1‑94 (MQ=60)
                                    CCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAC                                                           >  SRR4302547.6087713/1‑94 (MQ=60)
                                    CCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAC                                                           >  SRR4302547.5324700/1‑94 (MQ=60)
                                     CATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACC                                                          <  SRR4302547.968095/94‑1 (MQ=60)
                                     CATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACC                                                          <  SRR4302547.4066146/94‑1 (MQ=60)
                                     CATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACC                                                          <  SRR4302547.3658143/94‑1 (MQ=60)
                                     CATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACC                                                          <  SRR4302547.1339079/94‑1 (MQ=60)
                                      ATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCA                                                         >  SRR4302547.2408574/1‑94 (MQ=60)
                                      ATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCA                                                         <  SRR4302547.6582054/94‑1 (MQ=60)
                                       TTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAG                                                        <  SRR4302547.4357013/94‑1 (MQ=60)
                                        TTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGA                                                       <  SRR4302547.4727260/94‑1 (MQ=60)
                                        TTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGA                                                       >  SRR4302547.663312/1‑94 (MQ=60)
                                         TACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGAT                                                      <  SRR4302547.3036636/94‑1 (MQ=60)
                                          ACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATC                                                     <  SRR4302547.506216/94‑1 (MQ=60)
                                           CTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCT                                                    <  SRR4302547.6509895/94‑1 (MQ=60)
                                            TGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTT                                                   <  SRR4302547.5250615/94‑1 (MQ=60)
                                              GGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCC                                                 <  SRR4302547.5992055/94‑1 (MQ=60)
                                                CAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTG                                               <  SRR4302547.6056992/94‑1 (MQ=60)
                                                CAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTG                                               <  SRR4302547.137432/94‑1 (MQ=60)
                                                 AAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGA                                              >  SRR4302547.1251267/1‑94 (MQ=60)
                                                   AATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATT                                            >  SRR4302547.6391131/1‑94 (MQ=60)
                                                    ATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTA                                           >  SRR4302547.2981270/1‑94 (MQ=60)
                                                       CAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATG                                        <  SRR4302547.4837264/94‑1 (MQ=60)
                                                       CAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATG                                        <  SRR4302547.4712563/94‑1 (MQ=60)
                                                         GTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAA                                      >  SRR4302547.6193747/1‑94 (MQ=60)
                                                           CCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACG                                    >  SRR4302547.1642941/1‑94 (MQ=60)
                                                            CCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGC                                   >  SRR4302547.6329609/1‑94 (MQ=60)
                                                             CGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCC                                  >  SRR4302547.4026816/1‑94 (MQ=60)
                                                             CGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCC                                  <  SRR4302547.3954114/94‑1 (MQ=60)
                                                              GGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCA                                 >  SRR4302547.548708/1‑94 (MQ=60)
                                                                GTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGC                               <  SRR4302547.3603718/94‑1 (MQ=60)
                                                                  GCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCGTGATTAATGAACGCCAGCTC                             <  SRR4302547.6632130/94‑1 (MQ=60)
                                                                   CTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCT                            >  SRR4302547.4934768/1‑94 (MQ=60)
                                                                   CTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCT                            >  SRR4302547.826151/1‑94 (MQ=60)
                                                                    TGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTC                           >  SRR4302547.6488646/1‑94 (MQ=60)
                                                                    TGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTC                           >  SRR4302547.3839832/1‑94 (MQ=60)
                                                                     GAAAGCCGATCAGGAATTCGACGTTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCC                          >  SRR4302547.6027228/1‑94 (MQ=60)
                                                                     GAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCC                          <  SRR4302547.3736049/94‑1 (MQ=60)
                                                                     GAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCC                          <  SRR4302547.1169334/94‑1 (MQ=60)
                                                                     GAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCC                          >  SRR4302547.3577547/1‑94 (MQ=60)
                                                                      AAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCG                         <  SRR4302547.1968432/94‑1 (MQ=60)
                                                                      AAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCG                         >  SRR4302547.2813413/1‑94 (MQ=60)
                                                                      AAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCG                         >  SRR4302547.4220026/1‑94 (MQ=60)
                                                                       AAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGT                        >  SRR4302547.2551870/1‑94 (MQ=60)
                                                                       AAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCCGCTCTCCGT                        >  SRR4302547.4097251/1‑94 (MQ=60)
                                                                          CCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAA                     >  SRR4302547.4025305/1‑94 (MQ=60)
                                                                          CCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAA                     <  SRR4302547.6048538/94‑1 (MQ=60)
                                                                           CGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAAT                    <  SRR4302547.3052395/94‑1 (MQ=60)
                                                                           CGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAAT                    <  SRR4302547.4377894/94‑1 (MQ=60)
                                                                            GATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATC                   >  SRR4302547.5799373/1‑94 (MQ=60)
                                                                             ATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCA                  >  SRR4302547.4679679/1‑94 (MQ=60)
                                                                              TCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAA                 <  SRR4302547.6433709/94‑1 (MQ=60)
                                                                              TCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAA                 <  SRR4302547.1397027/94‑1 (MQ=60)
                                                                               CAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAAC                <  SRR4302547.707890/94‑1 (MQ=60)
                                                                               CAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAAC                >  SRR4302547.2027398/1‑94 (MQ=60)
                                                                               CAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAAC                >  SRR4302547.4238300/1‑94 (MQ=60)
                                                                                AGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACA               <  SRR4302547.5894318/94‑1 (MQ=60)
                                                                                AGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACA               >  SRR4302547.4380941/1‑94 (MQ=60)
                                                                                 GGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAA              <  SRR4302547.3804388/94‑1 (MQ=60)
                                                                                 GGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAA              <  SRR4302547.4967884/94‑1 (MQ=60)
                                                                                   AATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAA            <  SRR4302547.3838829/94‑1 (MQ=60)
                                                                                   AATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAA            >  SRR4302547.3053064/1‑94 (MQ=60)
                                                                                   AATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAA            >  SRR4302547.3572042/1‑94 (MQ=60)
                                                                                    ATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAC           >  SRR4302547.573921/1‑94 (MQ=60)
                                                                                    ATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAC           <  SRR4302547.1679875/94‑1 (MQ=60)
                                                                                    ATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAC           >  SRR4302547.6653489/1‑94 (MQ=60)
                                                                                     TTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACT          >  SRR4302547.5472486/1‑94 (MQ=60)
                                                                                     TTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACT          <  SRR4302547.6724903/94‑1 (MQ=60)
                                                                                      TCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTG         <  SRR4302547.4660425/94‑1 (MQ=60)
                                                                                      TCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTG         <  SRR4302547.2665137/94‑1 (MQ=60)
                                                                                        GACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGC       >  SRR4302547.4605639/1‑94 (MQ=60)
                                                                                        GACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGC       >  SRR4302547.3219241/1‑94 (MQ=60)
                                                                                         ACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCT      <  SRR4302547.4767656/94‑1 (MQ=60)
                                                                                         ACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCT      <  SRR4302547.71179/94‑1 (MQ=60)
                                                                                         ACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCT      <  SRR4302547.3367343/94‑1 (MQ=60)
                                                                                         ACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCT      <  SRR4302547.14975/94‑1 (MQ=60)
                                                                                         ACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCT      <  SRR4302547.2097816/94‑1 (MQ=60)
                                                                                          CGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTG     <  SRR4302547.4336763/94‑1 (MQ=60)
                                                                                          CGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTG     >  SRR4302547.720623/1‑94 (MQ=60)
                                                                                          CGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTG     <  SRR4302547.90743/94‑1 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGG    >  SRR4302547.5889159/1‑94 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGC    >  SRR4302547.4268680/1‑94 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGC    >  SRR4302547.6513267/1‑94 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGC    >  SRR4302547.2138888/1‑94 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGC    >  SRR4302547.1598714/1‑94 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGC    >  SRR4302547.3708944/1‑94 (MQ=60)
                                                                                           GGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCATGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGC    >  SRR4302547.4086116/1‑94 (MQ=60)
                                                                                            GTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG   <  SRR4302547.5331881/94‑1 (MQ=60)
                                                                                            GTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG   >  SRR4302547.121498/1‑94 (MQ=60)
                                                                                            GCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG   >  SRR4302547.1831550/1‑94 (MQ=60)
                                                                                             CCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGT  >  SRR4302547.3185747/1‑94 (MQ=60)
                                                                                             |                                                                                             
TATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGT  >  NC_000913/2621249‑2621435

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