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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,385,848 G→A D367D (GAC→GAT yjeP ← putative periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,385,8480GA95.9% 727.9 / 6.5 220D367D (GAC→GATyjePputative periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (2/7);  new base A (99/111);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 223 total aligned reads displayed.

GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA  >  NC_000913/4385758‑4385938
                                                                                          |                                                                                          
gCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATc                                                                                          >  2:830132/1‑93 (MQ=255)
gCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATc                                                                                          <  2:572559/93‑1 (MQ=255)
 ccACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCg                                                                                         <  2:554483/93‑1 (MQ=255)
  cACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCgg                                                                                        >  1:4602498/1‑93 (MQ=255)
  cACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCgg                                                                                        <  1:3262023/93‑1 (MQ=255)
   aCCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCCTCGGc                                                                                       >  2:826923/1‑93 (MQ=255)
   aCCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGc                                                                                       >  1:1010888/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  1:2813400/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  2:3602972/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      <  1:1728880/93‑1 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  1:1409786/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAAATCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  2:442846/1‑93 (MQ=255)
       tGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtga                                                                                   >  2:456720/1‑93 (MQ=255)
       tGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtga                                                                                   >  1:4048760/1‑93 (MQ=255)
        gCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtgat                                                                                  <  1:3119779/93‑1 (MQ=255)
        gCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGCCAGCGGCTGACCATCGGCCtgat                                                                                  >  1:4172361/1‑93 (MQ=255)
          aacaaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtgatga                                                                                >  1:1887419/1‑93 (MQ=255)
           acaaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGaa                                                                               >  1:4642854/1‑93 (MQ=255)
           acaaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGaa                                                                               <  1:2703781/93‑1 (MQ=255)
            caaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAt                                                                              <  1:1352048/93‑1 (MQ=255)
            caaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAt                                                                              >  2:924536/1‑93 (MQ=255)
            caaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGAACATCGGCCTGATGAAt                                                                              >  1:1378886/1‑93 (MQ=255)
              aTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAttt                                                                            >  2:1879480/1‑93 (MQ=255)
               tGAGTTCAGCAACTCACGCTGAGTCCGCAGTTGTGCTCCCAGAATACGGTTTTGCCCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTg                                                                           <  2:1018138/93‑1 (MQ=255)
                  gTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTTAGCGGCTGACCATCGGCCTGATGAATTTGCCg                                                                        <  1:3283229/93‑1 (MQ=255)
                    tcatcaACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTa                                                                      <  2:3313634/93‑1 (MQ=255)
                    tCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTa                                                                      <  1:2192791/93‑1 (MQ=255)
                      aGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTagc                                                                    <  1:3952895/93‑1 (MQ=255)
                       gCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTagca                                                                   <  1:3674675/93‑1 (MQ=255)
                       gCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTagca                                                                   >  2:3701412/1‑93 (MQ=255)
                          aCTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCg                                                                >  1:4652638/1‑93 (MQ=255)
                           cTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCgg                                                               >  2:3041260/1‑93 (MQ=255)
                           cTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCgg                                                               >  1:4375676/1‑93 (MQ=255)
                           cTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCgg                                                               >  1:4314465/1‑93 (MQ=255)
                            tCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGc                                                              >  2:1359834/1‑93 (MQ=255)
                             cACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCt                                                             <  2:4061282/93‑1 (MQ=255)
                                gCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGtt                                                          >  2:1492420/1‑93 (MQ=255)
                                 cTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGttt                                                         >  2:4687706/1‑93 (MQ=255)
                                 cTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGttt                                                         <  2:3458773/93‑1 (MQ=255)
                                   taGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTttat                                                       <  2:3651904/92‑1 (MQ=255)
                                    aGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTttatt                                                      <  1:1655271/93‑1 (MQ=255)
                                    aGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTttatt                                                      >  1:1146568/1‑93 (MQ=255)
                                      tGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTaa                                                    <  1:1079913/93‑1 (MQ=255)
                                      tGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTaa                                                    <  1:4193573/93‑1 (MQ=255)
                                       gcgcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAg                                                   <  2:2494830/93‑1 (MQ=255)
                                       gcgcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAg                                                   >  1:2817489/1‑93 (MQ=255)
                                       gcgcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGGTTATTAAg                                                   >  2:203381/1‑93 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  1:1630289/93‑1 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  2:1335080/93‑1 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  1:680261/93‑1 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  2:822492/93‑1 (MQ=255)
                                          cAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAg                                                <  2:4264843/93‑1 (MQ=255)
                                          cAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAg                                                <  2:37237/93‑1 (MQ=255)
                                             ttGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATc                                             >  1:3699739/1‑93 (MQ=255)
                                             ttGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATc                                             >  1:1224811/1‑93 (MQ=255)
                                             ttGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATc                                             <  1:3979022/93‑1 (MQ=255)
                                              tgtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATcc                                            <  1:2842838/93‑1 (MQ=255)
                                               gtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCt                                           <  2:1897949/93‑1 (MQ=255)
                                               gtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCt                                           <  1:468288/93‑1 (MQ=255)
                                               gtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCACATCCt                                           >  2:684249/1‑93 (MQ=255)
                                                   ttCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATa                                       >  1:3486563/1‑93 (MQ=255)
                                                    tCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATaa                                      <  2:2944777/93‑1 (MQ=255)
                                                     ccAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAAc                                     <  2:1431386/93‑1 (MQ=255)
                                                      cAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACg                                    >  1:1797473/1‑93 (MQ=255)
                                                       aGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGt                                   >  1:2542280/1‑93 (MQ=255)
                                                       aGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGt                                   >  1:1994593/1‑93 (MQ=255)
                                                        gAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTa                                  >  1:2405596/1‑93 (MQ=255)
                                                         aaTACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTaa                                 <  1:4155987/93‑1 (MQ=255)
                                                         aaTACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTaa                                 >  2:504973/1‑93 (MQ=255)
                                                          aTACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAc                                <  1:1840370/93‑1 (MQ=255)
                                                          aTACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAc                                <  2:2631402/93‑1 (MQ=255)
                                                          aTACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAc                                <  2:211421/93‑1 (MQ=255)
                                                           tACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAcc                               >  1:3871530/1‑93 (MQ=255)
                                                             cGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGt                             <  2:2193868/93‑1 (MQ=255)
                                                              ggTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGtt                            <  1:215633/93‑1 (MQ=255)
                                                               gTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGGTg                           >  2:3261404/1‑93 (MQ=255)
                                                                ttttGCTCTGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGc                          <  2:1264732/93‑1 (MQ=255)
                                                                ttttGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGACGTAGCAGCGGCTGTTTATTAAGCAGATACTCATAACGTAACCGTTGc                          >  1:3398963/1‑93 (MQ=255)
                                                                   tGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGcaca                       <  2:3394836/93‑1 (MQ=255)
                                                                      tCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGc                    >  2:1626311/1‑93 (MQ=255)
                                                                      tCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGc                    <  2:3763511/93‑1 (MQ=255)
                                                                        ggCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCaa                  >  2:1629668/1‑93 (MQ=255)
                                                                         gCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATCAAGTAACCGTTTCACCCGCAAc                 >  2:3549336/1‑93 (MQ=255)
                                                                         gCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAc                 >  2:2272761/1‑93 (MQ=255)
                                                                         gCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAc                 <  1:2528306/93‑1 (MQ=255)
                                                                          cAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACt                <  1:2635094/93‑1 (MQ=255)
                                                                            gTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgc              <  2:2318150/93‑1 (MQ=255)
                                                                            gTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgc              <  1:4343663/93‑1 (MQ=255)
                                                                            gTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgc              <  2:268468/93‑1 (MQ=255)
                                                                              cAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgcgc            >  2:4402126/1‑93 (MQ=255)
                                                                               aGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGcc           <  1:36193/93‑1 (MQ=255)
                                                                               aGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGcc           <  1:2513614/93‑1 (MQ=255)
                                                                                 cGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAt         <  1:293155/93‑1 (MQ=255)
                                                                                 cGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAt         >  1:3344094/1‑93 (MQ=255)
                                                                                   gCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAttt       >  1:4391576/1‑93 (MQ=255)
                                                                                   gCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAttt       >  2:1078076/1‑93 (MQ=255)
                                                                                    cTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTc      >  2:1468743/1‑93 (MQ=255)
                                                                                     tGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCg     >  1:1269600/1‑93 (MQ=255)
                                                                                      gACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCgg    <  1:1488850/93‑1 (MQ=255)
                                                                                      gACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCgg    <  2:1702087/93‑1 (MQ=255)
                                                                                          |                                                                                          
GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA  >  NC_000913/4385758‑4385938

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ 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
TCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCA  >  NC_000913/4385756‑4385941
                                                                                            |                                                                                             
TCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCAT                                                                                              >  SRR4302535.344834/1‑94 (MQ=60)
 CGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATC                                                                                             >  SRR4302535.2724307/1‑94 (MQ=60)
  GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCG                                                                                            <  SRR4302535.4201016/94‑1 (MQ=60)
  GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCG                                                                                            >  SRR4302535.368756/1‑94 (MQ=60)
  GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCG                                                                                            >  SRR4302535.1524499/1‑94 (MQ=60)
  GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCG                                                                                            <  SRR4302535.4697484/94‑1 (MQ=60)
   CCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGG                                                                                           >  SRR4302535.4633797/1‑94 (MQ=60)
    CACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGC                                                                                          >  SRR4302535.1018068/1‑94 (MQ=60)
    CACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGC                                                                                          <  SRR4302535.3284682/94‑1 (MQ=60)
     ACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCC                                                                                         >  SRR4302535.2833394/1‑94 (MQ=60)
     ACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCC                                                                                         >  SRR4302535.1419852/1‑94 (MQ=60)
      CCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCT                                                                                        <  SRR4302535.1741373/94‑1 (MQ=60)
        CTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGA                                                                                      >  SRR4302535.4076475/1‑94 (MQ=60)
         TGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGAT                                                                                     >  SRR4302535.265313/1‑94 (MQ=60)
         TGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGAT                                                                                     <  SRR4302535.1741299/94‑1 (MQ=60)
         TGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGCCAGCGGCTGACCATCGGCCTGAT                                                                                     >  SRR4302535.4200921/1‑94 (MQ=60)
          GCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATG                                                                                    <  SRR4302535.3141500/94‑1 (MQ=60)
           CAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGA                                                                                   >  SRR4302535.1901015/1‑94 (MQ=60)
           CAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGA                                                                                   >  SRR4302535.4056064/1‑94 (MQ=60)
            AACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTGCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAA                                                                                  >  SRR4302535.3325608/1‑94 (MQ=60)
            AACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAA                                                                                  <  SRR4302535.1862500/94‑1 (MQ=60)
            AACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAA                                                                                  >  SRR4302535.4674323/1‑94 (MQ=60)
             ACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAAT                                                                                 >  SRR4302535.4243459/1‑94 (MQ=60)
             ACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAT                                                                                 <  SRR4302535.2723097/94‑1 (MQ=60)
             ACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGAACATCGGCCTGATGAAT                                                                                 >  SRR4302535.1388721/1‑94 (MQ=60)
              CAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATT                                                                                <  SRR4302535.1361725/94‑1 (MQ=60)
               AATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTT                                                                               <  SRR4302535.4457572/94‑1 (MQ=60)
                 TGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATAAATTTGC                                                                             >  SRR4302535.341660/1‑94 (MQ=60)
                  GAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCGCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCC                                                                            <  SRR4302535.260502/94‑1 (MQ=60)
                   AGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCG                                                                           >  SRR4302535.501436/1‑94 (MQ=60)
                    GTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTTAGCGGCTGACCATCGGCCTGATGAATTTGCCGT                                                                          <  SRR4302535.3306033/94‑1 (MQ=60)
                    GTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGT                                                                          <  SRR4302535.1303889/94‑1 (MQ=60)
                    GTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGT                                                                          >  SRR4302535.3704469/1‑94 (MQ=60)
                      TCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAG                                                                        <  SRR4302535.4306750/94‑1 (MQ=60)
                      TCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAG                                                                        <  SRR4302535.2208418/94‑1 (MQ=60)
                      TCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCGCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAG                                                                        <  SRR4302535.2539773/94‑1 (MQ=60)
                        AGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCA                                                                      <  SRR4302535.3979882/94‑1 (MQ=60)
                         GCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAG                                                                     <  SRR4302535.3700119/94‑1 (MQ=60)
                          CAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGC                                                                    >  SRR4302535.3993476/1‑94 (MQ=60)
                           AACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCG                                                                   <  SRR4302535.2833635/94‑1 (MQ=60)
                           AACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCG                                                                   >  SRR4302535.4684140/1‑94 (MQ=60)
                            ACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGG                                                                  >  SRR4302535.3052363/1‑94 (MQ=60)
                            ACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGG                                                                  >  SRR4302535.2677474/1‑94 (MQ=60)
                            ACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGG                                                                  >  SRR4302535.4343979/1‑94 (MQ=60)
                            ACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGG                                                                  >  SRR4302535.4405615/1‑94 (MQ=60)
                              TCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCT                                                                <  SRR4302535.4680516/94‑1 (MQ=60)
                                  GCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCTGCCTGATGAATTTGCCGTAGCAGCGGCTGTTT                                                            <  SRR4302535.4642166/94‑1 (MQ=60)
                                  GCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTT                                                            <  SRR4302535.110407/94‑1 (MQ=60)
                                  GCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGCTT                                                            >  SRR4302535.3454973/1‑94 (MQ=60)
                                    TGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTAT                                                          <  SRR4302535.1295749/94‑1 (MQ=60)
                                     GAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATT                                                         >  SRR4302535.1154737/1‑94 (MQ=60)
                                      AGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTA                                                        <  SRR4302535.1667243/94‑1 (MQ=60)
                                      AGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTAACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTA                                                        <  SRR4302535.3322997/94‑1 (MQ=60)
                                        TGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAG                                                      <  SRR4302535.1087597/94‑1 (MQ=60)
                                        TGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAG                                                      <  SRR4302535.4222271/94‑1 (MQ=60)
                                        TGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAG                                                      <  SRR4302535.3965236/94‑1 (MQ=60)
                                        TGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAG                                                      >  SRR4302535.2837500/1‑94 (MQ=60)
                                           GCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAG                                                   <  SRR4302535.685027/94‑1 (MQ=60)
                                           GCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAG                                                   <  SRR4302535.1642061/94‑1 (MQ=60)
                                             AGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGAT                                                 >  SRR4302535.1970493/1‑94 (MQ=60)
                                              GTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATC                                                <  SRR4302535.3820185/94‑1 (MQ=60)
                                              GTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATC                                                >  SRR4302535.1233570/1‑94 (MQ=60)
                                              GTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATC                                                >  SRR4302535.3608747/1‑94 (MQ=60)
                                              GTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATC                                                >  SRR4302535.3725366/1‑94 (MQ=60)
                                               TTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCC                                               <  SRR4302535.4006211/94‑1 (MQ=60)
                                               TTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCC                                               >  SRR4302535.975683/1‑94 (MQ=60)
                                                TGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCT                                              <  SRR4302535.2862995/94‑1 (MQ=60)
                                                TGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCT                                              <  SRR4302535.1769693/94‑1 (MQ=60)
                                                 GTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTC                                             >  SRR4302535.3286361/1‑94 (MQ=60)
                                                 GTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTC                                             <  SRR4302535.471562/94‑1 (MQ=60)
                                                  TGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCA                                            <  SRR4302535.4192787/94‑1 (MQ=60)
                                                    CTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATA                                          >  SRR4302535.3510769/1‑94 (MQ=60)
                                                     TTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAA                                         <  SRR4302535.2513320/94‑1 (MQ=60)
                                                      TCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAAC                                        <  SRR4302535.1862692/94‑1 (MQ=60)
                                                       CCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACG                                       >  SRR4302535.1810467/1‑94 (MQ=60)
                                                        CAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGT                                      >  SRR4302535.1066111/1‑94 (MQ=60)
                                                        CAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGT                                      >  SRR4302535.2008958/1‑94 (MQ=60)
                                                        CAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGT                                      >  SRR4302535.1404991/1‑94 (MQ=60)
                                                        CAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGT                                      >  SRR4302535.2560589/1‑94 (MQ=60)
                                                         AGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTA                                     >  SRR4302535.2422908/1‑94 (MQ=60)
                                                           AATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAC                                   <  SRR4302535.4184438/94‑1 (MQ=60)
                                                            ATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACC                                  <  SRR4302535.1853645/94‑1 (MQ=60)
                                                            ATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACC                                  <  SRR4302535.4012762/94‑1 (MQ=60)
                                                            ATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACC                                  <  SRR4302535.2555841/94‑1 (MQ=60)
                                                            ATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACC                                  >  SRR4302535.3898135/1‑94 (MQ=60)
                                                                GGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTG                              >  SRR4302535.4394980/1‑94 (MQ=60)
                                                                GGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTG                              <  SRR4302535.217111/94‑1 (MQ=60)
                                                                 GTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGACGTAGCAGCGGCTGTTTATTAAGCAGATACTCATAACGTAACCGTTGC                             >  SRR4302535.3422613/1‑94 (MQ=60)
                                                                  TTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCA                            >  SRR4302535.3189338/1‑94 (MQ=60)
                                                                  TTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCA                            <  SRR4302535.719661/94‑1 (MQ=60)
                                                                   TTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCAC                           <  SRR4302535.4472296/94‑1 (MQ=60)
                                                                     TGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCGCAC                         <  SRR4302535.3657474/94‑1 (MQ=60)
                                                                      GCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACG                        <  SRR4302535.2846369/94‑1 (MQ=60)
                                                                       CTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCATATCCTCATAACGTAACCGTTGCACACGC                       >  SRR4302535.4365953/1‑94 (MQ=60)
                                                                          GGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAC                    <  SRR4302535.2914057/94‑1 (MQ=60)
                                                                          GGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAC                    <  SRR4302535.4416039/94‑1 (MQ=60)
                                                                          GGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAC                    >  SRR4302535.1839946/1‑94 (MQ=60)
                                                                          GGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAC                    <  SRR4302535.2215253/94‑1 (MQ=60)
                                                                           GCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACT                   <  SRR4302535.2546520/94‑1 (MQ=60)
                                                                           GCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACT                   <  SRR4302535.2064211/94‑1 (MQ=60)
                                                                            CAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTG                  <  SRR4302535.2653977/94‑1 (MQ=60)
                                                                             AGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGC                 >  SRR4302535.4091548/1‑94 (MQ=60)
                                                                              GTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCG                >  SRR4302535.1660683/1‑94 (MQ=60)
                                                                              GTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCG                <  SRR4302535.4373381/94‑1 (MQ=60)
                                                                              GTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCG                <  SRR4302535.1518424/94‑1 (MQ=60)
                                                                                CAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCC              >  SRR4302535.2180957/1‑94 (MQ=60)
                                                                                 AGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCA             <  SRR4302535.2531737/94‑1 (MQ=60)
                                                                                 AGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCA             <  SRR4302535.36472/94‑1 (MQ=60)
                                                                                  GCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAT            >  SRR4302535.3367390/1‑94 (MQ=60)
                                                                                   CGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATT           <  SRR4302535.1278372/94‑1 (MQ=60)
                                                                                   CGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATT           <  SRR4302535.295205/94‑1 (MQ=60)
                                                                                    GGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTT          >  SRR4302535.120720/1‑94 (MQ=60)
                                                                                    GGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTT          >  SRR4302535.4421622/1‑94 (MQ=60)
                                                                                     GCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTC         <  SRR4302535.1373431/94‑1 (MQ=60)
                                                                                      CTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCG        <  SRR4302535.2409870/94‑1 (MQ=60)
                                                                                      CTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCG        >  SRR4302535.1278674/1‑94 (MQ=60)
                                                                                        GACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGT      <  SRR4302535.1499525/94‑1 (MQ=60)
                                                                                        GACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGT      >  SRR4302535.2924070/1‑94 (MQ=60)
                                                                                         ACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA     <  SRR4302535.1465772/94‑1 (MQ=60)
                                                                                         ACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA     >  SRR4302535.4718335/1‑94 (MQ=60)
                                                                                         ACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA     >  SRR4302535.2568345/1‑94 (MQ=60)
                                                                                          CCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTAT    <  SRR4302535.4545295/94‑1 (MQ=60)
                                                                                          CCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTAT    <  SRR4302535.626716/94‑1 (MQ=60)
                                                                                          CCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTAT    <  SRR4302535.3785685/94‑1 (MQ=60)
                                                                                          CCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTAT    <  SRR4302535.2371953/94‑1 (MQ=60)
                                                                                           CATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATC   >  SRR4302535.3735847/1‑94 (MQ=60)
                                                                                            ATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCA  >  SRR4302535.1482613/1‑94 (MQ=60)
                                                                                            ATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCA  >  SRR4302535.1539515/1‑94 (MQ=60)
                                                                                            |                                                                                             
TCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCA  >  NC_000913/4385756‑4385941

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