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 3,950,557 2 bp→GC coding (451‑452/930 nt) ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,5570CG100.0% 684.7 / NA 213R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (99/114);  total (99/114)
*NC_0009133,950,5580GC99.5% 664.3 / ‑5.9 213R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (98/114);  total (99/114)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 220 total aligned reads displayed.

GATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT  >  NC_000913/3950466‑3950643
                                                                                           ||                                                                                     
    aTCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCa                                                                                   >  1:3630922/1‑93 (MQ=255)
    aTCGCTGCTGTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCa                                                                                   <  2:3161389/93‑1 (MQ=255)
       gctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                >  1:3793108/1‑93 (MQ=255)
        ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                                               >  2:766991/1‑93 (MQ=255)
         tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                              <  1:453290/93‑1 (MQ=255)
         tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                              >  1:4222280/1‑93 (MQ=255)
          gctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                                             <  1:3769707/93‑1 (MQ=255)
          gctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                                             >  2:634028/1‑93 (MQ=255)
           ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                            <  2:4675144/93‑1 (MQ=255)
           ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                            <  2:331697/93‑1 (MQ=255)
              tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                         >  1:1853729/1‑93 (MQ=255)
              tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                         <  1:4018898/93‑1 (MQ=255)
                ccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                                       >  1:3273075/1‑93 (MQ=255)
                 cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                      >  1:3253108/1‑93 (MQ=255)
                  gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                     >  2:1224717/1‑93 (MQ=255)
                  gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACAATCCCGag                                                                     >  2:3438152/1‑92 (MQ=255)
                    ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcgg                                                                   >  2:3233490/1‑93 (MQ=255)
                         gCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                                              >  1:1797520/1‑93 (MQ=255)
                           gTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                                            <  1:3347926/93‑1 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           >  1:428489/1‑93 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           <  2:3709931/93‑1 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           <  2:1719194/93‑1 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           >  1:1220855/1‑93 (MQ=255)
                             aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                          >  1:4469600/1‑93 (MQ=255)
                             aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                          <  1:854384/93‑1 (MQ=255)
                              tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                         <  1:4566552/93‑1 (MQ=255)
                              tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                         <  1:1138614/93‑1 (MQ=255)
                               cTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGcc                                                        <  1:926008/93‑1 (MQ=255)
                                tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                       <  1:330135/93‑1 (MQ=255)
                                 gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                      <  1:2401915/93‑1 (MQ=255)
                                  ggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggt                                                     <  1:4171729/93‑1 (MQ=255)
                                   gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                    <  2:1508311/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  <  2:817427/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  <  1:2588695/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  <  1:3092722/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  >  2:2226349/1‑93 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                <  2:4296306/93‑1 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                >  1:584703/1‑93 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                >  1:1422034/1‑93 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCGGGAACCGCGCAGCACCAACCGCCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                <  2:3937953/93‑1 (MQ=255)
                                        gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                               >  2:4068268/1‑93 (MQ=255)
                                        gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                               <  2:3288365/93‑1 (MQ=255)
                                         aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACt                                              <  1:2670748/93‑1 (MQ=255)
                                           gcgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAATCACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAc                                            <  2:1277327/93‑1 (MQ=255)
                                            cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                                           >  1:1587594/1‑93 (MQ=255)
                                             gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                          >  1:2887508/1‑93 (MQ=255)
                                             gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                          <  2:1543477/93‑1 (MQ=255)
                                             gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGATGGCCAAAGCCGGTGGTAACTACct                                          >  1:1557930/1‑93 (MQ=255)
                                              cTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                         <  1:768096/93‑1 (MQ=255)
                                               tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                        >  2:2465628/1‑93 (MQ=255)
                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  1:1053495/93‑1 (MQ=255)
                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  1:3964715/93‑1 (MQ=255)
                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  1:4642508/93‑1 (MQ=255)
                                                gaatcAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAACCACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  2:1769184/89‑1 (MQ=255)
                                                     aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                  >  2:2165645/1‑93 (MQ=255)
                                                     aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                  <  2:2890317/93‑1 (MQ=255)
                                                     aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                  <  1:3391756/93‑1 (MQ=255)
                                                      ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                 <  1:3061242/93‑1 (MQ=255)
                                                      ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                 >  2:2148768/1‑93 (MQ=255)
                                                       gggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                >  2:1862917/1‑93 (MQ=255)
                                                         gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                              <  1:1677842/93‑1 (MQ=255)
                                                         gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                              <  2:1819514/93‑1 (MQ=255)
                                                         gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGGGGTAACTACCTCTCTTCCctgct                              >  2:1871939/1‑93 (MQ=255)
                                                           tCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCTCTGCtgg                            >  1:3214348/1‑93 (MQ=255)
                                                           tCGATGCGATGGATTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCCGCtgg                            >  2:3574171/1‑93 (MQ=255)
                                                            cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                           <  2:3776091/93‑1 (MQ=255)
                                                            cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                           <  2:4215466/93‑1 (MQ=255)
                                                              atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                         >  1:81548/1‑93 (MQ=255)
                                                              atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                         >  1:3575207/1‑93 (MQ=255)
                                                                gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                       <  1:2473064/93‑1 (MQ=255)
                                                                gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                       >  1:3892885/1‑93 (MQ=255)
                                                                  gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                     <  2:1438595/93‑1 (MQ=255)
                                                                   atgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGc                    >  2:2698361/1‑93 (MQ=255)
                                                                   atgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAACTCTCTTCCCTGCTGGTGGGTAGc                    <  2:4519138/93‑1 (MQ=255)
                                                                     ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                  <  2:1015877/93‑1 (MQ=255)
                                                                     ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                  <  2:2769543/93‑1 (MQ=255)
                                                                      gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                 <  1:3104862/93‑1 (MQ=255)
                                                                      gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                 >  1:4313523/1‑93 (MQ=255)
                                                                      gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                 <  2:1171757/93‑1 (MQ=255)
                                                                       tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCACAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                >  1:2123091/1‑93 (MQ=255)
                                                                       tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                >  2:2022485/1‑93 (MQ=255)
                                                                          cctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgc             >  2:3894911/1‑93 (MQ=255)
                                                                            tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc           >  1:3500224/1‑93 (MQ=255)
                                                                            tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc           >  1:519012/1‑93 (MQ=255)
                                                                             cctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgcc          >  2:895485/1‑93 (MQ=255)
                                                                              ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg         >  2:132603/1‑93 (MQ=255)
                                                                              ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccc         >  2:3598298/1‑92 (MQ=255)
                                                                               tGGAGCCGCGCAGCACCAAACCCCATCCCGACGGCGGCAAAAGCCGGTTGTAACTGCCGCTCTTCCCGGCTGGTGGGTAGCGAAGCGcgccgc        <  2:722535/93‑1 (MQ=255)
                                                                               tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc        >  2:4032006/1‑93 (MQ=255)
                                                                                ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGGGAAGCGCgccgcc       >  2:1664976/1‑93 (MQ=255)
                                                                                ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCgccgcc       >  2:4182072/1‑93 (MQ=255)
                                                                                 gAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa      <  2:2376840/93‑1 (MQ=255)
                                                                                 gAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa      >  1:1351296/1‑93 (MQ=255)
                                                                                  aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc     <  1:4320229/93‑1 (MQ=255)
                                                                                  aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc     >  2:4128853/1‑93 (MQ=255)
                                                                                  aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc     <  1:3603246/93‑1 (MQ=255)
                                                                                   aCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACg    >  2:2289274/1‑93 (MQ=255)
                                                                                    ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg   >  1:3474187/1‑93 (MQ=255)
                                                                                    ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg   <  2:2589237/93‑1 (MQ=255)
                                                                                    ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg   <  2:3152105/93‑1 (MQ=255)
                                                                                           ||                                                                                     
GATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT  >  NC_000913/3950466‑3950643

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  >  NC_000913/3950469‑3950646
                                                                                        |                                                                                         
TATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCA                                                                                      >  SRR4302535.3656093/1‑94 (MQ=60)
 ATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAA                                                                                     >  SRR4302535.3451052/1‑94 (MQ=60)
 ATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAA                                                                                     >  SRR4302535.4504479/1‑94 (MQ=60)
  TCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAA                                                                                    <  SRR4302535.2713545/94‑1 (MQ=60)
   CGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAAC                                                                                   >  SRR4302535.3819259/1‑94 (MQ=60)
   CGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAAC                                                                                   <  SRR4302535.2428708/94‑1 (MQ=60)
   CGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAAC                                                                                   >  SRR4302535.1997284/1‑94 (MQ=60)
    GCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACA                                                                                  <  SRR4302535.2567071/94‑1 (MQ=60)
     CTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAC                                                                                 >  SRR4302535.4251170/1‑94 (MQ=60)
      TGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACC                                                                                <  SRR4302535.456467/94‑1 (MQ=60)
       GCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCA                                                                               <  SRR4302535.3795724/94‑1 (MQ=60)
       GCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCA                                                                               <  SRR4302535.2363702/94‑1 (MQ=60)
         TTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATC                                                                             <  SRR4302535.1408040/94‑1 (MQ=60)
         TTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATC                                                                             <  SRR4302535.3265569/94‑1 (MQ=60)
          TTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGGGCAGGGGATCGATGCGATGGTTTCTTCCTGGAACCGCGCGGCACCAAACACCATCC                                                                            >  SRR4302535.3810786/1‑94 (MQ=60)
          TTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCC                                                                            >  SRR4302535.1867081/1‑94 (MQ=60)
          TTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCC                                                                            <  SRR4302535.1078473/94‑1 (MQ=60)
           TCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCC                                                                           <  SRR4302535.4046396/94‑1 (MQ=60)
            CCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCG                                                                          >  SRR4302535.3295809/1‑94 (MQ=60)
             CCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGA                                                                         >  SRR4302535.3275711/1‑94 (MQ=60)
              CGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAC                                                                        <  SRR4302535.4346279/94‑1 (MQ=60)
               GTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACG                                                                       >  SRR4302535.1818173/1‑94 (MQ=60)
                 GGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGC                                                                     <  SRR4302535.567402/94‑1 (MQ=60)
                 GGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGC                                                                     >  SRR4302535.1546761/1‑94 (MQ=60)
                   GGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAAACGCGCAGCACCAAACACCATCCCGACGGCGG                                                                   >  SRR4302535.3645173/1‑94 (MQ=60)
                    GAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGC                                                                  <  SRR4302535.1755957/94‑1 (MQ=60)
                     AGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCA                                                                 >  SRR4302535.1810514/1‑94 (MQ=60)
                        GTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAA                                                              >  SRR4302535.431468/1‑94 (MQ=60)
                        GTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAA                                                              <  SRR4302535.3371247/94‑1 (MQ=60)
                        GTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAA                                                              >  SRR4302535.830455/1‑94 (MQ=60)
                        GTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAA                                                              >  SRR4302535.1229604/1‑94 (MQ=60)
                        GTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAA                                                              >  SRR4302535.4615538/1‑94 (MQ=60)
                         TATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAG                                                             >  SRR4302535.4500268/1‑94 (MQ=60)
                          ATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGC                                                            <  SRR4302535.860404/94‑1 (MQ=60)
                           TCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCC                                                           <  SRR4302535.4597703/94‑1 (MQ=60)
                           TCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCC                                                           <  SRR4302535.3612794/94‑1 (MQ=60)
                           TCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCC                                                           <  SRR4302535.1146731/94‑1 (MQ=60)
                            CTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCG                                                          <  SRR4302535.932569/94‑1 (MQ=60)
                             TGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGG                                                         <  SRR4302535.779280/94‑1 (MQ=60)
                             TGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGG                                                         <  SRR4302535.332422/94‑1 (MQ=60)
                              GGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGT                                                        <  SRR4302535.2419192/94‑1 (MQ=60)
                               GGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTG                                                       <  SRR4302535.3124859/94‑1 (MQ=60)
                               GGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTG                                                       <  SRR4302535.4200288/94‑1 (MQ=60)
                               GGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTG                                                       <  SRR4302535.4227156/94‑1 (MQ=60)
                                 CGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGT                                                     <  SRR4302535.1496347/94‑1 (MQ=60)
                                  GCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTA                                                    <  SRR4302535.3114257/94‑1 (MQ=60)
                                  GCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTA                                                    <  SRR4302535.2607286/94‑1 (MQ=60)
                                   CAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAA                                                   >  SRR4302535.588803/1‑94 (MQ=60)
                                   CAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAA                                                   >  SRR4302535.1432186/1‑94 (MQ=60)
                                    AGAAGCGCTGGAGCAGGGGATCGATGCGATGTTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  <  SRR4302535.3269100/94‑1 (MQ=60)
                                    AGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  >  SRR4302535.3559719/1‑94 (MQ=60)
                                     GAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACT                                                 <  SRR4302535.3313450/94‑1 (MQ=60)
                                      AAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTA                                                <  SRR4302535.2689857/94‑1 (MQ=60)
                                      AAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTA                                                >  SRR4302535.4060374/1‑94 (MQ=60)
                                       AGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAC                                               >  SRR4302535.535570/1‑94 (MQ=60)
                                        GCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACC                                              >  SRR4302535.1599051/1‑94 (MQ=60)
                                        GCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACC                                              <  SRR4302535.4686981/94‑1 (MQ=60)
                                         CGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCT                                             >  SRR4302535.2907938/1‑94 (MQ=60)
                                         CGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGATGGCCAAAGCCGGTGGTAACTACCT                                             >  SRR4302535.1569162/1‑94 (MQ=60)
                                           CTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCT                                           <  SRR4302535.773503/94‑1 (MQ=60)
                                            TGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTC                                          <  SRR4302535.3837439/94‑1 (MQ=60)
                                             GGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCT                                         <  SRR4302535.1573782/94‑1 (MQ=60)
                                             GGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCT                                         <  SRR4302535.3991799/94‑1 (MQ=60)
                                             GGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCT                                         <  SRR4302535.4673975/94‑1 (MQ=60)
                                             GGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCT                                         <  SRR4302535.1060987/94‑1 (MQ=60)
                                                GCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCC                                      >  SRR4302535.4018274/1‑94 (MQ=60)
                                                GCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCC                                      <  SRR4302535.4323966/94‑1 (MQ=60)
                                                 CAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCC                                     <  SRR4302535.3387636/94‑1 (MQ=60)
                                                  AGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCT                                    <  SRR4302535.3415354/94‑1 (MQ=60)
                                                  AGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCT                                    <  SRR4302535.3103707/94‑1 (MQ=60)
                                                  AGGAGATCGATGCGATGGTGTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCT                                    <  SRR4302535.659909/94‑1 (MQ=60)
                                                   GGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTG                                   <  SRR4302535.3082651/94‑1 (MQ=60)
                                                    GGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGC                                  <  SRR4302535.4269056/94‑1 (MQ=60)
                                                     GGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCT                                 <  SRR4302535.1605801/94‑1 (MQ=60)
                                                     GGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCT                                 >  SRR4302535.3771910/1‑94 (MQ=60)
                                                      GATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTG                                >  SRR4302535.3345047/1‑94 (MQ=60)
                                                      GATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTG                                <  SRR4302535.1689992/94‑1 (MQ=60)
                                                       ATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCTCTGCTGG                               >  SRR4302535.3236611/1‑94 (MQ=60)
                                                       ATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGG                               >  SRR4302535.3980143/1‑94 (MQ=60)
                                                       ATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGG                               <  SRR4302535.3054292/94‑1 (MQ=60)
                                                         CGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTG                             <  SRR4302535.962443/94‑1 (MQ=60)
                                                          GATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGG                            >  SRR4302535.3170072/1‑94 (MQ=60)
                                                          GATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGG                            >  SRR4302535.3484592/1‑94 (MQ=60)
                                                          GATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGG                            >  SRR4302535.3599994/1‑94 (MQ=60)
                                                          GATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGG                            >  SRR4302535.82127/1‑94 (MQ=60)
                                                           ATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCCAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGG                           >  SRR4302535.1167006/1‑94 (MQ=60)
                                                            TGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGT                          <  SRR4302535.3520005/94‑1 (MQ=60)
                                                            TGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGT                          >  SRR4302535.3919572/1‑94 (MQ=60)
                                                             GCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTA                         <  SRR4302535.2490901/94‑1 (MQ=60)
                                                             GCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTA                         >  SRR4302535.288172/1‑94 (MQ=60)
                                                              CGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAG                        >  SRR4302535.1002561/1‑94 (MQ=60)
                                                                ATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCG                      >  SRR4302535.4362459/1‑94 (MQ=60)
                                                                ATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCG                      <  SRR4302535.1348587/94‑1 (MQ=60)
                                                                 TGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGA                     >  SRR4302535.3286119/1‑94 (MQ=60)
                                                                  GGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAA                    >  SRR4302535.4343028/1‑94 (MQ=60)
                                                                   GTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCACAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   >  SRR4302535.2138393/1‑94 (MQ=60)
                                                                   GTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   <  SRR4302535.4127704/94‑1 (MQ=60)
                                                                   GTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   <  SRR4302535.3126479/94‑1 (MQ=60)
                                                                       CCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCG               >  SRR4302535.927719/1‑94 (MQ=60)
                                                                        CTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAACAGCCGGTGGTAACTACCTCTCTTCCCTGCTAGTGGGTAGCGAAGCGCGC              >  SRR4302535.2847364/1‑94 (MQ=60)
                                                                        CTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGC              >  SRR4302535.522644/1‑94 (MQ=60)
                                                                        CTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGC              >  SRR4302535.3524559/1‑94 (MQ=60)
                                                                        CTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGC              >  SRR4302535.4578898/1‑94 (MQ=60)
                                                                         TCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCC             >  SRR4302535.2634609/1‑94 (MQ=60)
                                                                          CCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCG            >  SRR4302535.4407329/1‑94 (MQ=60)
                                                                             GGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCA         >  SRR4302535.1360970/1‑94 (MQ=60)
                                                                               AACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACG       <  SRR4302535.3628249/94‑1 (MQ=60)
                                                                               AACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACG       <  SRR4302535.4349783/94‑1 (MQ=60)
                                                                                ACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGG      >  SRR4302535.3498326/1‑94 (MQ=60)
                                                                                 CCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT     <  SRR4302535.3775240/94‑1 (MQ=60)
                                                                                  CGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT    <  SRR4302535.1894030/94‑1 (MQ=60)
                                                                                    CGCAGCACCAAACACCCTCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  <  SRR4302535.1061846/94‑1 (MQ=60)
                                                                                    CGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  <  SRR4302535.1770246/94‑1 (MQ=60)
                                                                                    CGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  >  SRR4302535.3114782/1‑94 (MQ=60)
                                                                                        |                                                                                         
TATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  >  NC_000913/3950469‑3950646

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