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 2,222,925 T→C intergenic (‑26/+139) pbpG ← / ← yohC penicillin‑binding protein 7/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,222,9250TC96.5% 413.8 / 1.1 143intergenic (‑26/+139)pbpG/yohCpenicillin‑binding protein 7/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (73/64);  total (76/67)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

AAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGT  >  NC_000913/2222835‑2223008
                                                                                          |                                                                                   
aaaGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTg                                                                                   >  1:3337952/1‑93 (MQ=255)
  aGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTgcg                                                                                 >  2:582512/1‑93 (MQ=255)
         gccagcATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCAccgcc                                                                          >  1:89094/1‑93 (MQ=255)
         gccagcATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCAccgcc                                                                          >  2:1410718/1‑93 (MQ=255)
           cggcATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCgg                                                                        >  2:4616157/3‑93 (MQ=255)
             gcATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAg                                                                      >  2:1235369/1‑93 (MQ=255)
               aTCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTa                                                                    >  2:3897413/1‑93 (MQ=255)
                tCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTaa                                                                   >  1:508281/1‑93 (MQ=255)
                tCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTaa                                                                   >  1:1376897/1‑93 (MQ=255)
                tCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTaa                                                                   >  2:2266098/1‑93 (MQ=255)
                   gggCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGTGCACCGCCGGAGTAAGGa                                                                >  1:3023280/1‑93 (MQ=255)
                    ggCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGAt                                                               <  1:4572350/93‑1 (MQ=255)
                      gtaGGCTAAATAAAGAAACTCGAAATTGCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGAttt                                                             <  1:102982/91‑1 (MQ=255)
                       cAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTa                                                            >  2:3175902/1‑93 (MQ=255)
                       cAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCAACGCCGGAGGAAGGATTTa                                                            >  2:756480/1‑93 (MQ=255)
                       cAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGATCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTa                                                            >  2:410279/1‑93 (MQ=255)
                        aGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTAc                                                           <  2:3470527/93‑1 (MQ=255)
                        aGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTAc                                                           <  1:1204769/93‑1 (MQ=255)
                         ggCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACt                                                          <  1:3982479/93‑1 (MQ=255)
                         ggCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACt                                                          >  1:3966466/1‑93 (MQ=255)
                         ggCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACt                                                          >  1:2759282/1‑93 (MQ=255)
                         ggCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACt                                                          >  1:2800550/1‑93 (MQ=255)
                            taaataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAg                                                       >  2:4160016/1‑93 (MQ=255)
                            taaataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGAGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAg                                                       >  1:663515/1‑93 (MQ=255)
                              aataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGc                                                     <  2:4342457/93‑1 (MQ=255)
                              aataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGc                                                     >  1:3408813/1‑93 (MQ=255)
                              aataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGc                                                     <  2:425086/93‑1 (MQ=255)
                              aataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGc                                                     <  2:1938564/93‑1 (MQ=255)
                              aataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGGAAGGATTTACTGAGGc                                                     >  2:470589/1‑93 (MQ=255)
                               ataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCt                                                    >  1:1781037/1‑93 (MQ=255)
                               ataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCt                                                    <  1:4254485/93‑1 (MQ=255)
                               ataaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCt                                                    >  1:4357210/1‑93 (MQ=255)
                                taaaGACCCTCGAAATCTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCGCCGCCGGAGTAAGGATTTACTGAGGCTa                                                   <  2:3334828/93‑1 (MQ=255)
                                taaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTa                                                   <  2:1451353/93‑1 (MQ=255)
                                taaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTa                                                   <  2:847721/93‑1 (MQ=255)
                                 aaaGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAg                                                  <  1:1433989/93‑1 (MQ=255)
                                 aaaGAAACACGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAg                                                  <  1:1200869/93‑1 (MQ=255)
                                    gAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGa                                               >  2:3136842/1‑93 (MQ=255)
                                    gAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGa                                               >  2:1870776/1‑93 (MQ=255)
                                    gAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGa                                               <  2:1549864/93‑1 (MQ=255)
                                    gAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGa                                               >  1:3125716/1‑93 (MQ=255)
                                     aaaCTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGAc                                              <  1:1472257/93‑1 (MQ=255)
                                      aaCTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACg                                             >  2:291498/1‑93 (MQ=255)
                                      aaCTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACg                                             >  2:658728/1‑93 (MQ=255)
                                      aaCTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACg                                             >  1:3204314/1‑93 (MQ=255)
                                         tCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCca                                          >  2:1988061/1‑93 (MQ=255)
                                         tCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCca                                          <  2:3544915/93‑1 (MQ=255)
                                          cGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcat                                         <  2:3039808/93‑1 (MQ=255)
                                           gAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcatc                                        >  2:1820106/1‑93 (MQ=255)
                                           gAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcatc                                        <  1:698642/93‑1 (MQ=255)
                                            aaaTTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcatca                                       <  2:3273817/93‑1 (MQ=255)
                                            aaaTTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcatca                                       <  2:2460326/93‑1 (MQ=255)
                                            aaaTTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcatca                                       <  2:2520715/93‑1 (MQ=255)
                                             aaTTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCcatcat                                      >  2:1118959/1‑93 (MQ=255)
                                                ttCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAAc                                   <  2:4498499/93‑1 (MQ=255)
                                                 tCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATaaag                                  <  2:3779034/93‑3 (MQ=255)
                                                  cGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGa                                 <  2:1078727/93‑1 (MQ=255)
                                                    gCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGc                               >  2:2576521/1‑93 (MQ=255)
                                                     cATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAAAGAGCa                              >  2:10039/1‑93 (MQ=255)
                                                        gatgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCaaca                           >  1:1879591/1‑91 (MQ=255)
                                                        gatgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCaaaa                           <  2:2176093/93‑1 (MQ=255)
                                                        gatgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCaaaa                           >  2:318492/1‑93 (MQ=255)
                                                        gatgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCaaaa                           >  1:46625/1‑93 (MQ=255)
                                                         atgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCaaaaa                          >  1:261417/1‑93 (MQ=255)
                                                         atgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCaaaaa                          >  2:2891106/1‑93 (MQ=255)
                                                          tgaGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAg                         >  1:4298353/1‑93 (MQ=255)
                                                             gCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGc                      >  1:3447124/1‑93 (MQ=255)
                                                             gCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGc                      >  2:834729/1‑93 (MQ=255)
                                                               aTTCAGATAGTGGTGCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGa                    <  1:1124599/93‑1 (MQ=255)
                                                                ttCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGAt                   <  2:1563384/93‑1 (MQ=255)
                                                                ttCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGAt                   <  1:4155717/93‑1 (MQ=255)
                                                                ttCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGAt                   <  2:4004064/93‑1 (MQ=255)
                                                                ttCAGATAGTGGTTCACGCGCACGGGATGCGCACCGCCGGAGTAAGGATTTAATGAGGCTAGCGACGCCATCATAAAGCACAAACAGTGCGAt                   >  2:2843889/1‑93 (MQ=255)
                                                                 tCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATc                  <  1:1674744/93‑1 (MQ=255)
                                                                   aGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGt                >  1:2056962/1‑93 (MQ=255)
                                                                    gATAGTGGTTCACGGGGCCGCGCTGCGGACCAGCGGGGAGAGGACTTACTGAGGCTAGTGACGACAACATAACGAGCAAAAAGCGCGATCGTc               >  2:1107846/1‑93 (MQ=255)
                                                                     aTAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCa              <  1:4544536/93‑1 (MQ=255)
                                                                     aTAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCa              >  1:362441/1‑93 (MQ=255)
                                                                      tAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCaa             >  1:1348725/1‑93 (MQ=255)
                                                                      tAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCaa             >  1:1727005/1‑93 (MQ=255)
                                                                      tAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCaa             >  2:2656169/1‑93 (MQ=255)
                                                                       aGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCaaa            >  1:879828/1‑93 (MQ=255)
                                                                       aGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCaaa            >  1:1523053/1‑93 (MQ=255)
                                                                        gTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAg           <  2:4301735/93‑1 (MQ=255)
                                                                          ggTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGga         <  2:4370995/93‑1 (MQ=255)
                                                                          ggTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGga         <  1:2203383/93‑1 (MQ=255)
                                                                           gTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGgag        >  2:1147326/1‑93 (MQ=255)
                                                                             tCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGaa      <  1:1646605/93‑1 (MQ=255)
                                                                             tCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGaa      <  1:2566404/93‑1 (MQ=255)
                                                                             tCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGaa      <  2:3502075/93‑1 (MQ=255)
                                                                              cACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAAt     <  2:2833784/93‑1 (MQ=255)
                                                                              cACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAAt     <  1:1986045/93‑1 (MQ=255)
                                                                               aCGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATc    <  2:3236374/93‑1 (MQ=255)
                                                                               aCGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATc    <  2:627532/93‑1 (MQ=255)
                                                                                cgcgCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCg   >  2:3897140/1‑93 (MQ=255)
                                                                                cgcgCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCg   <  2:4505072/93‑1 (MQ=255)
                                                                                cgcgCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCg   >  1:3614329/1‑93 (MQ=255)
                                                                                 gcgcACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGt  <  1:1564500/93‑1 (MQ=255)
                                                                                 gcgcACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGt  <  2:3343041/93‑1 (MQ=255)
                                                                                 gcgcACGGGCTGCGCACCGCAGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGt  <  2:3553856/93‑1 (MQ=255)
                                                                                          |                                                                                   
AAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGT  >  NC_000913/2222835‑2223008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ 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
AAAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATG  >  NC_000913/2222834‑2223017
                                                                                           |                                                                                            
AAAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTG                                                                                            >  SRR4302535.3361202/1‑94 (MQ=60)
         AGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCC                                                                                   >  SRR4302535.3483961/1‑94 (MQ=60)
         AGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCC                                                                                   >  SRR4302535.89724/1‑94 (MQ=60)
          GCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCG                                                                                  <  SRR4302535.4031764/94‑1 (MQ=60)
                ATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAA                                                                            >  SRR4302535.1386726/1‑94 (MQ=60)
                ATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAA                                                                            >  SRR4302535.511837/1‑94 (MQ=60)
                 TCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAG                                                                           <  SRR4302535.393753/94‑1 (MQ=60)
                   AGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGTGCACCGCCGGAGTAAGGA                                                                         >  SRR4302535.3044447/1‑94 (MQ=60)
                     GGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATT                                                                       <  SRR4302535.4603521/94‑1 (MQ=60)
                       gtAGGCTAAATAAAGAAACTCGAAATTGCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTA                                                                     <  SRR4302535.103712/92‑1 (MQ=60)
                       CCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTA                                                                     >  SRR4302535.4624008/1‑94 (MQ=60)
                        CAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTAC                                                                    <  SRR4302535.4445029/94‑1 (MQ=60)
                         AGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACT                                                                   <  SRR4302535.1213420/94‑1 (MQ=60)
                         AGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACT                                                                   >  SRR4302535.2820472/1‑94 (MQ=60)
                         AGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACT                                                                   >  SRR4302535.2778944/1‑94 (MQ=60)
                         AGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACT                                                                   >  SRR4302535.3993559/1‑94 (MQ=60)
                          GGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTG                                                                  <  SRR4302535.1176367/94‑1 (MQ=60)
                          GGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTG                                                                  <  SRR4302535.4009689/94‑1 (MQ=60)
                           GCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGA                                                                 <  SRR4302535.924476/94‑1 (MQ=60)
                            CTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAG                                                                >  SRR4302535.2710098/1‑94 (MQ=60)
                            CTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGAGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAG                                                                >  SRR4302535.668151/1‑94 (MQ=60)
                              AAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGC                                                              >  SRR4302535.3432538/1‑94 (MQ=60)
                               AATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCT                                                             >  SRR4302535.1793925/1‑94 (MQ=60)
                               AATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCT                                                             >  SRR4302535.4387039/1‑94 (MQ=60)
                                ATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTA                                                            <  SRR4302535.4283632/94‑1 (MQ=60)
                                  AAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGG                                                          >  SRR4302535.4309513/1‑94 (MQ=60)
                                  AAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGC                                                          <  SRR4302535.1444227/94‑1 (MQ=60)
                                  AAAGAAACACGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGC                                                          <  SRR4302535.1209472/94‑1 (MQ=60)
                                    AGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGA                                                        >  SRR4302535.3147487/1‑94 (MQ=60)
                                      AAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACG                                                      <  SRR4302535.1482798/94‑1 (MQ=60)
                                      AAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACG                                                      >  SRR4302535.3226505/1‑94 (MQ=60)
                                        ACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCC                                                    >  SRR4302535.3405475/1‑94 (MQ=60)
                                            GAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCA                                                <  SRR4302535.703543/94‑1 (MQ=60)
                                               ATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAA                                             <  SRR4302535.1843360/94‑1 (MQ=60)
                                                   CGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAG                                         >  SRR4302535.283953/1‑94 (MQ=60)
                                                     GCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCA                                       <  SRR4302535.129193/94‑1 (MQ=60)
                                                        TGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAACA                                    >  SRR4302535.1893134/1‑94 (MQ=60)
                                                        TGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAA                                    >  SRR4302535.46977/1‑94 (MQ=60)
                                                        TGATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGGAAGGATTTACTGAGGCTAGCGACGCACTCATAACGAGCAAAA                                    >  SRR4302535.1504100/1‑94 (MQ=60)
                                                         GATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAA                                   >  SRR4302535.263229/1‑94 (MQ=60)
                                                          ATGAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAG                                  >  SRR4302535.4327767/1‑94 (MQ=60)
                                                            GAGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTG                                <  SRR4302535.1396023/94‑1 (MQ=60)
                                                             AGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGC                               <  SRR4302535.3154459/94‑1 (MQ=60)
                                                             AGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGC                               >  SRR4302535.3471093/1‑94 (MQ=60)
                                                             AGCATTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGC                               <  SRR4302535.1726625/94‑1 (MQ=60)
                                                               CATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGA                             >  SRR4302535.2352525/1‑94 (MQ=60)
                                                                ATTCAGATAGTGGTGCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGAT                            <  SRR4302535.1132623/94‑1 (MQ=60)
                                                                 TTCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATC                           <  SRR4302535.4184166/94‑1 (MQ=60)
                                                                  TCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCG                          <  SRR4302535.1686870/94‑1 (MQ=60)
                                                                  TCAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCG                          <  SRR4302535.3705463/94‑1 (MQ=60)
                                                                   CAGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGT                         >  SRR4302535.2071781/1‑94 (MQ=60)
                                                                    AGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTC                        >  SRR4302535.2064998/1‑94 (MQ=60)
                                                                    AGATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTC                        <  SRR4302535.1645184/94‑1 (MQ=60)
                                                                     GATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCA                       >  SRR4302535.364941/1‑94 (MQ=60)
                                                                      ATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAA                      <  SRR4302535.4575617/94‑1 (MQ=60)
                                                                      ATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAA                      >  SRR4302535.1358358/1‑94 (MQ=60)
                                                                      ATAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAA                      >  SRR4302535.1739490/1‑94 (MQ=60)
                                                                       TAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAA                     >  SRR4302535.1534024/1‑94 (MQ=60)
                                                                       TAGTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAA                     >  SRR4302535.886057/1‑94 (MQ=60)
                                                                         GTGGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGG                   <  SRR4302535.555670/94‑1 (MQ=60)
                                                                           GGTTCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAG                 <  SRR4302535.2219083/94‑1 (MQ=60)
                                                                              TCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAAT              <  SRR4302535.1658505/94‑1 (MQ=60)
                                                                              TCACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAAT              <  SRR4302535.2584867/94‑1 (MQ=60)
                                                                               CACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATC             <  SRR4302535.2000348/94‑1 (MQ=60)
                                                                                ACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCG            >  SRR4302535.3639385/1‑94 (MQ=60)
                                                                                ACGCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCG            <  SRR4302535.838812/94‑1 (MQ=60)
                                                                                  GCGCACGGGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTG          <  SRR4302535.1575796/94‑1 (MQ=60)
                                                                                         GGCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAAT   <  SRR4302535.277251/94‑1 (MQ=60)
                                                                                          GCTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATG  <  SRR4302535.1666086/94‑1 (MQ=60)
                                                                                           |                                                                                            
AAAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATG  >  NC_000913/2222834‑2223017

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