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 4,001,486 2 bp→GC coding (316‑317/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4860CG99.6% 783.1 / ‑5.6 246R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (114/131);  minor base A (1/0);  total (115/131)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
*NC_0009134,001,4870GC100.0% 797.7 / NA 246R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (115/131);  total (115/131)

Only 100 of 247 total aligned reads displayed.

AGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAA  >  NC_000913/4001399‑4001573
                                                                                       ||                                                                                      
agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                                                    >  2:4797083/1‑93 (MQ=255)
   caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                 >  2:4121137/1‑93 (MQ=255)
    aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                                                >  1:2461903/1‑93 (MQ=255)
      tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                              <  1:5122292/93‑1 (MQ=255)
      tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                              >  2:4711250/1‑93 (MQ=255)
       gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                             <  2:5591458/93‑1 (MQ=255)
         aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                                           >  2:2997475/1‑93 (MQ=255)
         aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                                           >  2:6244454/1‑93 (MQ=255)
         aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                                           <  2:1838664/93‑1 (MQ=255)
          tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                                          >  2:1177915/1‑93 (MQ=255)
          tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                                          <  1:521755/93‑1 (MQ=255)
             ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt                                                                       >  1:201358/1‑93 (MQ=255)
             ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt                                                                       >  2:5614452/1‑93 (MQ=255)
              gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                                      <  2:4195643/93‑1 (MQ=255)
              gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                                      >  2:3797305/1‑93 (MQ=255)
              gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                                      <  2:2674966/93‑1 (MQ=255)
              gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                                      >  2:3346880/1‑93 (MQ=255)
               cggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTg                                                                     <  1:5816503/93‑1 (MQ=255)
                 gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt                                                                   <  2:5879062/93‑1 (MQ=255)
                 gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt                                                                   <  1:3494961/93‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                >  1:4269956/5‑93 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                >  1:5419174/5‑93 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                <  1:63059/89‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                <  1:3181138/89‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                >  2:1341771/5‑93 (MQ=255)
                      gctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCg                                                              >  1:2165657/3‑93 (MQ=255)
                        tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                            >  2:3875953/1‑93 (MQ=255)
                        tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                            >  1:4725252/1‑93 (MQ=255)
                         cGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                                           >  2:1869142/1‑93 (MQ=255)
                          gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                                          <  2:3626685/93‑1 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                        <  2:711710/93‑1 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                        >  1:1263059/1‑93 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                        <  2:1557505/93‑1 (MQ=255)
                             gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                                       <  2:1430178/93‑1 (MQ=255)
                               ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                     >  2:2355775/1‑93 (MQ=255)
                               ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                     >  1:137959/1‑93 (MQ=255)
                                gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                    >  1:4408835/1‑93 (MQ=255)
                                gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                    <  1:5564534/93‑1 (MQ=255)
                                gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                    >  1:132379/1‑93 (MQ=255)
                                gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                    >  2:397547/1‑93 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                  >  2:751075/1‑93 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                  >  2:3182348/1‑93 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                  <  2:5727846/93‑1 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                  <  2:3136766/93‑1 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                  <  1:3323090/93‑1 (MQ=255)
                                     tcatcaCCCGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                                               >  2:6107558/1‑93 (MQ=255)
                                     tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                                               <  1:6165799/93‑1 (MQ=255)
                                     tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                                               <  2:5316604/93‑1 (MQ=255)
                                      catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                                              <  1:5964779/93‑1 (MQ=255)
                                      catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                                              >  1:2785376/1‑93 (MQ=255)
                                      catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                                              <  2:2492498/93‑1 (MQ=255)
                                       atcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGc                                             <  2:4943946/93‑1 (MQ=255)
                                        tcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                                            >  1:1601012/1‑93 (MQ=255)
                                         ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCa                                           <  1:4555240/93‑1 (MQ=255)
                                          cccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                                          >  1:2532751/1‑93 (MQ=255)
                                          cccAGCACAATGGTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                                          >  1:3405366/1‑93 (MQ=255)
                                           ccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGt                                         >  1:3981726/1‑93 (MQ=255)
                                              gCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                                      >  1:6221888/1‑93 (MQ=255)
                                                acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                    >  2:2111319/1‑93 (MQ=255)
                                                 caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                                   >  1:1095514/1‑93 (MQ=255)
                                                  aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                  >  2:5148629/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                 >  2:1640144/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                 >  1:2050210/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                 >  1:2605303/1‑93 (MQ=255)
                                                    tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCaa                                <  1:4793726/93‑1 (MQ=255)
                                                      ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc                              <  1:4129488/93‑1 (MQ=255)
                                                      ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc                              <  2:2802367/93‑1 (MQ=255)
                                                       tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                             >  2:103085/1‑93 (MQ=255)
                                                       tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                             >  2:4399208/1‑93 (MQ=255)
                                                        cACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAg                            <  1:3834003/93‑1 (MQ=255)
                                                        cACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAg                            <  2:2255416/93‑1 (MQ=255)
                                                         aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                           >  2:535675/1‑93 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                          >  2:1241072/1‑93 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                          >  2:5483429/1‑93 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                          <  2:1363512/93‑1 (MQ=255)
                                                           cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGTCAGTTGCCACCAATCAGCAc                         <  2:3586794/93‑1 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                     <  1:4606462/93‑1 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                     >  1:2321566/1‑93 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATTAGCACGGCag                    <  2:1479494/93‑1 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag                    <  1:4960504/93‑1 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                  >  1:5875562/1‑93 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                  <  1:5951955/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                <  1:5470870/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                <  1:3893496/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                <  1:191243/93‑1 (MQ=255)
                                                                      taaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTg              >  1:4618660/1‑93 (MQ=255)
                                                                       aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCGGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGc             <  1:1012866/93‑1 (MQ=255)
                                                                       aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGc             <  1:2977841/93‑1 (MQ=255)
                                                                        aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGcc            <  1:1836569/93‑1 (MQ=255)
                                                                         aGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCa           <  2:6224589/93‑1 (MQ=255)
                                                                          gTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg          <  1:770172/93‑1 (MQ=255)
                                                                           tAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGc         <  2:6134593/93‑1 (MQ=255)
                                                                            aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCa        <  1:856062/93‑1 (MQ=255)
                                                                             aCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCAt       <  1:3094129/93‑1 (MQ=255)
                                                                               cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaa     <  1:4678224/93‑1 (MQ=255)
                                                                               cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaa     >  1:965992/1‑93 (MQ=255)
                                                                                aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACCGCAGAGACTGCCAGCATaaa    >  1:2875549/1‑93 (MQ=255)
                                                                                 aGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaa   <  1:1451390/93‑1 (MQ=255)
                                                                                  ggCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaaa  >  1:2842309/1‑93 (MQ=255)
                                                                                  ggCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaaa  <  2:2885613/93‑1 (MQ=255)
                                                                                       ||                                                                                      
AGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAA  >  NC_000913/4001399‑4001573

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATT  >  NC_000913/4001393‑4001575
                                                                                             |                                                                                         
CCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCG                                                                                           <  SRR4302533.1692261/94‑1 (MQ=60)
     CAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCC                                                                                      >  SRR4302533.856190/1‑94 (MQ=60)
      AGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCA                                                                                     >  SRR4302533.6268655/1‑94 (MQ=60)
       GCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAG                                                                                    <  SRR4302533.5162980/94‑1 (MQ=60)
         CATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCA                                                                                  >  SRR4302533.2477092/1‑94 (MQ=60)
          ATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCAT                                                                                 <  SRR4302533.1958512/94‑1 (MQ=60)
           TTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATA                                                                                <  SRR4302533.4643289/94‑1 (MQ=60)
            TGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATAT                                                                               <  SRR4302533.5153239/94‑1 (MQ=60)
              CATCCGGCGGAATCGCTTGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGG                                                                             >  SRR4302533.3484628/1‑94 (MQ=60)
              CATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGG                                                                             <  SRR4302533.4256834/94‑1 (MQ=60)
               ATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGT                                                                            >  SRR4302533.4393085/1‑94 (MQ=60)
               ATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGT                                                                            <  SRR4302533.1890037/94‑1 (MQ=60)
                TCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTG                                                                           <  SRR4302533.524851/94‑1 (MQ=60)
                 CCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGA                                                                          <  SRR4302533.5922435/94‑1 (MQ=60)
                  CGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACTAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAT                                                                         <  SRR4302533.4115321/94‑1 (MQ=60)
                  CGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAT                                                                         <  SRR4302533.3478797/94‑1 (MQ=60)
                  CGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAT                                                                         >  SRR4302533.202638/1‑94 (MQ=60)
                   GGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATT                                                                        <  SRR4302533.3020907/94‑1 (MQ=60)
                   GGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATT                                                                        <  SRR4302533.2874708/94‑1 (MQ=60)
                   GGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATT                                                                        >  SRR4302533.6175630/1‑94 (MQ=60)
                    GCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTG                                                                       <  SRR4302533.2431461/94‑1 (MQ=60)
                     CGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGT                                                                      >  SRR4302533.1398535/1‑94 (MQ=60)
                     CGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGT                                                                      <  SRR4302533.5851633/94‑1 (MQ=60)
                     CGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATAGTGATTGT                                                                      >  SRR4302533.1312277/1‑94 (MQ=60)
                      GGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTT                                                                     >  SRR4302533.2424617/1‑94 (MQ=60)
                       GAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTC                                                                    <  SRR4302533.3516515/94‑1 (MQ=60)
                       GAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTC                                                                    <  SRR4302533.2033885/94‑1 (MQ=60)
                       GAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTC                                                                    <  SRR4302533.4574311/94‑1 (MQ=60)
                        AATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCA                                                                   >  SRR4302533.1314569/1‑94 (MQ=60)
                         atcgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAC                                                                  >  SRR4302533.4295659/6‑94 (MQ=60)
                         atcgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAC                                                                  >  SRR4302533.5451775/6‑94 (MQ=60)
                         atcgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAC                                                                  <  SRR4302533.1171701/89‑1 (MQ=60)
                          tcgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACC                                                                 <  SRR4302533.63473/90‑1 (MQ=60)
                          tcgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACC                                                                 <  SRR4302533.3200709/90‑1 (MQ=60)
                           cgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCG                                                                <  SRR4302533.2584346/91‑1 (MQ=60)
                           cgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCG                                                                <  SRR4302533.5895378/91‑1 (MQ=60)
                           cgcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCG                                                                >  SRR4302533.2178979/4‑94 (MQ=60)
                            gcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGC                                                               <  SRR4302533.2920479/92‑1 (MQ=60)
                            gcTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGC                                                               >  SRR4302533.166363/3‑94 (MQ=60)
                             CTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCC                                                              >  SRR4302533.4753902/1‑94 (MQ=60)
                             CTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCC                                                              <  SRR4302533.1414040/94‑1 (MQ=60)
                               CGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCA                                                            <  SRR4302533.1091985/94‑1 (MQ=60)
                                 CCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGA                                                          >  SRR4302533.3891838/1‑94 (MQ=60)
                                 CCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGA                                                          >  SRR4302533.1270615/1‑94 (MQ=60)
                                 CCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGA                                                          >  SRR4302533.4691613/1‑94 (MQ=60)
                                  CGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAT                                                         >  SRR4302533.3957490/1‑94 (MQ=60)
                                  CGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAT                                                         <  SRR4302533.5985666/94‑1 (MQ=60)
                                  CGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAT                                                         >  SRR4302533.5187593/1‑94 (MQ=60)
                                   GAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATA                                                        >  SRR4302533.3774819/1‑94 (MQ=60)
                                    AGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAA                                                       >  SRR4302533.138834/1‑94 (MQ=60)
                                     GGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAA                                                      >  SRR4302533.4435460/1‑94 (MQ=60)
                                     GGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAA                                                      >  SRR4302533.133211/1‑94 (MQ=60)
                                      GAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAAC                                                     <  SRR4302533.5598169/94‑1 (MQ=60)
                                        AAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAG                                                   <  SRR4302533.3343613/94‑1 (MQ=60)
                                        AAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCACAGATAAAAAG                                                   >  SRR4302533.2584432/1‑94 (MQ=60)
                                           TCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAG                                                >  SRR4302533.2802491/1‑94 (MQ=60)
                                           TCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAG                                                <  SRR4302533.4321966/94‑1 (MQ=60)
                                           TCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAG                                                <  SRR4302533.6202860/94‑1 (MQ=60)
                                           TCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAG                                                <  SRR4302533.4751128/94‑1 (MQ=60)
                                            CATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTTGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGC                                               >  SRR4302533.6249153/1‑94 (MQ=60)
                                            CATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGC                                               <  SRR4302533.6000875/94‑1 (MQ=60)
                                             ATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCC                                              >  SRR4302533.481303/1‑94 (MQ=60)
                                             ATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCC                                              >  SRR4302533.1610702/1‑94 (MQ=60)
                                             ATCCCCAGCACAATGTTCACCAGCGGGGTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGAC                                              >  SRR4302533.5664045/1‑94 (MQ=60)
                                               CCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAG                                            <  SRR4302533.2190096/94‑1 (MQ=60)
                                               CCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAG                                            <  SRR4302533.4582868/94‑1 (MQ=60)
                                               CCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAG                                            >  SRR4302533.2548342/1‑94 (MQ=60)
                                               CCCCAGCACAATGGTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAG                                            >  SRR4302533.3426398/1‑94 (MQ=60)
                                                CCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGT                                           <  SRR4302533.3251713/94‑1 (MQ=60)
                                                CCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGT                                           >  SRR4302533.4005968/1‑94 (MQ=60)
                                                CCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGT                                           >  SRR4302533.2829554/1‑94 (MQ=60)
                                                CCAAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGT                                           >  SRR4302533.5221819/1‑94 (MQ=60)
                                                  CAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTG                                         >  SRR4302533.4958347/1‑94 (MQ=60)
                                                  CAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTG                                         <  SRR4302533.4748143/94‑1 (MQ=60)
                                                   AGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGC                                        >  SRR4302533.6259180/1‑94 (MQ=60)
                                                   AGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGC                                        >  SRR4302533.3131243/1‑94 (MQ=60)
                                                     CACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCA                                      >  SRR4302533.1610954/1‑94 (MQ=60)
                                                     CACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCA                                      >  SRR4302533.2104217/1‑94 (MQ=60)
                                                      ACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCAC                                     >  SRR4302533.1102000/1‑94 (MQ=60)
                                                        AATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCA                                   >  SRR4302533.4033358/1‑94 (MQ=60)
                                                        AATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCA                                   >  SRR4302533.2621306/1‑94 (MQ=60)
                                                        AATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCA                                   >  SRR4302533.2062859/1‑94 (MQ=60)
                                                         ATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAA                                  <  SRR4302533.4614852/94‑1 (MQ=60)
                                                          TGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAT                                 <  SRR4302533.4822767/94‑1 (MQ=60)
                                                            TTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCA                               <  SRR4302533.6272106/94‑1 (MQ=60)
                                                            TTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCA                               <  SRR4302533.3236299/94‑1 (MQ=60)
                                                            TTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCA                               <  SRR4302533.3346612/94‑1 (MQ=60)
                                                            TTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCA                               <  SRR4302533.4154555/94‑1 (MQ=60)
                                                             TCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAG                              <  SRR4302533.3007496/94‑1 (MQ=60)
                                                              CACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGC                             <  SRR4302533.3857434/94‑1 (MQ=60)
                                                              CACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGC                             >  SRR4302533.5992457/1‑94 (MQ=60)
                                                               ACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCA                            <  SRR4302533.2128170/94‑1 (MQ=60)
                                                                 CAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACG                          >  SRR4302533.2426850/1‑94 (MQ=60)
                                                                 CAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACG                          <  SRR4302533.1550892/94‑1 (MQ=60)
                                                                    CGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCA                       >  SRR4302533.2335935/1‑94 (MQ=60)
                                                                     GGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAG                      <  SRR4302533.6180328/94‑1 (MQ=60)
                                                                     GGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAG                      <  SRR4302533.869259/94‑1 (MQ=60)
                                                                     GGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAG                      <  SRR4302533.4634365/94‑1 (MQ=60)
                                                                     GGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAG                      <  SRR4302533.2892216/94‑1 (MQ=60)
                                                                      GGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGA                     <  SRR4302533.1004550/94‑1 (MQ=60)
                                                                      GGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGA                     <  SRR4302533.4990583/94‑1 (MQ=60)
                                                                      GGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGA                     <  SRR4302533.773479/94‑1 (MQ=60)
                                                                       GTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAG                    <  SRR4302533.3885994/94‑1 (MQ=60)
                                                                       GTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAG                    <  SRR4302533.5367977/94‑1 (MQ=60)
                                                                       GTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAG                    <  SRR4302533.2659628/94‑1 (MQ=60)
                                                                       GTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAG                    >  SRR4302533.5911081/1‑94 (MQ=60)
                                                                       GTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAG                    <  SRR4302533.5714931/94‑1 (MQ=60)
                                                                        TTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGA                   <  SRR4302533.5987959/94‑1 (MQ=60)
                                                                          AATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACT                 <  SRR4302533.3917256/94‑1 (MQ=60)
                                                                          AATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACT                 <  SRR4302533.192444/94‑1 (MQ=60)
                                                                          AATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACT                 <  SRR4302533.5503841/94‑1 (MQ=60)
                                                                          AATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACT                 >  SRR4302533.131167/1‑94 (MQ=60)
                                                                           ATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTG                >  SRR4302533.4646649/1‑94 (MQ=60)
                                                                           ATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTG                <  SRR4302533.1160109/94‑1 (MQ=60)
                                                                             AAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCGGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCC              <  SRR4302533.1018796/94‑1 (MQ=60)
                                                                             AAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCC              >  SRR4302533.4617900/1‑94 (MQ=60)
                                                                             AAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCC              <  SRR4302533.2996060/94‑1 (MQ=60)
                                                                              AAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCA             <  SRR4302533.1847771/94‑1 (MQ=60)
                                                                               AGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAG            >  SRR4302533.1346255/1‑94 (MQ=60)
                                                                               AGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAG            <  SRR4302533.4888485/94‑1 (MQ=60)
                                                                               AGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAG            <  SRR4302533.2002773/94‑1 (MQ=60)
                                                                                GTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGC           <  SRR4302533.774705/94‑1 (MQ=60)
                                                                                GTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGC           >  SRR4302533.2540559/1‑94 (MQ=60)
                                                                                  AACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCAT         >  SRR4302533.5504279/1‑94 (MQ=60)
                                                                                  AACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCAT         >  SRR4302533.5468232/1‑94 (MQ=60)
                                                                                  AACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCAT         <  SRR4302533.861077/94‑1 (MQ=60)
                                                                                   ACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATA        <  SRR4302533.3113152/94‑1 (MQ=60)
                                                                                    CCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAA       >  SRR4302533.971660/1‑94 (MQ=60)
                                                                                     CAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAA      <  SRR4302533.4706575/94‑1 (MQ=60)
                                                                                     CAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAA      <  SRR4302533.5452139/94‑1 (MQ=60)
                                                                                     CAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACCGCAGAGACTGCCAGCATAAA      >  SRR4302533.2893243/1‑94 (MQ=60)
                                                                                      AAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAA     <  SRR4302533.5121338/94‑1 (MQ=60)
                                                                                      AAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAA     >  SRR4302533.1318057/1‑94 (MQ=60)
                                                                                       AGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAA    >  SRR4302533.2859804/1‑94 (MQ=60)
                                                                                       AGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAA    <  SRR4302533.1460148/94‑1 (MQ=60)
                                                                                       AGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAA    >  SRR4302533.5549422/1‑94 (MQ=60)
                                                                                        GGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAT   >  SRR4302533.4531322/1‑94 (MQ=60)
                                                                                        GGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAT   <  SRR4302533.864487/94‑1 (MQ=60)
                                                                                         GCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATT  >  SRR4302533.6064672/1‑94 (MQ=60)
                                                                                             |                                                                                         
ACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATT  >  NC_000913/4001393‑4001575

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