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 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,494,539 T→C F222S (TTC→TCC)  nanX → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,494,5390TC97.3% 235.2 / ‑3.1 75F222S (TTC→TCC) nanXMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (33/40);  total (34/41)
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.98e-01

AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG  >  NZ_CP009273/4494452‑4494625
                                                                                       |                                                                                      
aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                      >  2:704810/1‑90 (MQ=255)
aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                      <  2:619588/90‑1 (MQ=255)
aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                      <  1:35086/90‑1 (MQ=255)
    gTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGt                                                                                  >  1:323762/1‑90 (MQ=255)
          tGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttt                                                                            <  2:1552748/90‑1 (MQ=255)
           gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGtttt                                                                           >  1:1469848/1‑90 (MQ=255)
               aTAAATGTAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTc                                                                       <  1:570121/90‑1 (MQ=255)
                tAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCt                                                                      >  1:29467/1‑90 (MQ=255)
                  aaTATAAGGCTAAATGAACATTTTTGTCTGTCTTCAGAAAACCAAATCTTGCAATCTCTATGATCGTTTCCCTGTTCTGTTTTTGTCTAt                                                                    >  1:732484/1‑90 (MQ=255)
                   atatAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTAtt                                                                   <  2:559546/90‑1 (MQ=255)
                    tataAGGATAAATCAACATTTTTGTCTGTCTTCAGAGAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTAttt                                                                  <  1:1310616/90‑1 (MQ=255)
                        aGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTg                                                              >  2:416695/1‑90 (MQ=255)
                          gATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCa                                                            >  2:17353/1‑90 (MQ=255)
                           aTAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaa                                                           >  1:1393251/1‑90 (MQ=255)
                            tAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaaa                                                          >  2:650620/1‑90 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttttt                                                                          >  1:1452315/1‑73 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttttt                                                                          <  2:1452315/73‑1 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         <  1:1039169/90‑1 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         >  2:1452075/1‑90 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         >  2:1005286/1‑90 (MQ=255)
                                   aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                   >  2:225165/1‑90 (MQ=255)
                                   aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                   >  1:871670/1‑90 (MQ=255)
                                   aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                   >  1:647973/1‑90 (MQ=255)
                                        tttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATaaa                                              >  1:1535015/1‑90 (MQ=255)
                                        tttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATaaa                                              >  2:513027/1‑90 (MQ=255)
                                         ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                             <  2:1560440/90‑1 (MQ=255)
                                             ctgtctTCAGAAACCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  1:1291622/90‑1 (MQ=255)
                                             ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  1:458730/90‑1 (MQ=255)
                                             ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  2:1013665/90‑1 (MQ=255)
                                              tgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGact                                        >  2:271609/1‑90 (MQ=255)
                                                 ctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGactact                                     <  1:158230/90‑1 (MQ=255)
                                                  ttCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  1:768161/90‑1 (MQ=255)
                                                  ttCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  2:1022622/90‑1 (MQ=255)
                                                  ttCAGAAAACCACATCTTTCAATCTCTAAGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  1:704810/90‑1 (MQ=255)
                                                      gAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActtcctt                                <  2:1409075/90‑1 (MQ=255)
                                                      gAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActtcctt                                <  2:389659/90‑1 (MQ=255)
                                                             acaTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCt                         >  1:1094028/1‑90 (MQ=255)
                                                               aTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTgg                       <  2:871670/90‑1 (MQ=255)
                                                                 cTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCa                     <  1:246617/90‑1 (MQ=255)
                                                                  tttCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAg                    >  2:128725/1‑90 (MQ=255)
                                                                     cAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATa                 >  2:119184/1‑90 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGGCTACTTCCTTCCTACCTGGCAGATAATg              <  2:1099230/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  2:1094028/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  2:922563/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:984864/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:128725/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1439377/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1452075/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1551310/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:162034/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:78254/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:873382/90‑1 (MQ=255)
                                                                          tctATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:678543/90‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         >  2:1258513/1‑61 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  1:1258513/61‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  1:1145631/90‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:851440/90‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:1250067/90‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:296232/90‑1 (MQ=255)
                                                                           ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:1375991/1‑90 (MQ=255)
                                                                           ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  2:1408662/1‑90 (MQ=255)
                                                                           ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:596260/1‑90 (MQ=255)
                                                                           ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  2:1065878/1‑90 (MQ=255)
                                                                           ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:833982/1‑90 (MQ=255)
                                                                           ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCCGGCAGATAATGGAg           >  1:1030930/1‑90 (MQ=255)
                                                                             aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt         <  1:1161139/90‑1 (MQ=255)
                                                                             aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt         >  2:618734/1‑90 (MQ=255)
                                                                                  cGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCCGATAATGGAGTTAATAc    <  1:1576399/90‑1 (MQ=255)
                                                                                  cGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAc    <  2:1178305/90‑1 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTa                            <  2:28421/65‑1 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTa                            >  1:28421/1‑65 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              >  2:468043/1‑79 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1290541/79‑1 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:468043/79‑1 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              >  2:1290541/1‑79 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  2:429824/1‑90 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  1:1485912/1‑90 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  2:217319/1‑90 (MQ=255)
                                                                                    tttcCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAg  <  1:8079/90‑1 (MQ=255)
                                                                                       |                                                                                      
AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG  >  NZ_CP009273/4494452‑4494625

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  NZ_CP009273/4494446‑4494635
                                                                                             |                                                                                                
AAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGT                                                                                            >  SRR3722078.327727/1‑100 (MQ=60)
      AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTT                                                                                      <  SRR3722078.35515/100‑1 (MQ=60)
       GGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTT                                                                                     >  SRR3722078.1493817/1‑100 (MQ=60)
            GGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCT                                                                                >  SRR3722078.29832/1‑100 (MQ=60)
              GTTGAAGATAAATATAAGGCTAAATGAACATTTTTGTCTGTCTTCAGAAAACCAAATCTTGCAATCTCTATGATCGTTTCCCTGTTCTGTTTTTGTCTAT                                                                              >  SRR3722078.741656/1‑100 (MQ=60)
                     gatctacacgcccgccatcgtcggcagcgtcagatgtgtataagagacagCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGC                                                                       <  SRR3722078.177372/50‑1 (MQ=60)
                     ATAAATGTAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGC                                                                       <  SRR3722078.576799/100‑1 (MQ=60)
                       AAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAA                                                                     >  SRR3722078.1415764/1‑100 (MQ=60)
                         ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGctgtctc                                                                   >  SRR3722078.1475995/1‑93 (MQ=60)
                          TATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAGAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACT                                                                  <  SRR3722078.1331540/100‑1 (MQ=60)
                               GGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCG                                                             >  SRR3722078.883349/1‑100 (MQ=60)
                               GGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCG                                                             >  SRR3722078.655418/1‑100 (MQ=60)
                                   AAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAA                                                         <  SRR3722078.1054663/100‑1 (MQ=60)
                                    AATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAA                                                        >  SRR3722078.1560036/1‑100 (MQ=60)
                                                   agatgtgtataagagacagTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTC                                         <  SRR3722078.1278463/81‑1 (MQ=60)
                                                   CTGTCTTCAGAAACCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTC                                         <  SRR3722078.1312183/100‑1 (MQ=60)
                                                   CTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTC                                         <  SRR3722078.464208/100‑1 (MQ=60)
                                                       CTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTAC                                     <  SRR3722078.160077/100‑1 (MQ=60)
                                                        TTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACC                                    <  SRR3722078.778069/100‑1 (MQ=60)
                                                        TTCAGAAAACCACATCTTTCAATCTCTAAGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACC                                    <  SRR3722078.713333/100‑1 (MQ=60)
                                                         TCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCT                                   >  SRR3722078.1110654/1‑100 (MQ=60)
                                                                       CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                     >  SRR3722078.845072/1‑100 (MQ=60)
                                                                       CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                     <  SRR3722078.249628/100‑1 (MQ=60)
                                                                       CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                     >  SRR3722078.1398174/1‑100 (MQ=60)
                                                                       CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                     >  SRR3722078.603176/1‑100 (MQ=60)
                                                                       CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCCGGCAGATAATGGAG                     >  SRR3722078.1046235/1‑100 (MQ=60)
                                                                              gccCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.473616/97‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.885096/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.79179/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.999108/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.163924/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.130212/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.1462804/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.1576579/100‑1 (MQ=60)
                                                                              TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.1475751/100‑1 (MQ=60)
                                                                              GCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722078.1311082/100‑1 (MQ=60)
                                                                               CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACctgtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatctcgg             >  SRR3722078.1606469/1‑46 (MQ=60)
                                                                               CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTctgtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatctcg             >  SRR3722078.99708/1‑47 (MQ=60)
                                                                               CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATctgtctcttatacac             >  SRR3722078.28770/1‑85 (MQ=60)
                                                                               CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA             >  SRR3722078.1510138/1‑100 (MQ=60)
                                                                                TCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG            <  SRR3722078.1163273/100‑1 (MQ=60)
                                                                                 CTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGT           >  SRR3722078.1268593/1‑100 (MQ=60)
                                                                                 CTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGT           >  SRR3722078.1093275/1‑100 (MQ=60)
                                                                                   ATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGG         <  SRR3722078.1179132/100‑1 (MQ=60)
                                                                                       TCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCAT     >  SRR3722078.1115079/1‑100 (MQ=60)
                                                                                        CGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCCGATAATGGAGTTAATACAGTGGTCATT    <  SRR3722078.1602055/100‑1 (MQ=60)
                                                                                         GTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACctgtctcttat   >  SRR3722078.663632/1‑89 (MQ=60)
                                                                                          TTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  <  SRR3722078.8163/100‑1 (MQ=60)
                                                                                             |                                                                                                
AAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  NZ_CP009273/4494446‑4494635

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