Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,766,427 2 bp→AC coding (697‑698/864 nt) b2632 → putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,766,4270TA100.0% 618.0 / NA 171W233T (TGG→ACG) b2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (75/96);  total (75/96)
*NC_0009132,766,4280GC99.4% 589.2 / ‑5.5 171W233T (TGG→ACG) b2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (75/95);  minor base T (0/1);  total (75/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 175 total aligned reads displayed.

AGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTC  >  NC_000913/2766340‑2766511
                                                                                       ||                                                                                   
aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 >  2:253105/1‑93 (MQ=21)
aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 >  2:2498837/1‑93 (MQ=21)
aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 <  2:1339623/93‑1 (MQ=25)
aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 <  2:787267/93‑1 (MQ=25)
   gTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                                              >  2:3814565/1‑93 (MQ=21)
     gATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCg                                                                            <  2:5243566/93‑1 (MQ=25)
            gagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAg                                                                     <  2:1942291/93‑1 (MQ=25)
              gCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGc                                                                   <  2:928674/93‑1 (MQ=25)
               cGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGcc                                                                  >  2:135671/1‑92 (MQ=25)
                ggATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                                                 >  2:1021563/1‑93 (MQ=255)
                 gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  1:1574544/1‑93 (MQ=255)
                 gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                <  1:322590/93‑1 (MQ=255)
                  atgatTAAGTGTCTGGCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               <  2:4362507/93‑1 (MQ=255)
                  atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               <  1:4036375/93‑1 (MQ=255)
                  atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               >  1:6048635/1‑93 (MQ=255)
                  atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               <  2:317288/93‑1 (MQ=255)
                  atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               <  2:75814/93‑1 (MQ=255)
                  atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               >  1:3604980/1‑93 (MQ=255)
                   tgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTg                                                              >  2:4045824/1‑93 (MQ=255)
                    gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGa                                                             >  1:3018661/1‑93 (MQ=255)
                    gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGa                                                             >  1:6130032/1‑93 (MQ=255)
                    gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGa                                                             <  2:2456460/93‑1 (MQ=255)
                     atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCGCGTGAc                                                            <  1:938479/93‑1 (MQ=255)
                     atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAc                                                            <  1:719481/93‑1 (MQ=255)
                     atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAc                                                            >  2:162256/1‑93 (MQ=255)
                     atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAc                                                            <  2:2496036/93‑1 (MQ=255)
                       tAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGa                                                          <  1:5179828/93‑1 (MQ=255)
                       tAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGa                                                          >  1:2364494/1‑93 (MQ=255)
                        aaGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAt                                                         >  1:5846997/1‑93 (MQ=255)
                           tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                                      <  2:468357/93‑1 (MQ=255)
                           tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                                      <  1:5960550/93‑1 (MQ=255)
                           tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                                      >  1:1681329/1‑93 (MQ=255)
                           tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                                      <  2:3881887/93‑1 (MQ=255)
                           tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                                      <  1:3637485/93‑1 (MQ=255)
                             tCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCgg                                                    >  2:601940/1‑93 (MQ=255)
                             tCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCgg                                                    <  2:2402378/93‑1 (MQ=255)
                              cTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                                                   <  1:5098849/93‑1 (MQ=255)
                               tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTg                                                  <  1:4063450/93‑1 (MQ=255)
                               tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTg                                                  >  2:2029277/1‑93 (MQ=255)
                               tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTAGGTg                                                  >  2:459024/1‑93 (MQ=255)
                               tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGAAGCCCGGGAGCAGGCACGTGACGATTTCGGta                                                  >  2:5005338/1‑92 (MQ=255)
                                gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 >  1:1607/1‑93 (MQ=255)
                                     gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCg                                            >  1:2365576/1‑93 (MQ=255)
                                      tGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCga                                           >  1:77101/1‑92 (MQ=255)
                                      tGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCga                                           >  1:1321641/1‑92 (MQ=255)
                                       gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  2:1801435/93‑1 (MQ=255)
                                       gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  1:533644/93‑1 (MQ=255)
                                       gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  1:678932/93‑1 (MQ=255)
                                       gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  1:1065204/93‑1 (MQ=255)
                                        aggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAg                                         <  2:2789843/90‑1 (MQ=255)
                                        aggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAg                                         <  2:29492/90‑1 (MQ=255)
                                        aggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAg                                         <  2:121128/90‑1 (MQ=255)
                                         ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                                        <  1:3908537/91‑2 (MQ=255)
                                          gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                       <  2:4002648/92‑3 (MQ=255)
                                              gcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCa                                   >  2:3173815/1‑93 (MQ=255)
                                               cagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                  <  2:3372806/93‑1 (MQ=255)
                                               cagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                  <  1:288997/93‑1 (MQ=255)
                                                cgcCCGGTAGGAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATa                                 <  1:5269436/92‑1 (MQ=255)
                                                 gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                >  2:290327/1‑93 (MQ=255)
                                                 gcCCGGTAGGAGCGCTTCTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                <  2:1509657/93‑1 (MQ=255)
                                                   ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                              <  1:4243038/93‑1 (MQ=255)
                                                   ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                              >  1:5705598/1‑93 (MQ=255)
                                                   ccGGTAGTAGCGCTAATACAACATCCCTTTCGAACAACGGAAGCCCGAGAGCAGGCACGTGACGATTTCGGTGAAACCGTAAAAACAATACTc                              <  1:5493557/93‑1 (MQ=255)
                                                    cGGTAGTAGCGCTACTGCAACATTCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                             >  1:4751612/1‑93 (MQ=255)
                                                    cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                             >  1:1907847/1‑93 (MQ=255)
                                                       tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                          <  1:2941609/93‑1 (MQ=255)
                                                       tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                          <  2:2556315/93‑1 (MQ=255)
                                                        agtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                         <  2:3921130/93‑1 (MQ=255)
                                                        agtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                         <  1:4529823/93‑1 (MQ=255)
                                                         gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATaag                        >  2:4489866/1‑91 (MQ=255)
                                                         gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                        >  2:2990869/1‑93 (MQ=255)
                                                         gtagCGCTACTGCAACATCACTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                        <  2:5457595/93‑1 (MQ=255)
                                                          tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg                       <  1:1879839/93‑1 (MQ=255)
                                                           agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      <  1:2206692/93‑1 (MQ=255)
                                                           agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      <  2:2070126/93‑1 (MQ=255)
                                                            gcgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTa                     <  1:5214098/93‑1 (MQ=255)
                                                            gcgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTa                     >  1:4339456/1‑93 (MQ=255)
                                                              gcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAg                   <  1:4728745/93‑1 (MQ=255)
                                                               cttcTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACATGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                  <  1:5390741/90‑1 (MQ=255)
                                                               cTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                  >  1:1470620/1‑93 (MQ=255)
                                                               cTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                  <  1:5598504/93‑1 (MQ=255)
                                                               cTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                  >  1:2406166/1‑93 (MQ=255)
                                                               cTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                  >  2:741122/1‑93 (MQ=255)
                                                                tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAgac                 >  2:3738784/1‑91 (MQ=255)
                                                                caCTGCAACAACCCTTTCGAACAACGGTAGCCCGGGAGGAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCa                 <  2:5211891/92‑1 (MQ=255)
                                                                  cTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAcc               <  1:1787159/93‑1 (MQ=255)
                                                                   tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCt              <  2:1291658/93‑1 (MQ=255)
                                                                    gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             <  1:4248524/93‑1 (MQ=255)
                                                                    gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATGTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             >  1:5475248/1‑93 (MQ=255)
                                                                      aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTc           >  1:3088522/1‑93 (MQ=255)
                                                                      aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTc           >  1:4039005/1‑93 (MQ=255)
                                                                      aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTc           <  2:1439473/93‑1 (MQ=255)
                                                                         aTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc        <  1:1958718/93‑1 (MQ=255)
                                                                          tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCt       <  1:742845/93‑1 (MQ=255)
                                                                          tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCt       <  1:376346/93‑1 (MQ=255)
                                                                           cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt      >  2:4548087/1‑93 (MQ=255)
                                                                           cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt      <  2:1704827/93‑1 (MQ=255)
                                                                           cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt      >  1:1168427/1‑93 (MQ=255)
                                                                            ccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTa     <  2:3614490/93‑1 (MQ=255)
                                                                               ttCGAACAATGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTc  <  1:4163377/93‑1 (MQ=255)
                                                                                       ||                                                                                   
AGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTC  >  NC_000913/2766340‑2766511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
AGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTT  >  NC_000913/2766340‑2766507
                                                                                       |                                                                                
AGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAG                                                                            <  SRR4302532.3795288/94‑1 (MQ=60)
 GGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGC                                                                           >  SRR4302532.4422915/1‑94 (MQ=60)
    TGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCG                                                                        >  SRR4302532.5929251/1‑94 (MQ=60)
       TGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGA                                                                     <  SRR4302532.3271438/94‑1 (MQ=60)
             AGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGC                                                               <  SRR4302532.4197782/94‑1 (MQ=60)
              GCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCA                                                              >  SRR4302532.15968/1‑94 (MQ=60)
                GGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACG                                                            <  SRR4302532.2333287/94‑1 (MQ=60)
                GGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACG                                                            >  SRR4302532.1627845/1‑94 (MQ=60)
                 GATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGT                                                           >  SRR4302532.6257496/1‑94 (MQ=60)
                 GATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGT                                                           <  SRR4302532.333799/94‑1 (MQ=60)
                 GATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGT                                                           >  SRR4302532.3725474/1‑94 (MQ=60)
                  ATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTG                                                          <  SRR4302532.4172344/94‑1 (MQ=60)
                   TGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGA                                                         >  SRR4302532.3118253/1‑94 (MQ=60)
                   TGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGA                                                         >  SRR4302532.6341549/1‑94 (MQ=60)
                     ATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCGCGTGACG                                                       <  SRR4302532.970643/94‑1 (MQ=60)
                     ATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACG                                                       <  SRR4302532.744042/94‑1 (MQ=60)
                      TTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGA                                                      >  SRR4302532.2443367/1‑94 (MQ=60)
                       TAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAT                                                     <  SRR4302532.5358224/94‑1 (MQ=60)
                       TAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAT                                                     >  SRR4302532.6049889/1‑94 (MQ=60)
                        AAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATT                                                    >  SRR4302532.4567625/1‑94 (MQ=60)
                          GTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTC                                                  >  SRR4302532.1738153/1‑94 (MQ=60)
                           TGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCG                                                 <  SRR4302532.3759521/94‑1 (MQ=60)
                           TGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCG                                                 <  SRR4302532.6166684/94‑1 (MQ=60)
                              CTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTG                                              <  SRR4302532.5274161/94‑1 (MQ=60)
                               TGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGA                                             <  SRR4302532.4200001/94‑1 (MQ=60)
                               TGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGA                                             >  SRR4302532.1667/1‑94 (MQ=60)
                                GCCGCGTGAGTCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAA                                            <  SRR4302532.395483/94‑1 (MQ=60)
                                GCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAA                                            >  SRR4302532.2472574/1‑94 (MQ=60)
                                 CCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAA                                           >  SRR4302532.3413058/1‑94 (MQ=60)
                                   GCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACC                                         <  SRR4302532.775211/94‑1 (MQ=60)
                                    CGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCG                                        >  SRR4302532.2444482/1‑94 (MQ=60)
                                    CGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCG                                        >  SRR4302532.5383170/1‑94 (MQ=60)
                                     GTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGA                                       >  SRR4302532.79852/1‑94 (MQ=60)
                                     GTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGA                                       >  SRR4302532.1365810/1‑94 (MQ=60)
                                       gaggctatCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGTGAGCAGTCACGTGACGATATCGGTGGAACCGTAG                                     >  SRR4302532.5677174/9‑94 (MQ=60)
                                       GAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAG                                     <  SRR4302532.1101034/94‑1 (MQ=60)
                                       GAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAG                                     <  SRR4302532.878777/94‑1 (MQ=60)
                                       GAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAG                                     <  SRR4302532.552246/94‑1 (MQ=60)
                                       GAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAG                                     <  SRR4302532.702284/94‑1 (MQ=60)
                                        AGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGG                                    <  SRR4302532.5543346/94‑1 (MQ=60)
                                         GGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGT                                   <  SRR4302532.4040331/94‑1 (MQ=60)
                                         GGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGT                                   <  SRR4302532.4049574/94‑1 (MQ=60)
                                         GGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGT                                   <  SRR4302532.1908130/94‑1 (MQ=60)
                                          GCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGgtg                                  <  SRR4302532.5847794/94‑4 (MQ=60)
                                              GCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAT                              <  SRR4302532.1066291/94‑1 (MQ=60)
                                               CAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATA                             <  SRR4302532.298965/94‑1 (MQ=60)
                                                CGCCCGGTAGGAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAC                            <  SRR4302532.5451136/94‑1 (MQ=60)
                                                AGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAC                            >  SRR4302532.4260095/1‑94 (MQ=60)
                                                  CCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTC                          >  SRR4302532.5904530/1‑94 (MQ=60)
                                                   CCGGTAGTAGCGCTACTGCAACATTCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCG                         >  SRR4302532.4914073/1‑94 (MQ=60)
                                                   CCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCG                         >  SRR4302532.1971164/1‑94 (MQ=60)
                                                   CCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCG                         <  SRR4302532.4385590/94‑1 (MQ=60)
                                                   CCGGTAGTAGCGCTAATACAACATCCCTTTCGAACAACGGAAGCCCGAGAGCAGGCACGTGACGATTTCGGTGAAACCGTAAAAACAATACTCG                         <  SRR4302532.5684058/94‑1 (MQ=46)
                                                   CCGGTAGTAGAGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCG                         >  SRR4302532.1752700/1‑94 (MQ=60)
                                                    CGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGA                        <  SRR4302532.4312797/94‑1 (MQ=60)
                                                       TAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAC                     <  SRR4302532.3039281/94‑1 (MQ=60)
                                                       TAGTAGCGCTACGGCAACATCCCGTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAC                     <  SRR4302532.200195/94‑1 (MQ=60)
                                                        AGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACG                    <  SRR4302532.4684419/94‑1 (MQ=60)
                                                        AGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACG                    <  SRR4302532.1499171/94‑1 (MQ=60)
                                                         GTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGG                   >  SRR4302532.2027963/1‑94 (MQ=60)
                                                          TAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGT                  <  SRR4302532.1942142/94‑1 (MQ=60)
                                                           AGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTA                 >  SRR4302532.4486033/1‑94 (MQ=60)
                                                           AGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTA                 <  SRR4302532.2280061/94‑1 (MQ=60)
                                                           AGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTA                 <  SRR4302532.617394/94‑1 (MQ=60)
                                                           AGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTA                 <  SRR4302532.3136599/94‑1 (MQ=60)
                                                            GCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAA                <  SRR4302532.5393910/94‑1 (MQ=60)
                                                             CGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAG               <  SRR4302532.767207/94‑1 (MQ=60)
                                                              GCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGC              >  SRR4302532.1520170/1‑94 (MQ=60)
                                                              GCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGC              >  SRR4302532.2486240/1‑94 (MQ=60)
                                                              GCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGC              >  SRR4302532.171173/1‑94 (MQ=60)
                                                              GCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGC              <  SRR4302532.4890142/94‑1 (MQ=60)
                                                              GCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCCGTGAAACCGTctgtgacatacacattacggtacga              >  SRR4302532.171174/1‑69 (MQ=60)
                                                               CTTCTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACATGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCA             <  SRR4302532.5577498/94‑1 (MQ=60)
                                                               CTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCC             >  SRR4302532.5013908/1‑94 (MQ=60)
                                                               CTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCA             <  SRR4302532.2989146/94‑1 (MQ=60)
                                                               CTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCA             <  SRR4302532.5792850/94‑1 (MQ=60)
                                                                TACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAC            <  SRR4302532.1315441/94‑1 (MQ=60)
                                                                  CTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGGCATACTCGATACGGTAAGCAACT          >  SRR4302532.2618278/1‑94 (MQ=60)
                                                                  CTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCT          <  SRR4302532.1846892/94‑1 (MQ=60)
                                                                   TGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATGTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTT         >  SRR4302532.5665084/1‑94 (MQ=60)
                                                                   GGCAACATCCCTTTCGAACAACGGGAGCGCGGGAGCAGGCATGTGCCGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAATCACCTT         <  SRR4302532.5578733/94‑1 (MQ=40)
                                                                    GCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTT        <  SRR4302532.4391270/94‑1 (MQ=60)
                                                                     CAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTC       >  SRR4302532.5734464/1‑94 (MQ=60)
                                                                     CAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTC       >  SRR4302532.4175062/1‑94 (MQ=60)
                                                                     CAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTC       >  SRR4302532.3190819/1‑94 (MQ=60)
                                                                      AACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCC      >  SRR4302532.529024/1‑94 (MQ=60)
                                                                        CATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCC    >  SRR4302532.4882309/1‑94 (MQ=60)
                                                                        CATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCC    >  SRR4302532.6359876/1‑94 (MQ=60)
                                                                         ATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCT   >  SRR4302532.6323497/1‑94 (MQ=60)
                                                                         ATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCT   >  SRR4302532.5521231/1‑94 (MQ=60)
                                                                         ATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCT   <  SRR4302532.2023787/94‑1 (MQ=60)
                                                                          TCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTT  >  SRR4302532.1207465/1‑94 (MQ=60)
                                                                          TCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTT  <  SRR4302532.389418/94‑1 (MQ=60)
                                                                          TCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTT  <  SRR4302532.768230/94‑1 (MQ=60)
                                                                                       |                                                                                
AGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTT  >  NC_000913/2766340‑2766507

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