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 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,429 C→T A186A (GCC→GCT wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,4290CT99.6% 946.4 / ‑5.3 253A186A (GCC→GCTwecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (130/122);  total (131/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 260 total aligned reads displayed.

TGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAA  >  NC_000913/3969340‑3969517
                                                                                         |                                                                                        
tGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTcgc                                                                                       >  1:5563290/1‑93 (MQ=255)
tGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTcgc                                                                                       <  1:4287254/93‑1 (MQ=255)
   ccGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGcc                                                                                    <  1:3731133/93‑1 (MQ=255)
   ccGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGcc                                                                                    <  1:3086165/93‑1 (MQ=255)
    cGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCa                                                                                   >  2:1335280/1‑93 (MQ=255)
     ggAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAg                                                                                  <  2:4790026/93‑1 (MQ=255)
     ggAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAg                                                                                  <  1:4848880/93‑1 (MQ=255)
      gAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGc                                                                                 >  2:3273034/1‑93 (MQ=255)
       aaCGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCg                                                                                >  2:2092377/1‑93 (MQ=255)
        aCGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGa                                                                               <  2:3647221/93‑1 (MQ=255)
          gcgTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAAc                                                                             >  1:5820504/1‑93 (MQ=255)
           cgTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACt                                                                            >  2:4267696/1‑93 (MQ=255)
            gtgtTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTg                                                                           <  2:4906242/93‑1 (MQ=255)
             tgtTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGt                                                                          <  1:1786146/93‑1 (MQ=255)
              gtTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTa                                                                         >  2:5349151/1‑93 (MQ=255)
               ttACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTAc                                                                        <  2:2463676/93‑1 (MQ=255)
               ttACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTAc                                                                        >  2:1843468/1‑93 (MQ=255)
                tACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACa                                                                       <  1:3520158/93‑1 (MQ=255)
                  ccAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACaaa                                                                     <  2:1245868/93‑1 (MQ=255)
                   cAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACaaaa                                                                    <  1:5401586/93‑1 (MQ=255)
                   cAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACaaaa                                                                    >  2:3795000/1‑93 (MQ=255)
                      gACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAAttt                                                                 >  1:1896407/1‑93 (MQ=255)
                      gACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAAttt                                                                 <  2:5246921/93‑1 (MQ=255)
                        cAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTc                                                               >  1:4649351/1‑93 (MQ=255)
                         aGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTcc                                                              <  2:1891275/93‑1 (MQ=255)
                         aGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTcc                                                              <  1:1152564/93‑1 (MQ=255)
                          ggTAATGGTCGAGCTGATTAAAAACGATCGCGTTGTTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCt                                                             >  2:5193283/1‑93 (MQ=255)
                          ggTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCt                                                             <  1:4359723/93‑1 (MQ=255)
                           gTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTc                                                            <  1:1568962/93‑1 (MQ=255)
                           gTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTc                                                            >  2:1433433/1‑93 (MQ=255)
                            tAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCg                                                           >  2:5466760/1‑93 (MQ=255)
                              aTGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGaa                                                         <  1:2397277/93‑1 (MQ=255)
                              aTGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGaa                                                         <  1:473388/93‑1 (MQ=255)
                              aTGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGaa                                                         <  1:5951362/93‑1 (MQ=255)
                               tGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAg                                                        <  2:3587706/93‑1 (MQ=255)
                               tGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAg                                                        >  1:1223816/1‑93 (MQ=255)
                               tGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAg                                                        <  2:570316/93‑1 (MQ=255)
                                tgTCCAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAgg                                                       <  1:1679117/92‑1 (MQ=255)
                                 gTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGt                                                      >  1:3677736/1‑93 (MQ=255)
                                   cGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGa                                                    >  1:1586994/1‑93 (MQ=255)
                                   cGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTGGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGa                                                    <  2:446274/93‑1 (MQ=255)
                                     aGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgt                                                  >  2:4298169/1‑93 (MQ=255)
                                     aGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGCGAgg                                                  >  2:5067643/1‑92 (MQ=255)
                                      gCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtg                                                 >  2:4137144/1‑93 (MQ=255)
                                      gCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtg                                                 >  1:2804861/1‑93 (MQ=255)
                                      gCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtg                                                 >  1:491058/1‑93 (MQ=255)
                                      gCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtg                                                 >  1:2060430/1‑93 (MQ=255)
                                       cTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgt                                                <  1:242265/93‑1 (MQ=255)
                                        tGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgtg                                               >  1:1409588/1‑93 (MQ=255)
                                        tGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgtg                                               >  2:2324398/1‑93 (MQ=255)
                                         gATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgtgt                                              <  1:548765/93‑1 (MQ=255)
                                         gATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgtgt                                              >  1:4600292/1‑93 (MQ=255)
                                         gATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgtgt                                              >  2:1089878/1‑93 (MQ=255)
                                          aTTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTgtc                                             <  2:4603267/93‑1 (MQ=255)
                                             aaaaaCGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTgtcgtc                                          <  1:855373/93‑1 (MQ=255)
                                             aaaaaCGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTgtcgtc                                          <  1:3457919/93‑1 (MQ=255)
                                              aaaaCGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCa                                         >  2:2810093/1‑93 (MQ=255)
                                              aaaaCGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCa                                         >  1:3572767/1‑93 (MQ=255)
                                               aaaCGATCGCGTGATTGGTGGTATGACGCCGGTTTGTGCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCAc                                        <  1:3236193/93‑1 (MQ=255)
                                               aaaCGATCGCGGGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCAc                                        >  1:2453611/1‑93 (MQ=255)
                                                  cGATCGCGTGATTGGTGGTATGACGACGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTaa                                     <  2:5993972/93‑1 (MQ=255)
                                                    aTCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACt                                   >  1:5420299/1‑93 (MQ=255)
                                                     tCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTCCAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTc                                  <  1:2486578/93‑1 (MQ=255)
                                                       gcgTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTcgc                                >  1:5221135/1‑93 (MQ=255)
                                                       gcgTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTcgc                                >  1:4384652/1‑93 (MQ=255)
                                                       gcgTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTcgc                                >  1:2715162/1‑93 (MQ=255)
                                                        cgTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTcgcg                               >  1:562545/1‑93 (MQ=255)
                                                         gTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCgg                              >  2:1708650/1‑93 (MQ=255)
                                                          tGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCggc                             >  1:4421934/1‑92 (MQ=255)
                                                            aTTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAcg                           <  2:2902394/93‑1 (MQ=255)
                                                             ttGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAcgg                          >  1:3110698/1‑93 (MQ=255)
                                                             ttGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAcgg                          >  1:3997821/1‑93 (MQ=255)
                                                             ttGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAcgg                          >  2:3842406/1‑93 (MQ=255)
                                                              tggtggTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAcggc                         >  1:32239/1‑93 (MQ=255)
                                                              tggtggTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAcggc                         <  2:2621242/93‑1 (MQ=255)
                                                                 tggtATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGa                      >  1:4221239/1‑93 (MQ=255)
                                                                   gtATGACGCCGTTTTGTTCGGCTCGCGCAAGCGAACTGTACAAAATTTTCCTCGAAGGTGGGTGTGTGGTCACTAACTCGCGGACGGCGGaaa                    <  2:693606/93‑1 (MQ=255)
                                                                   gtATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGaaa                    <  1:1108388/93‑1 (MQ=255)
                                                                     aTGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAtg                  <  2:620612/93‑1 (MQ=255)
                                                                     aTGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAtg                  <  2:3494067/93‑1 (MQ=255)
                                                                     aTGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAtg                  <  1:258345/93‑1 (MQ=255)
                                                                      tGACGCCTGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAtgt                 <  2:749179/93‑1 (MQ=255)
                                                                      tGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAtgt                 >  2:1859729/1‑93 (MQ=255)
                                                                        aCGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATgtgt               <  1:1323748/93‑1 (MQ=255)
                                                                        aCGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATgtgt               <  1:5322916/93‑1 (MQ=255)
                                                                         cGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTa              >  1:964852/1‑93 (MQ=255)
                                                                         cGCCGGGTTGGTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGGGTCGTCACTAACTCGCGGACGGCGGAAATGTGTa              <  2:1614117/93‑1 (MQ=255)
                                                                          gCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTaa             <  2:2059293/93‑1 (MQ=255)
                                                                           ccGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAg            >  2:3193021/1‑93 (MQ=255)
                                                                            cGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGc           <  2:5775550/93‑1 (MQ=255)
                                                                             ggTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCt          <  2:5918224/93‑1 (MQ=255)
                                                                             ggTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCt          >  1:2921146/1‑93 (MQ=255)
                                                                              gTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTc         >  1:448837/1‑93 (MQ=255)
                                                                              gTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTc         >  2:5580662/1‑93 (MQ=255)
                                                                                ttGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCAc       <  2:101994/93‑1 (MQ=255)
                                                                                ttGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCAc       <  1:4960132/93‑1 (MQ=255)
                                                                                ttGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCAc       <  1:2010480/93‑1 (MQ=255)
                                                                                 tGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCAcc      >  1:6131343/1‑93 (MQ=255)
                                                                                 tGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCAcc      >  2:5015825/1‑93 (MQ=255)
                                                                                  gTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCg     >  1:5270891/1‑93 (MQ=255)
                                                                                         |                                                                                        
TGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAA  >  NC_000913/3969340‑3969517

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACA  >  NC_000913/3969336‑3969520
                                                                                             |                                                                                           
TTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCT                                                                                             <  SRR4302533.5118702/94‑1 (MQ=60)
 TACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTC                                                                                            >  SRR4302533.4866615/1‑94 (MQ=60)
 TACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTC                                                                                            >  SRR4302533.3568048/1‑94 (MQ=60)
   CTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGC                                                                                          >  SRR4302533.4380902/1‑94 (MQ=60)
   CTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGC                                                                                          >  SRR4302533.5596916/1‑94 (MQ=60)
    TGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCG                                                                                         <  SRR4302533.4313035/94‑1 (MQ=60)
       CCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCA                                                                                      <  SRR4302533.3754019/94‑1 (MQ=60)
       CCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCA                                                                                      <  SRR4302533.3105137/94‑1 (MQ=60)
         GGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGC                                                                                    <  SRR4302533.4878265/94‑1 (MQ=60)
         GGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGC                                                                                    >  SRR4302533.2121174/1‑94 (MQ=60)
         GGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGC                                                                                    <  SRR4302533.3235275/94‑1 (MQ=60)
          GAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCG                                                                                   >  SRR4302533.542086/1‑94 (MQ=60)
          GAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCG                                                                                   <  SRR4302533.1609048/94‑1 (MQ=60)
           AACGCGTGTTACCAGGACAGGTAATGGTCGCGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGA                                                                                  >  SRR4302533.1297221/1‑94 (MQ=60)
           AACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGA                                                                                  >  SRR4302533.6054726/1‑94 (MQ=60)
             CGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAAC                                                                                >  SRR4302533.5855664/1‑94 (MQ=60)
             CGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAAC                                                                                >  SRR4302533.1733428/1‑94 (MQ=60)
             CGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCACGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAAC                                                                                <  SRR4302533.978163/94‑1 (MQ=60)
               CGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTG                                                                              >  SRR4302533.2769096/1‑94 (MQ=60)
                GTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGT                                                                             >  SRR4302533.259495/1‑94 (MQ=60)
                 TGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTA                                                                            <  SRR4302533.1797006/94‑1 (MQ=60)
                  GTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTAC                                                                           >  SRR4302533.1486630/1‑94 (MQ=60)
                  GTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTAC                                                                           >  SRR4302533.3874426/1‑94 (MQ=60)
                    TACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAA                                                                         <  SRR4302533.3541859/94‑1 (MQ=60)
                    TACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAA                                                                         >  SRR4302533.4457957/1‑94 (MQ=60)
                       CAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAAT                                                                      >  SRR4302533.3788462/1‑94 (MQ=60)
                       CAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAAT                                                                      <  SRR4302533.2865599/94‑1 (MQ=60)
                       CAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAAT                                                                      <  SRR4302533.5434045/94‑1 (MQ=60)
                        AGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATT                                                                     >  SRR4302533.4360163/1‑94 (MQ=60)
                         GGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTT                                                                    >  SRR4302533.2828943/1‑94 (MQ=60)
                         GGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTT                                                                    >  SRR4302533.1907983/1‑94 (MQ=60)
                         GGACAGGTAATGGTCGAGCTGATTAAAAACGATCCCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTT                                                                    <  SRR4302533.1947871/94‑1 (MQ=60)
                          GACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTT                                                                   <  SRR4302533.1557532/94‑1 (MQ=60)
                           ACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTC                                                                  >  SRR4302533.4677520/1‑94 (MQ=60)
                            CAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCC                                                                 >  SRR4302533.6168151/1‑94 (MQ=60)
                            CAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCC                                                                 >  SRR4302533.2690483/1‑94 (MQ=60)
                            CAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCC                                                                 <  SRR4302533.2412499/94‑1 (MQ=60)
                            CAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCC                                                                 >  SRR4302533.1392077/1‑94 (MQ=60)
                             AGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCT                                                                >  SRR4302533.4705527/1‑94 (MQ=60)
                             AGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCT                                                                >  SRR4302533.1593085/1‑94 (MQ=60)
                             AGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCT                                                                <  SRR4302533.1159467/94‑1 (MQ=60)
                              GGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTC                                                               <  SRR4302533.4386039/94‑1 (MQ=60)
                               GTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCG                                                              <  SRR4302533.1578410/94‑1 (MQ=60)
                                  ATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAG                                                           >  SRR4302533.1231155/1‑94 (MQ=60)
                                  ATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAG                                                           <  SRR4302533.2412076/94‑1 (MQ=60)
                                  ATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAG                                                           <  SRR4302533.476217/94‑1 (MQ=60)
                                  ATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAG                                                           <  SRR4302533.4434179/94‑1 (MQ=60)
                                  ATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAG                                                           <  SRR4302533.5987359/94‑1 (MQ=60)
                                    tgtccAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGT                                                         <  SRR4302533.1689287/89‑1 (MQ=60)
                                    GGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGT                                                         <  SRR4302533.3708731/94‑1 (MQ=60)
                                    GGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGT                                                         >  SRR4302533.3700343/1‑94 (MQ=60)
                                    GGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGT                                                         <  SRR4302533.5344227/94‑1 (MQ=60)
                                    GGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGT                                                         <  SRR4302533.3101196/94‑1 (MQ=60)
                                    GGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGCACAAAATTTTCCTCGAAGGG                                                         >  SRR4302533.484002/1‑94 (MQ=60)
                                      TCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGA                                                       >  SRR4302533.1596566/1‑94 (MQ=60)
                                       CGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAG                                                      >  SRR4302533.1188598/1‑94 (MQ=60)
                                       CGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAG                                                      >  SRR4302533.4642482/1‑94 (MQ=60)
                                         AGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTG                                                    >  SRR4302533.2822125/1‑94 (MQ=60)
                                         AGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTG                                                    >  SRR4302533.493992/1‑94 (MQ=60)
                                         AGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTG                                                    >  SRR4302533.2073136/1‑94 (MQ=60)
                                         AGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTG                                                    <  SRR4302533.1891150/94‑1 (MQ=60)
                                          GCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGT                                                   <  SRR4302533.4485728/94‑1 (MQ=60)
                                           CTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTG                                                  <  SRR4302533.4572421/94‑1 (MQ=60)
                                           CTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTG                                                  >  SRR4302533.1418110/1‑94 (MQ=60)
                                           CTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTG                                                  >  SRR4302533.5730044/1‑94 (MQ=60)
                                           CTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTG                                                  <  SRR4302533.243766/94‑1 (MQ=60)
                                            TGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGT                                                 >  SRR4302533.4628161/1‑94 (MQ=60)
                                            TGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGT                                                 <  SRR4302533.3319028/94‑1 (MQ=60)
                                            TGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGT                                                 >  SRR4302533.1888235/1‑94 (MQ=60)
                                             GATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTC                                                >  SRR4302533.1384295/1‑94 (MQ=60)
                                             GATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTC                                                <  SRR4302533.552020/94‑1 (MQ=60)
                                              ATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCG                                               <  SRR4302533.5272681/94‑1 (MQ=60)
                                              ATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCG                                               >  SRR4302533.5437970/1‑94 (MQ=60)
                                                 AAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCA                                            >  SRR4302533.4818229/1‑94 (MQ=60)
                                                 AAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCA                                            <  SRR4302533.860386/94‑1 (MQ=60)
                                                 AAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCA                                            >  SRR4302533.3254721/1‑94 (MQ=60)
                                                 AAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCA                                            >  SRR4302533.3594780/1‑94 (MQ=60)
                                                 AAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCA                                            <  SRR4302533.3479265/94‑1 (MQ=60)
                                                  AAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCAC                                           >  SRR4302533.3129574/1‑94 (MQ=60)
                                                  AAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCAC                                           >  SRR4302533.5523379/1‑94 (MQ=60)
                                                  AAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGGGAGTGTGTCGTCAC                                           >  SRR4302533.5432750/1‑94 (MQ=60)
                                                  AAAACGATCGCGGGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCAC                                           >  SRR4302533.2468751/1‑94 (MQ=60)
                                                   AAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACT                                          <  SRR4302533.5059590/94‑1 (MQ=60)
                                                   AAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTGCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACT                                          <  SRR4302533.3256135/94‑1 (MQ=60)
                                                    AACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTA                                         >  SRR4302533.3696786/1‑94 (MQ=60)
                                                      CGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAAT                                       >  SRR4302533.4442758/1‑94 (MQ=60)
                                                       GATCGCGTGATTGGTGGTATTACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACT                                      >  SRR4302533.5701267/1‑94 (MQ=60)
                                                       GATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAATT                                      >  SRR4302533.1657742/1‑94 (MQ=60)
                                                       GATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACT                                      >  SRR4302533.4942825/1‑94 (MQ=60)
                                                       GATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACT                                      >  SRR4302533.5452915/1‑94 (MQ=60)
                                                       GATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACT                                      <  SRR4302533.3483170/94‑1 (MQ=60)
                                                       GATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACT                                      >  SRR4302533.3436330/1‑94 (MQ=60)
                                                        ATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTC                                     <  SRR4302533.2940030/94‑1 (MQ=60)
                                                        ATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTC                                     <  SRR4302533.5653312/94‑1 (MQ=60)
                                                         TCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTCCAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCG                                    <  SRR4302533.2501897/94‑1 (MQ=60)
                                                         TCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCG                                    <  SRR4302533.2095047/94‑1 (MQ=60)
                                                          CGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGC                                   >  SRR4302533.2731833/1‑94 (MQ=60)
                                                          CGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGC                                   >  SRR4302533.4411125/1‑94 (MQ=60)
                                                          CGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGC                                   >  SRR4302533.5252554/1‑94 (MQ=60)
                                                           GCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCG                                  >  SRR4302533.565881/1‑94 (MQ=60)
                                                           GCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCG                                  <  SRR4302533.4179073/94‑1 (MQ=60)
                                                           GCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCG                                  <  SRR4302533.5969221/94‑1 (MQ=60)
                                                            CGTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGG                                 >  SRR4302533.4400285/1‑94 (MQ=60)
                                                             GTGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGC                                >  SRR4302533.4448670/1‑94 (MQ=60)
                                                              TGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAC                               <  SRR4302533.49365/94‑1 (MQ=60)
                                                              TGATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGAC                               <  SRR4302533.4665615/94‑1 (MQ=60)
                                                                ATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGG                             >  SRR4302533.982465/1‑94 (MQ=60)
                                                                ATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGG                             >  SRR4302533.3129820/1‑94 (MQ=60)
                                                                ATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGG                             >  SRR4302533.4022161/1‑94 (MQ=60)
                                                                ATTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGG                             >  SRR4302533.3188187/1‑94 (MQ=60)
                                                                 TTGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGC                            >  SRR4302533.32461/1‑94 (MQ=60)
                                                                  TGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCG                           >  SRR4302533.2498573/1‑94 (MQ=60)
                                                                  TGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCG                           <  SRR4302533.2646403/94‑1 (MQ=60)
                                                                  TGGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCG                           <  SRR4302533.3574017/94‑1 (MQ=60)
                                                                   GGTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGG                          >  SRR4302533.806204/1‑94 (MQ=60)
                                                                    GTGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGA                         >  SRR4302533.4246766/1‑94 (MQ=60)
                                                                     TGGTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAA                        >  SRR4302533.5020257/1‑94 (MQ=60)
                                                                       GTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAT                      <  SRR4302533.1114962/94‑1 (MQ=60)
                                                                       GTATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAT                      >  SRR4302533.248487/1‑94 (MQ=60)
                                                                       GTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAT                      >  SRR4302533.2468617/1‑94 (MQ=60)
                                                                        TATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATG                     >  SRR4302533.563260/1‑94 (MQ=60)
                                                                         ATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGT                    >  SRR4302533.1959779/1‑94 (MQ=60)
                                                                         ATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGT                    <  SRR4302533.259917/94‑1 (MQ=60)
                                                                         ATGACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGG                    >  SRR4302533.122571/1‑94 (MQ=60)
                                                                            ACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTA                 <  SRR4302533.1331686/94‑1 (MQ=60)
                                                                            ACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTA                 <  SRR4302533.5354814/94‑1 (MQ=60)
                                                                            ACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTA                 >  SRR4302533.4033047/1‑94 (MQ=60)
                                                                            ACGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTA                 >  SRR4302533.970515/1‑94 (MQ=60)
                                                                             CGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAA                <  SRR4302533.4892053/94‑1 (MQ=60)
                                                                             CGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAA                <  SRR4302533.928541/94‑1 (MQ=60)
                                                                             CGCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCGCGAAGGGGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAA                <  SRR4302533.517059/94‑1 (MQ=60)
                                                                              GCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAG               <  SRR4302533.3529809/94‑1 (MQ=60)
                                                                              GCCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAG               <  SRR4302533.898572/94‑1 (MQ=60)
                                                                               CCGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGC              <  SRR4302533.5453191/94‑1 (MQ=60)
                                                                                CGGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCT             >  SRR4302533.2939039/1‑94 (MQ=60)
                                                                                 GGTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTC            >  SRR4302533.451521/1‑94 (MQ=60)
                                                                                  GTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCA           <  SRR4302533.546439/94‑1 (MQ=60)
                                                                                  GTTTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCA           <  SRR4302533.1125935/94‑1 (MQ=60)
                                                                                    TTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACC         <  SRR4302533.2461068/94‑1 (MQ=60)
                                                                                    TTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACC         <  SRR4302533.2022865/94‑1 (MQ=60)
                                                                                    TTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACC         >  SRR4302533.6168287/1‑94 (MQ=60)
                                                                                    TTGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACC         <  SRR4302533.4990180/94‑1 (MQ=60)
                                                                                     TGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCG        <  SRR4302533.5311372/94‑1 (MQ=60)
                                                                                     TGTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCG        >  SRR4302533.5302504/1‑94 (MQ=60)
                                                                                      GTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGA       >  SRR4302533.2046302/1‑94 (MQ=60)
                                                                                      GTTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGA       >  SRR4302533.3677052/1‑94 (MQ=60)
                                                                                       TTCGGCTCGCGCCAGCGAACTGTACAAAATTTTCCACGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAA      <  SRR4302533.4310416/94‑1 (MQ=60)
                                                                                          GGCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAAC   >  SRR4302533.3001314/1‑94 (MQ=60)
                                                                                           GCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACA  >  SRR4302533.3534799/1‑94 (MQ=60)
                                                                                           GCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACA  <  SRR4302533.1212090/94‑1 (MQ=60)
                                                                                           GCTCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACA  <  SRR4302533.962451/94‑1 (MQ=60)
                                                                                             |                                                                                           
TTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACA  >  NC_000913/3969336‑3969520

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