Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,847,298 A→G K146K (AAA→AAG sppA → protease IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,847,2980AG99.3% 438.7 / NA 136K146K (AAA→AAGsppAprotease IV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (58/77);  total (59/77)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387
                                                                                          |                                                                                         
 gtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTGTGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCg                                                                                        >  2:2257409/1‑93 (MQ=255)
 gtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCg                                                                                        <  2:793682/93‑1 (MQ=255)
 gtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCg                                                                                        >  2:2257410/1‑93 (MQ=255)
   gATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGt                                                                                      <  1:1361913/93‑1 (MQ=255)
   gATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGt                                                                                      >  2:3062596/1‑93 (MQ=255)
   gATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGt                                                                                      <  1:2214334/93‑1 (MQ=255)
    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  1:687870/93‑1 (MQ=255)
     tGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGttt                                                                                    <  1:1464356/93‑1 (MQ=255)
      ggATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  2:370090/93‑1 (MQ=255)
      gcagcTGAAAAGCTTCGCAGGCGGCGTCCAACCGTCGATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  2:370089/89‑1 (MQ=255)
        aTCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATg                                                                                 <  2:1489838/93‑1 (MQ=255)
        aTCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATg                                                                                 >  1:138700/1‑93 (MQ=255)
          cTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGcc                                                                               <  2:2747078/93‑1 (MQ=255)
           tGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCg                                                                              >  2:975452/1‑93 (MQ=255)
            gAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGt                                                                             >  1:2563922/1‑93 (MQ=255)
              aaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTg                                                                           <  1:3572373/93‑1 (MQ=255)
              aaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTg                                                                           <  2:3865338/93‑1 (MQ=255)
               aaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTgg                                                                          >  1:1919008/1‑93 (MQ=255)
               aaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGGTgg                                                                          >  1:4298950/1‑93 (MQ=255)
                aaCTTCGCAGGCGGAGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGAGGGAAGCCGGTTTATGCCGTTGGc                                                                         <  2:2533799/93‑1 (MQ=255)
                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCga                                                                       <  2:4631955/93‑1 (MQ=255)
                   ttCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCgag                                                                      >  2:1615261/1‑93 (MQ=255)
                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   <  1:3873562/93‑1 (MQ=255)
                      gCAGGCGGCGACCAACCATCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   <  1:4146637/93‑1 (MQ=255)
                       cAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACt                                                                  >  1:2782046/1‑93 (MQ=255)
                         ggcggcggcCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                <  1:3300036/93‑1 (MQ=255)
                         ggcggcGACCACCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                <  2:4415734/93‑1 (MQ=255)
                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                <  1:2063228/93‑1 (MQ=255)
                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               >  2:3591933/1‑93 (MQ=255)
                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  2:1169265/93‑1 (MQ=255)
                           cggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAg                                                              <  1:609451/93‑1 (MQ=255)
                            ggcgACCAACCGTCTATGCAGTACATCGGCAAAGGTCTGAAAAAGTTTAGTGACAGAGGGAAACCGGTTTATGCCATTGGCTAGAACTACAGc                                                             >  2:2523798/1‑93 (MQ=255)
                            ggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGc                                                             <  1:4559781/93‑1 (MQ=255)
                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGcc                                                            <  1:3213142/93‑1 (MQ=255)
                                aCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                         <  2:1125356/93‑1 (MQ=255)
                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        <  1:2131383/93‑1 (MQ=255)
                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        >  2:155246/1‑93 (MQ=255)
                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        <  1:1882838/93‑1 (MQ=255)
                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        <  1:2055281/93‑1 (MQ=255)
                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCAGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        >  1:1658568/1‑93 (MQ=255)
                                   aaCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                      >  2:822651/1‑93 (MQ=255)
                                   aaCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                      >  2:3877901/1‑93 (MQ=255)
                                   aaCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                      <  2:98048/93‑1 (MQ=255)
                                     ccGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCaa                                                    <  2:2468187/93‑1 (MQ=255)
                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  <  1:1084096/93‑1 (MQ=255)
                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  <  1:1398276/93‑1 (MQ=255)
                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 <  1:1764918/93‑1 (MQ=255)
                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  2:736056/1‑93 (MQ=255)
                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 <  2:1376029/93‑1 (MQ=255)
                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAAAGCCAGGGGCAatat                                                 >  2:283512/1‑93 (MQ=255)
                                          tATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatta                                               <  2:546400/93‑1 (MQ=255)
                                          tATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatta                                               <  1:3959066/93‑1 (MQ=255)
                                            tGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtc                                             <  1:3562913/93‑1 (MQ=255)
                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            <  1:85251/93‑1 (MQ=255)
                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            <  2:407950/93‑1 (MQ=255)
                                              cAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtctc                                           <  1:3220717/93‑1 (MQ=255)
                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         >  1:1348719/1‑93 (MQ=255)
                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         <  1:3628463/93‑1 (MQ=255)
                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         >  1:1106284/1‑93 (MQ=255)
                                                 tACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGcc                                        <  2:1087020/93‑1 (MQ=255)
                                                 tACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGcc                                        <  1:621134/93‑1 (MQ=255)
                                                 tACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGcc                                        <  1:3365611/93‑1 (MQ=255)
                                                  aCATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCa                                       <  2:1345289/93‑1 (MQ=255)
                                                  aCATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCa                                       <  2:4345901/93‑1 (MQ=255)
                                                   cATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAg                                      <  2:4552465/93‑1 (MQ=255)
                                                    aTCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGt                                     >  1:1620471/1‑93 (MQ=255)
                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGTGGCAATATTATCTCGCCAGTTTc                                  >  2:3393872/1‑93 (MQ=255)
                                                        gCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCg                                 >  2:2127155/1‑93 (MQ=255)
                                                           aaGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCa                              <  2:3218616/93‑1 (MQ=255)
                                                           aaGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCa                              >  1:2359857/1‑93 (MQ=255)
                                                            aGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCaa                             >  1:2093232/1‑93 (MQ=255)
                                                            aGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCaa                             >  1:574285/1‑93 (MQ=255)
                                                             gCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAAt                            >  2:1153612/1‑93 (MQ=255)
                                                              ctctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATa                           >  1:3408005/1‑93 (MQ=255)
                                                              ctctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATa                           <  2:2059577/93‑1 (MQ=255)
                                                                 tGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGACGTTGGCGAGAACTACAGCCACGGGCAATATTATCTCGTCAGTTTCGCCAATaaaa                        >  2:255332/1‑93 (MQ=255)
                                                                   aaaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGGAATATTATATCGCAAGTTTGGCCAATAAAAtt                      >  1:1900535/1‑93 (MQ=255)
                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     >  1:2559375/1‑93 (MQ=255)
                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     <  1:3673710/93‑1 (MQ=255)
                                                                      gagTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTgg                   <  2:4527972/93‑1 (MQ=255)
                                                                        gTTTGGTGACAGCGGGAATGCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGCCAATATTATCTCGCCAGACTCGCCAATAAAATTTGGCt                 >  2:872450/1‑93 (MQ=255)
                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 <  2:924448/93‑1 (MQ=255)
                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 >  2:872449/1‑93 (MQ=255)
                                                                         tttCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTg                <  1:1158618/93‑1 (MQ=255)
                                                                         tttCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTg                >  2:4561195/1‑93 (MQ=255)
                                                                           tCGTGGCAACGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtc              <  2:2841101/93‑1 (MQ=255)
                                                                           tCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtc              <  1:4679467/93‑1 (MQ=255)
                                                                            cGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtct             >  1:366981/1‑93 (MQ=255)
                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            <  2:1498222/93‑1 (MQ=255)
                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            <  1:4283653/93‑1 (MQ=255)
                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            >  2:2154747/1‑93 (MQ=255)
                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            >  1:3011788/1‑93 (MQ=255)
                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           <  2:1599049/93‑1 (MQ=255)
                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           <  2:4660233/93‑1 (MQ=255)
                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           >  2:1895453/1‑93 (MQ=255)
                                                                                  aGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCaa       >  1:1217244/1‑93 (MQ=255)
                                                                                  aGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCaa       <  1:4607139/93‑1 (MQ=255)
                                                                                     gggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGc    <  2:1143909/93‑1 (MQ=255)
                                                                                     gggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGc    <  1:1439836/93‑1 (MQ=255)
                                                                                      ggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCg   >  2:61658/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGT  >  NC_000913/1847205‑1847390
                                                                                             |                                                                                            
TATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAG                                                                                              >  SRR4302544.1729132/1‑94 (MQ=60)
 ATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGC                                                                                             >  SRR4302544.123482/1‑94 (MQ=60)
  TTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCC                                                                                            <  SRR4302544.1012688/94‑1 (MQ=60)
      GATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTT                                                                                        <  SRR4302544.2235166/94‑1 (MQ=60)
      GATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTT                                                                                        <  SRR4302544.1375436/94‑1 (MQ=60)
       ATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTT                                                                                       <  SRR4302544.694715/94‑1 (MQ=60)
        TGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTA                                                                                      <  SRR4302544.1478869/94‑1 (MQ=60)
          GATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATG                                                                                    >  SRR4302544.140013/1‑94 (MQ=60)
              TGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGT                                                                                >  SRR4302544.2588555/1‑94 (MQ=60)
                 AAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGG                                                                             <  SRR4302544.672148/94‑1 (MQ=60)
                 AAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGG                                                                             >  SRR4302544.1937549/1‑94 (MQ=60)
                 AAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGG                                                                             <  SRR4302544.3606629/94‑1 (MQ=60)
                 AAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGGTGG                                                                             >  SRR4302544.4342820/1‑94 (MQ=60)
                    ACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGA                                                                          >  SRR4302544.4762372/1‑94 (MQ=60)
                     CTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAG                                                                         <  SRR4302544.2362815/94‑1 (MQ=60)
                        CGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAC                                                                      >  SRR4302544.3815439/1‑94 (MQ=60)
                         GCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACT                                                                     <  SRR4302544.3912031/94‑1 (MQ=60)
                         GCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACT                                                                     >  SRR4302544.2808822/1‑94 (MQ=60)
                         GCAGGCGGCGACCAACCATCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACT                                                                     <  SRR4302544.4187691/94‑1 (MQ=60)
                            GGCGGCGGCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACA                                                                  <  SRR4302544.3331468/94‑1 (MQ=60)
                            GGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACA                                                                  <  SRR4302544.2083064/94‑1 (MQ=60)
                            GGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACA                                                                  <  SRR4302544.1211063/94‑1 (MQ=60)
                              CGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGC                                                                <  SRR4302544.615437/94‑1 (MQ=60)
                               GGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCC                                                               <  SRR4302544.4605260/94‑1 (MQ=60)
                                GCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCA                                                              <  SRR4302544.3243786/94‑1 (MQ=60)
                                  GACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGG                                                            <  SRR4302544.1809710/94‑1 (MQ=60)
                                   ACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCAGGTTTATGCCGTTGGCGAGAACTACAGCCAGGG                                                           >  SRR4302544.1674905/1‑94 (MQ=60)
                                    CCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGG                                                          <  SRR4302544.1901038/94‑1 (MQ=60)
                                    CCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGG                                                          <  SRR4302544.2075058/94‑1 (MQ=60)
                                    CCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGG                                                          <  SRR4302544.2151789/94‑1 (MQ=60)
                                          GTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATAT                                                    <  SRR4302544.1094880/94‑1 (MQ=60)
                                          GTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATAT                                                    <  SRR4302544.1412214/94‑1 (MQ=60)
                                           TCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATT                                                   <  SRR4302544.1782187/94‑1 (MQ=60)
                                             TATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAT                                                 <  SRR4302544.3998436/94‑1 (MQ=60)
                                               TGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCT                                               >  SRR4302544.1665172/1‑94 (MQ=60)
                                               TGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCT                                               <  SRR4302544.3597072/94‑1 (MQ=60)
                                                GCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTC                                              <  SRR4302544.86058/94‑1 (MQ=60)
                                                 CAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCG                                             <  SRR4302544.3251452/94‑1 (MQ=60)
                                                  AGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGTCGAGAACTACAGCCAGGGGCAATATTATCTCGC                                            >  SRR4302544.4611625/1‑94 (MQ=60)
                                                  AGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGC                                            >  SRR4302544.1117309/1‑94 (MQ=60)
                                                  AGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGC                                            >  SRR4302544.1362087/1‑94 (MQ=60)
                                                   GTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCC                                           <  SRR4302544.3664613/94‑1 (MQ=60)
                                                    TACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCA                                          <  SRR4302544.3397735/94‑1 (MQ=60)
                                                    TACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCA                                          <  SRR4302544.627228/94‑1 (MQ=60)
                                                      CATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGT                                        >  SRR4302544.1636432/1‑94 (MQ=60)
                                                         CGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTC                                     >  SRR4302544.3153067/1‑94 (MQ=60)
                                                         CGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTC                                     >  SRR4302544.2684136/1‑94 (MQ=60)
                                                            CAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGCTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCC                                  >  SRR4302544.138746/1‑94 (MQ=60)
                                                             AAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCA                                 >  SRR4302544.2382199/1‑94 (MQ=60)
                                                              AAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAA                                >  SRR4302544.2113382/1‑94 (MQ=60)
                                                              AAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAA                                >  SRR4302544.579899/1‑94 (MQ=60)
                                                              AAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAA                                >  SRR4302544.2113380/1‑94 (MQ=60)
                                                                GCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATA                              >  SRR4302544.3440592/1‑94 (MQ=60)
                                                                   CTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAA                           >  SRR4302544.4269171/1‑94 (MQ=60)
                                                                    TGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAT                          >  SRR4302544.3377162/1‑94 (MQ=60)
                                                                     GAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGGAATATTATATCGCAAGTTTGGCCAATAAAATT                         >  SRR4302544.1918897/1‑94 (MQ=60)
                                                                     GAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATT                         >  SRR4302544.2042743/1‑94 (MQ=60)
                                                                      AAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTT                        >  SRR4302544.2583955/1‑94 (MQ=60)
                                                                       AAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTG                       <  SRR4302544.3710350/94‑1 (MQ=60)
                                                                          AGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCT                    >  SRR4302544.3485567/1‑94 (MQ=60)
                                                                            TTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGT                  <  SRR4302544.1170151/94‑1 (MQ=60)
                                                                              TCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCT                >  SRR4302544.370734/1‑94 (MQ=60)
                                                                              TCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCT                <  SRR4302544.4725729/94‑1 (MQ=60)
                                                                              TCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCT                <  SRR4302544.885149/94‑1 (MQ=60)
                                                                               CGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTC               >  SRR4302544.3040653/1‑94 (MQ=60)
                                                                                GTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCC              <  SRR4302544.4327409/94‑1 (MQ=60)
                                                                                    CAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAA          >  SRR4302544.1229288/1‑94 (MQ=60)
                                                                                     AGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAG         <  SRR4302544.4652907/94‑1 (MQ=60)
                                                                                        GGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCG      <  SRR4302544.141637/94‑1 (MQ=60)
                                                                                        GGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCG      <  SRR4302544.1454165/94‑1 (MQ=60)
                                                                                          GAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTG    <  SRR4302544.3702828/94‑1 (MQ=60)
                                                                                           AAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGG   >  SRR4302544.1846487/1‑94 (MQ=60)
                                                                                           AAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGG   >  SRR4302544.869123/1‑94 (MQ=60)
                                                                                            AGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCTAATAAAATTTGGCTGTCTCCGCAAGGCGTGGT  >  SRR4302544.3449023/1‑94 (MQ=60)
                                                                                            AGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGT  >  SRR4302544.4629231/1‑94 (MQ=60)
                                                                                             |                                                                                            
TATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGT  >  NC_000913/1847205‑1847390

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