Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,832 A→G R146R (CGT→CGC yjjT ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,8320AG99.4% 555.9 / ‑4.6 178R146R (CGT→CGCyjjThypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (93/84);  total (93/85)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

Only 100 of 193 total aligned reads displayed.

AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604742‑4604921
                                                                                          |                                                                                         
aaaCACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcg                                                                                         >  1:3554393/1‑93 (MQ=255)
 aaCACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcga                                                                                        >  1:1305637/1‑93 (MQ=255)
 aaCACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcga                                                                                        >  1:1641635/1‑93 (MQ=255)
        cAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAt                                                                                 <  1:1515990/93‑1 (MQ=255)
        cAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAt                                                                                 <  2:4144683/93‑1 (MQ=255)
        cAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAt                                                                                 >  2:190949/1‑93 (MQ=255)
         aGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATa                                                                                <  1:1352098/93‑1 (MQ=255)
          ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               >  1:4486249/1‑93 (MQ=255)
            cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             >  1:4421018/1‑93 (MQ=255)
             aGCGTTTTGGCCGTCAGGCCATCGACGCTGTATTCACCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            <  2:350764/93‑1 (MQ=255)
             aGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            <  1:3674761/93‑1 (MQ=255)
              gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           >  2:529963/1‑93 (MQ=255)
                gTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCca                                                                         <  1:316975/93‑1 (MQ=255)
                gTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCca                                                                         >  1:3582397/1‑93 (MQ=255)
                 ttttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcac                                                                        >  2:285967/1‑93 (MQ=255)
                 ttttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaa                                                                        >  2:1231584/1‑92 (MQ=255)
                  tttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaca                                                                       >  1:3723509/1‑93 (MQ=255)
                    tGACCGTCATGCCATCGACGCTGTATGCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGc                                                                     <  2:1780799/93‑1 (MQ=255)
                    tGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGc                                                                     <  1:2652806/93‑1 (MQ=255)
                      aCCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcga                                                                   >  1:4308789/1‑93 (MQ=255)
                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  >  1:4567173/1‑93 (MQ=255)
                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  >  2:2816978/1‑93 (MQ=255)
                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  <  1:757501/93‑1 (MQ=255)
                         gTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCTAATACCGGCTGTTTTTCCATGCGACCAAAATAGAGGCCACAGcgacga                                                                >  2:4208748/1‑93 (MQ=255)
                         gTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgacga                                                                >  2:243923/1‑93 (MQ=255)
                          tCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAg                                                               <  2:122200/93‑1 (MQ=255)
                           cAGTCCATCCACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgc                                                              <  1:2643923/93‑1 (MQ=255)
                            aGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcg                                                             >  2:1623472/1‑93 (MQ=255)
                             ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            <  1:2583864/93‑1 (MQ=255)
                             ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            <  2:4405588/93‑1 (MQ=255)
                               ccATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTg                                                          >  1:4624622/1‑93 (MQ=255)
                                cATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGt                                                         >  2:3523571/1‑93 (MQ=255)
                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        <  1:31860/93‑1 (MQ=255)
                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        <  2:2772868/93‑1 (MQ=255)
                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  2:3591564/1‑93 (MQ=255)
                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  2:389395/1‑93 (MQ=255)
                                   cGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGa                                                      >  1:1061053/1‑93 (MQ=255)
                                    gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     >  2:2340569/1‑93 (MQ=255)
                                    gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     >  1:3919058/1‑93 (MQ=255)
                                    gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTACCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     <  2:3639244/93‑1 (MQ=255)
                                     aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  1:648960/93‑1 (MQ=255)
                                     aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  1:4450135/93‑1 (MQ=255)
                                       gCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTTCAGGCGACCAAAATAGAGGCCACAGCGACGACCGCGGCTGACttt                                                  >  2:2993086/1‑93 (MQ=255)
                                       gCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACttt                                                  >  1:3424396/1‑93 (MQ=255)
                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  2:895629/1‑93 (MQ=255)
                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  2:4105578/1‑93 (MQ=255)
                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  1:4473917/1‑93 (MQ=255)
                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 <  1:2607990/93‑1 (MQ=255)
                                          gTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAtt                                               >  2:3920301/1‑93 (MQ=255)
                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             <  2:4160039/93‑1 (MQ=255)
                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             >  1:2023613/1‑93 (MQ=255)
                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             >  2:4535615/1‑93 (MQ=255)
                                             ttCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCaa                                            <  1:3100909/93‑1 (MQ=255)
                                               cGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACg                                          <  1:3283058/93‑1 (MQ=255)
                                               cGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACg                                          <  2:1365069/93‑1 (MQ=255)
                                                 ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgc                                        <  2:61742/93‑1 (MQ=255)
                                                 ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgc                                        <  1:2419036/93‑1 (MQ=255)
                                                     aGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCAt                                    >  2:4634578/1‑93 (MQ=255)
                                                     aGAATTCCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCAt                                    <  1:543244/93‑1 (MQ=255)
                                                       aaTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATaa                                  <  1:3991445/93‑1 (MQ=255)
                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                >  1:841729/1‑93 (MQ=255)
                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                >  1:3133784/1‑93 (MQ=255)
                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                <  2:1059865/93‑1 (MQ=255)
                                                         ttgctcCGCATCAAATACCGGCTGTGTTGCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                <  1:2404994/90‑1 (MQ=255)
                                                             tcCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGcc                            >  1:2176641/1‑93 (MQ=255)
                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           <  2:1020807/93‑1 (MQ=255)
                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           <  2:3786660/93‑1 (MQ=255)
                                                                gCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGc                         <  1:4283035/93‑1 (MQ=255)
                                                                 cATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCa                        <  1:3404496/93‑1 (MQ=255)
                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCTACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  2:4547968/1‑93 (MQ=255)
                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  2:2792588/1‑93 (MQ=255)
                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  1:1772001/1‑93 (MQ=255)
                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  2:2264833/1‑93 (MQ=255)
                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  1:1376916/1‑93 (MQ=255)
                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  1:987126/1‑93 (MQ=255)
                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      <  1:4153494/93‑1 (MQ=255)
                                                                    cAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCt                     >  2:1279845/1‑93 (MQ=255)
                                                                     aaaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTg                    >  1:1191721/1‑93 (MQ=255)
                                                                     aaaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTg                    >  2:741804/1‑93 (MQ=255)
                                                                      aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   <  2:272272/93‑1 (MQ=255)
                                                                      aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   <  1:3783792/93‑1 (MQ=255)
                                                                      aaTACCGGCTGTTTTTCCAGGCGACAAAAATAGAGGCCACAGCGACGAACGCTGTCGACTTTATTCAACGGCGAATAATCTGCCAGCATCTGc                   >  1:41170/1‑93 (MQ=255)
                                                                       aTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCt                  >  2:4295481/1‑93 (MQ=255)
                                                                        ttctggCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  2:1431446/89‑1 (MQ=255)
                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  1:1104441/93‑1 (MQ=255)
                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 >  1:4666073/1‑93 (MQ=255)
                                                                          ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               <  2:4053336/93‑1 (MQ=255)
                                                                          ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               >  1:71946/1‑93 (MQ=255)
                                                                          ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               >  1:609301/1‑93 (MQ=255)
                                                                           cgtcTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              <  2:3705438/90‑1 (MQ=255)
                                                                           cGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              >  2:3842083/1‑93 (MQ=255)
                                                                            ggCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcg             >  1:1906526/1‑93 (MQ=255)
                                                                             gCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcgc            >  2:819846/1‑93 (MQ=255)
                                                                              cTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCt           >  1:195686/1‑93 (MQ=255)
                                                                              cTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCt           >  1:1148991/1‑93 (MQ=255)
                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  1:3336633/1‑93 (MQ=255)
                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  2:1887592/1‑93 (MQ=255)
                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  1:3872119/1‑93 (MQ=255)
                                                                                      ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   >  2:2586968/1‑93 (MQ=255)
                                                                                       cAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgcc  >  1:3474536/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604742‑4604921

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGCT  >  NC_000913/4604739‑4604924
                                                                                             |                                                                                            
GCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGG                                                                                              <  SRR4302535.826583/94‑1 (MQ=60)
  TAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCG                                                                                            >  SRR4302535.3579051/1‑94 (MQ=60)
   AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGA                                                                                           >  SRR4302535.1314957/1‑94 (MQ=60)
   AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGA                                                                                           >  SRR4302535.1653502/1‑94 (MQ=60)
        CGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAA                                                                                      <  SRR4302535.893997/94‑1 (MQ=60)
        CGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAA                                                                                      >  SRR4302535.2619787/1‑94 (MQ=60)
         GCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAA                                                                                     <  SRR4302535.1842909/94‑1 (MQ=60)
          CCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAT                                                                                    <  SRR4302535.2263630/94‑1 (MQ=60)
           CAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATA                                                                                   <  SRR4302535.3162308/94‑1 (MQ=60)
           CAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATA                                                                                   <  SRR4302535.1526909/94‑1 (MQ=60)
            AGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAG                                                                                  >  SRR4302535.4073151/1‑94 (MQ=60)
            AGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAG                                                                                  >  SRR4302535.4517044/1‑94 (MQ=60)
            AGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAG                                                                                  <  SRR4302535.1361775/94‑1 (MQ=60)
             GGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGA                                                                                 <  SRR4302535.72595/94‑1 (MQ=60)
              GCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAG                                                                                >  SRR4302535.4451292/1‑94 (MQ=60)
              GCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAG                                                                                <  SRR4302535.3352621/94‑1 (MQ=60)
                AGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGC                                                                              <  SRR4302535.3700205/94‑1 (MQ=60)
                 GCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCC                                                                             <  SRR4302535.3152095/94‑1 (MQ=60)
                  CGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCA                                                                            >  SRR4302535.3607245/1‑94 (MQ=60)
                   GTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCAC                                                                           <  SRR4302535.319176/94‑1 (MQ=60)
                   GTGTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTGTCCAGGCGACCAAAATAGAGGCCAC                                                                           >  SRR4302535.4561349/1‑94 (MQ=60)
                    TTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACA                                                                          >  SRR4302535.3749303/1‑94 (MQ=60)
                      TTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGC                                                                        <  SRR4302535.398799/94‑1 (MQ=60)
                       TGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCG                                                                       <  SRR4302535.2671810/94‑1 (MQ=60)
                       TGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCG                                                                       >  SRR4302535.3459150/1‑94 (MQ=60)
                       TGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCG                                                                       <  SRR4302535.113434/94‑1 (MQ=60)
                        GACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGA                                                                      >  SRR4302535.4338267/1‑94 (MQ=60)
                         ACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGAC                                                                     >  SRR4302535.4598326/1‑94 (MQ=60)
                          CCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACG                                                                    <  SRR4302535.762830/94‑1 (MQ=60)
                           CGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGA                                                                   >  SRR4302535.3705273/1‑94 (MQ=60)
                            GTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAG                                                                  <  SRR4302535.3663420/94‑1 (MQ=60)
                              CAGTCCATCCACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCG                                                                <  SRR4302535.2662859/94‑1 (MQ=60)
                                GGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCT                                                              <  SRR4302535.2602426/94‑1 (MQ=60)
                                GGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCT                                                              <  SRR4302535.4177506/94‑1 (MQ=60)
                                 GCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTG                                                             >  SRR4302535.4656024/1‑94 (MQ=60)
                                   CATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTC                                                           <  SRR4302535.3200596/94‑1 (MQ=60)
                                   CATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTC                                                           >  SRR4302535.3765200/1‑94 (MQ=60)
                                    ATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCG                                                          <  SRR4302535.32106/94‑1 (MQ=60)
                                     TCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGA                                                         >  SRR4302535.1068607/1‑94 (MQ=60)
                                      CGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAC                                                        >  SRR4302535.3945862/1‑94 (MQ=60)
                                        ACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTT                                                      <  SRR4302535.653483/94‑1 (MQ=60)
                                        ACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTT                                                      <  SRR4302535.4480657/94‑1 (MQ=60)
                                        ACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCGGTCGACTT                                                      >  SRR4302535.783059/1‑94 (MQ=60)
                                         CGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTT                                                     >  SRR4302535.3448210/1‑94 (MQ=60)
                                          GCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTA                                                    >  SRR4302535.4504616/1‑94 (MQ=60)
                                           CTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAT                                                   <  SRR4302535.2626704/94‑1 (MQ=60)
                                            TGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATT                                                  >  SRR4302535.3404822/1‑94 (MQ=60)
                                            TGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATT                                                  <  SRR4302535.2647931/94‑1 (MQ=60)
                                            TGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATT                                                  >  SRR4302535.4512619/1‑94 (MQ=60)
                                             GTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTC                                                 >  SRR4302535.4247487/1‑94 (MQ=60)
                                             GTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTC                                                 <  SRR4302535.1271360/94‑1 (MQ=60)
                                              TATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCA                                                >  SRR4302535.2038197/1‑94 (MQ=60)
                                               ATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGACGCCACCGCGACGAGCGCTGTCGACTTTATTCAC                                               >  SRR4302535.473918/1‑94 (MQ=60)
                                                TTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAAC                                              <  SRR4302535.3122496/94‑1 (MQ=60)
                                                TTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAAC                                              <  SRR4302535.4465378/94‑1 (MQ=60)
                                                 TCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACG                                             <  SRR4302535.1718822/94‑1 (MQ=60)
                                                  CGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGG                                            <  SRR4302535.3305861/94‑1 (MQ=60)
                                                  CGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGG                                            >  SRR4302535.327544/1‑94 (MQ=60)
                                                  CGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGG                                            <  SRR4302535.2966280/94‑1 (MQ=60)
                                                   GCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGC                                           <  SRR4302535.4338504/94‑1 (MQ=60)
                                                    CCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCG                                          <  SRR4302535.2436453/94‑1 (MQ=60)
                                                     CCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGC                                         >  SRR4302535.1178738/1‑94 (MQ=60)
                                                      CCAGAATTTCTCTGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCA                                        >  SRR4302535.2290070/1‑94 (MQ=60)
                                                        AGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATA                                      >  SRR4302535.5528/1‑94 (MQ=60)
                                                        AGAATTCCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATA                                      <  SRR4302535.547042/94‑1 (MQ=60)
                                                         GAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAA                                     >  SRR4302535.4337618/1‑94 (MQ=60)
                                                         GAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAA                                     <  SRR4302535.4660425/94‑1 (MQ=60)
                                                          AATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAT                                    <  SRR4302535.4018746/94‑1 (MQ=60)
                                                           ATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATC                                   >  SRR4302535.3155603/1‑94 (MQ=60)
                                                           ATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATC                                   >  SRR4302535.3888033/1‑94 (MQ=60)
                                                           ATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATC                                   >  SRR4302535.847634/1‑94 (MQ=60)
                                                            TTGCTCCGCATCAAATACCGGCTGTGTTGCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCT                                  <  SRR4302535.2422301/94‑1 (MQ=60)
                                                              TCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGC                                <  SRR4302535.1981035/94‑1 (MQ=60)
                                                               CTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCC                               >  SRR4302535.2192197/1‑94 (MQ=60)
                                                                  CGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACCTTATTCAACGGCGCATAATCTGCCAGC                            >  SRR4302535.1969617/1‑94 (MQ=60)
                                                                   GCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCA                           <  SRR4302535.4312337/94‑1 (MQ=60)
                                                                    CATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAGTTTATTCAACGGCGCATAATCTGCCAGCAT                          >  SRR4302535.2095588/1‑94 (MQ=60)
                                                                    CATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAT                          >  SRR4302535.1784842/1‑94 (MQ=60)
                                                                    CATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAT                          <  SRR4302535.3428187/94‑1 (MQ=60)
                                                                     ATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATC                         >  SRR4302535.994139/1‑94 (MQ=60)
                                                                     ATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATC                         <  SRR4302535.3598630/94‑1 (MQ=60)
                                                                     ATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATC                         >  SRR4302535.1386745/1‑94 (MQ=60)
                                                                     ATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATC                         >  SRR4302535.2665663/1‑94 (MQ=60)
                                                                      TCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCT                        <  SRR4302535.4181933/94‑1 (MQ=60)
                                                                       CAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTG                       >  SRR4302535.1200254/1‑94 (MQ=60)
                                                                       CAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTG                       >  SRR4302535.1564436/1‑94 (MQ=60)
                                                                        AAATACCGGCTGTTTTTCCAGGCGACAAAAATAGAGGCCACAGCGACGAACGCTGTCGACTTTATTCAACGGCGAATAATCTGCCAGCATCTGC                      >  SRR4302535.41482/1‑94 (MQ=60)
                                                                         AATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCT                     <  SRR4302535.3809899/94‑1 (MQ=60)
                                                                         AATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCT                     <  SRR4302535.1497563/94‑1 (MQ=60)
                                                                          ATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTC                    <  SRR4302535.131533/94‑1 (MQ=60)
                                                                          ATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTC                    >  SRR4302535.4697614/1‑94 (MQ=60)
                                                                          ATACCGGCTGTTTTGCCAGGCGACCAAAAGAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTC                    <  SRR4302535.136497/94‑1 (MQ=60)
                                                                           TACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCG                   <  SRR4302535.1112331/94‑1 (MQ=60)
                                                                           TACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCG                   <  SRR4302535.621314/94‑1 (MQ=60)
                                                                            ACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGG                  >  SRR4302535.72471/1‑94 (MQ=60)
                                                                            ACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGG                  >  SRR4302535.613549/1‑94 (MQ=60)
                                                                              CGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCG                >  SRR4302535.1920248/1‑94 (MQ=60)
                                                                              CGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCG                >  SRR4302535.668918/1‑94 (MQ=60)
                                                                               GGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGC               >  SRR4302535.4624520/1‑94 (MQ=60)
                                                                               GGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGC               >  SRR4302535.3927535/1‑94 (MQ=60)
                                                                                GCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCT              >  SRR4302535.1157174/1‑94 (MQ=60)
                                                                                GCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCT              >  SRR4302535.197033/1‑94 (MQ=60)
                                                                                GCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCT              <  SRR4302535.3431372/94‑1 (MQ=60)
                                                                                 CTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTG             >  SRR4302535.3898728/1‑94 (MQ=60)
                                                                                 CTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTG             >  SRR4302535.3359871/1‑94 (MQ=60)
                                                                                  TGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGC            >  SRR4302535.2533675/1‑94 (MQ=60)
                                                                                         CCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC     >  SRR4302535.3498678/1‑94 (MQ=60)
                                                                                           AGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGC   >  SRR4302535.3208122/1‑94 (MQ=60)
                                                                                           AGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTAATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGC   >  SRR4302535.2534448/1‑94 (MQ=60)
                                                                                            GGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGCT  <  SRR4302535.502322/94‑1 (MQ=60)
                                                                                            GGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGCT  <  SRR4302535.225525/94‑1 (MQ=60)
                                                                                            GGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGCT  <  SRR4302535.4179728/94‑1 (MQ=60)
                                                                                            GGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCACCATCTGCTCGGCGCTGCGCACGCCGCT  <  SRR4302535.4059966/94‑1 (MQ=60)
                                                                                            GGCGACAAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGCT  <  SRR4302535.126644/94‑1 (MQ=60)
                                                                                             |                                                                                            
GCTAAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCCGCT  >  NC_000913/4604739‑4604924

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