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 1624 250.6 12620048 98.0% 12367647 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,5570CG99.6% 808.3 / ‑5.5 251R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (114/136);  minor base A (1/0);  total (115/136)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
*NC_0009133,950,5580GC100.0% 798.0 / NA 251R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (115/136);  total (115/136)

Only 100 of 255 total aligned reads displayed.

ATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT  >  NC_000913/3950470‑3950644
                                                                                       ||                                                                                      
  cGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCaaa                                                                                  <  2:2320814/93‑1 (MQ=255)
    ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                                                >  2:4946679/1‑93 (MQ=255)
     tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                               >  2:3308961/1‑93 (MQ=255)
       ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                             <  1:3733482/93‑1 (MQ=255)
       ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                             <  2:1654566/93‑1 (MQ=255)
       ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                             <  1:3256741/93‑1 (MQ=255)
        tttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                            <  1:685658/93‑1 (MQ=255)
        tttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                            <  2:4699247/93‑1 (MQ=255)
        tttCCCGTGGGAAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                            <  2:621568/93‑1 (MQ=255)
         ttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATc                                                                           <  2:3196203/93‑1 (MQ=255)
          tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          >  1:5440924/1‑93 (MQ=255)
          tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          <  2:1594026/93‑1 (MQ=255)
          tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          <  2:2375888/93‑1 (MQ=255)
          tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          <  1:2047341/93‑1 (MQ=255)
          tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          <  1:5719968/93‑1 (MQ=255)
          tCCCGTGGCGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          <  1:961490/93‑1 (MQ=255)
             cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                       >  1:2282867/1‑93 (MQ=255)
             cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                       <  2:590491/93‑1 (MQ=255)
              gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                      >  1:4504632/1‑93 (MQ=255)
              gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                      >  2:6177454/1‑93 (MQ=255)
                ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcgg                                                                    >  2:553316/1‑93 (MQ=255)
                  ggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                                  <  1:4035247/93‑1 (MQ=255)
                  ggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                                  <  1:5134627/93‑1 (MQ=255)
                    aGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACggcggc                                                                <  2:6118435/93‑1 (MQ=255)
                      cGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaa                                                              <  2:729452/93‑1 (MQ=255)
                        tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                            >  1:3824493/1‑93 (MQ=255)
                        tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                            >  2:906315/1‑93 (MQ=255)
                          tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                          <  2:12808/93‑1 (MQ=255)
                           cTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGcc                                                         >  1:2710494/1‑93 (MQ=255)
                             gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                       <  2:4837081/93‑1 (MQ=255)
                             gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGATTCCTCCTGGAAGCGCGCAACACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                       >  2:4677464/1‑93 (MQ=255)
                               gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAACCACCATCCCGACGGCGGCAAAAGCCggtg                                                     <  2:6262673/93‑1 (MQ=255)
                                 gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                   <  1:1004798/93‑1 (MQ=255)
                                 gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                   >  2:504455/1‑93 (MQ=255)
                                   aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                 <  1:65028/93‑1 (MQ=255)
                                    gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                                >  1:135507/1‑93 (MQ=255)
                                      aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTa                                              >  2:2263260/1‑93 (MQ=255)
                                       gcgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAc                                             <  2:1927876/93‑1 (MQ=255)
                                        cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                                            <  2:3552654/93‑1 (MQ=255)
                                         gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                           <  2:4257645/93‑1 (MQ=255)
                                          gtGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                          <  2:4149111/92‑1 (MQ=255)
                                          cTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                          <  1:1681923/93‑1 (MQ=255)
                                          cTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                          >  2:3839985/1‑93 (MQ=255)
                                           tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                         >  2:2920856/1‑93 (MQ=255)
                                           tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                         <  1:3422794/93‑1 (MQ=255)
                                           tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                         >  1:4344193/1‑93 (MQ=255)
                                           tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                         <  1:3519271/93‑1 (MQ=255)
                                            ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                        <  2:1868964/93‑1 (MQ=255)
                                             gAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctct                                       <  1:4131783/93‑1 (MQ=255)
                                              aGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCtt                                      <  1:1529613/93‑1 (MQ=255)
                                               gCAGGGGGTCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                                     >  2:5986327/1‑93 (MQ=255)
                                               gCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                                     <  1:4140782/93‑1 (MQ=255)
                                                 aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTGCCTCTCTTccc                                   <  1:2150944/93‑1 (MQ=255)
                                                 aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                   >  2:4360673/1‑93 (MQ=255)
                                                 aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                   <  2:4382472/93‑1 (MQ=255)
                                                  ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                  <  1:5928118/93‑1 (MQ=255)
                                                  ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                  <  2:4859687/93‑1 (MQ=255)
                                                   gggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTATCTCTCTTCCctg                                 >  2:172959/1‑93 (MQ=255)
                                                   gggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                 >  2:2946896/1‑93 (MQ=255)
                                                     gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                               <  1:1047376/93‑1 (MQ=255)
                                                      aTCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgctg                              >  2:3509302/1‑93 (MQ=255)
                                                       tCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtgg                             >  2:77764/1‑93 (MQ=255)
                                                       tCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtgg                             >  1:5233768/1‑93 (MQ=255)
                                                       tCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtgg                             <  1:1818880/93‑1 (MQ=255)
                                                        cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAGCCCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTCCCTCTCTTCCCTGCtggt                            <  1:2789103/93‑1 (MQ=255)
                                                        cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                            <  2:2792475/93‑1 (MQ=255)
                                                        cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                            <  1:3029210/93‑1 (MQ=255)
                                                        cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                            <  1:3162787/93‑1 (MQ=255)
                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                           <  2:449816/93‑1 (MQ=255)
                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                           <  1:4727586/93‑1 (MQ=255)
                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                           >  2:5659659/1‑93 (MQ=255)
                                                          atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                          <  2:3824210/93‑1 (MQ=255)
                                                           tgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTggg                         >  1:2887762/1‑93 (MQ=255)
                                                              gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                      >  1:1438264/1‑93 (MQ=255)
                                                               atgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGc                     <  2:4291401/93‑1 (MQ=255)
                                                                tgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCg                    >  1:4815521/1‑93 (MQ=255)
                                                                  gTTTCCTCCTGGACCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                  <  1:17569/93‑1 (MQ=255)
                                                                  gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                  <  1:6092340/93‑1 (MQ=255)
                                                                   tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                 >  2:1400543/1‑93 (MQ=255)
                                                                   tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                 <  1:2614427/93‑1 (MQ=255)
                                                                    ttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgc                <  1:6267375/93‑1 (MQ=255)
                                                                      cctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgc              >  2:5101253/1‑93 (MQ=255)
                                                                      cctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGCGGGTAGAGAAgggc              >  2:2494099/1‑93 (MQ=255)
                                                                       ctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcg             >  1:2424685/1‑93 (MQ=255)
                                                                        tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc            <  2:3990918/93‑1 (MQ=255)
                                                                        tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc            <  2:5169858/93‑1 (MQ=255)
                                                                        tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc            <  2:550262/93‑1 (MQ=255)
                                                                           tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc         <  1:4917050/93‑1 (MQ=255)
                                                                           tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc         >  1:3941281/1‑93 (MQ=255)
                                                                            ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCgccgcc        >  1:5687947/1‑93 (MQ=255)
                                                                            ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCgccgcc        >  2:2128349/1‑93 (MQ=255)
                                                                             gAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa       <  2:1054167/93‑1 (MQ=255)
                                                                               aCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCCCg     >  2:2979078/1‑93 (MQ=255)
                                                                               aCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACg     >  1:4273539/1‑93 (MQ=255)
                                                                               aCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACg     >  2:3788573/1‑93 (MQ=255)
                                                                                ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTCCCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg    <  1:5613942/93‑1 (MQ=255)
                                                                                ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg    <  1:5760753/93‑1 (MQ=255)
                                                                                 cgcgCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGt   <  1:5507898/93‑1 (MQ=255)
                                                                                  gcgcAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGtt  >  1:243763/1‑93 (MQ=255)
                                                                                  gcgcAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGtt  >  1:4700490/1‑93 (MQ=255)
                                                                                       ||                                                                                      
ATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT  >  NC_000913/3950470‑3950644

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
GTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  >  NC_000913/3950464‑3950646
                                                                                             |                                                                                         
GTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAG                                                                                           <  SRR4302533.5732814/94‑1 (MQ=60)
     TATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCA                                                                                      >  SRR4302533.2131594/1‑94 (MQ=60)
       TCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAA                                                                                    <  SRR4302533.1730112/94‑1 (MQ=60)
        CGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAAC                                                                                   >  SRR4302533.653123/1‑94 (MQ=60)
        CGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAAC                                                                                   >  SRR4302533.317959/1‑94 (MQ=60)
         GCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACA                                                                                  <  SRR4302533.4037595/94‑1 (MQ=60)
         GCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACA                                                                                  >  SRR4302533.3896907/1‑94 (MQ=60)
         GCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACA                                                                                  >  SRR4302533.5371263/1‑94 (MQ=60)
          CTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAC                                                                                 >  SRR4302533.1696071/1‑94 (MQ=60)
           TGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACC                                                                                <  SRR4302533.4326856/94‑1 (MQ=60)
            GCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACTA                                                                               <  SRR4302533.5981450/94‑1 (MQ=60)
            GCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCA                                                                               >  SRR4302533.5191138/1‑94 (MQ=60)
             CTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAT                                                                              <  SRR4302533.3276826/94‑1 (MQ=60)
             CTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAT                                                                              >  SRR4302533.324263/1‑94 (MQ=60)
             CTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAT                                                                              <  SRR4302533.5708660/94‑1 (MQ=60)
             CTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAT                                                                              <  SRR4302533.2206373/94‑1 (MQ=60)
             CTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAT                                                                              <  SRR4302533.3756384/94‑1 (MQ=60)
              TTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATC                                                                             <  SRR4302533.689713/94‑1 (MQ=60)
               TTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCC                                                                            >  SRR4302533.5473691/1‑94 (MQ=60)
                TCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCC                                                                           <  SRR4302533.3516208/94‑1 (MQ=60)
                TCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCC                                                                           <  SRR4302533.2059966/94‑1 (MQ=60)
                TCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCC                                                                           <  SRR4302533.5754506/94‑1 (MQ=60)
                TCCCGTGGCGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCC                                                                           <  SRR4302533.967137/94‑1 (MQ=60)
                 CCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCG                                                                          >  SRR4302533.5018318/1‑94 (MQ=60)
                  CCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGC                                                                         >  SRR4302533.464717/1‑94 (MQ=60)
                  CCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGA                                                                         >  SRR4302533.2296998/1‑94 (MQ=60)
                   CGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAC                                                                        >  SRR4302533.4531919/1‑94 (MQ=60)
                    GTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACG                                                                       >  SRR4302533.2599203/1‑94 (MQ=60)
                    GTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACG                                                                       >  SRR4302533.2744996/1‑94 (MQ=60)
                    GTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACG                                                                       >  SRR4302533.57618/1‑94 (MQ=60)
                     TGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGG                                                                      <  SRR4302533.1463935/94‑1 (MQ=60)
                     TGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGG                                                                      >  SRR4302533.3027434/1‑94 (MQ=60)
                     TGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGG                                                                      >  SRR4302533.3638921/1‑94 (MQ=60)
                      GGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGC                                                                     <  SRR4302533.1146063/94‑1 (MQ=60)
                        GGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGG                                                                   <  SRR4302533.4059811/94‑1 (MQ=60)
                        GGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGG                                                                   <  SRR4302533.5165649/94‑1 (MQ=60)
                             GTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAA                                                              >  SRR4302533.3847876/1‑94 (MQ=60)
                              TATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAG                                                             >  SRR4302533.5674025/1‑94 (MQ=60)
                              TATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAG                                                             >  SRR4302533.2854780/1‑94 (MQ=60)
                               ATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGC                                                            >  SRR4302533.4730774/1‑94 (MQ=60)
                               ATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGC                                                            <  SRR4302533.754001/94‑1 (MQ=60)
                                TCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCC                                                           >  SRR4302533.2727133/1‑94 (MQ=60)
                                 CTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCG                                                          <  SRR4302533.3835755/94‑1 (MQ=60)
                                   GGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGT                                                        <  SRR4302533.5478350/94‑1 (MQ=60)
                                   GGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGT                                                        >  SRR4302533.447065/1‑94 (MQ=60)
                                      CGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGCCGGCAAAAGCCGGTGGA                                                     >  SRR4302533.5123723/1‑94 (MQ=60)
                                       GCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTA                                                    <  SRR4302533.1010675/94‑1 (MQ=60)
                                       GCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTA                                                    >  SRR4302533.182343/1‑94 (MQ=60)
                                         AGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  >  SRR4302533.6216093/1‑94 (MQ=60)
                                         AGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  <  SRR4302533.65448/94‑1 (MQ=60)
                                         AGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  >  SRR4302533.136365/1‑94 (MQ=60)
                                         AGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  >  SRR4302533.5584334/1‑94 (MQ=60)
                                         AGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  <  SRR4302533.3944646/94‑1 (MQ=60)
                                         AGAAGCGCTGGAGCAGGGGATCGATGCGATGGGTTCATCCTGGAACCGCGCATCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAC                                                  >  SRR4302533.756344/1‑94 (MQ=60)
                                              CGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCT                                             >  SRR4302533.5593599/1‑94 (MQ=60)
                                              CGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCT                                             <  SRR4302533.3336039/94‑1 (MQ=60)
                                                CTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCT                                           >  SRR4302533.4370374/1‑94 (MQ=60)
                                                CTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCT                                           <  SRR4302533.1692110/94‑1 (MQ=60)
                                                 TGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTC                                          <  SRR4302533.3443937/94‑1 (MQ=60)
                                                 TGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTC                                          <  SRR4302533.3540967/94‑1 (MQ=60)
                                                 TGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTC                                          >  SRR4302533.931005/1‑94 (MQ=60)
                                                   GAGGAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTT                                        <  SRR4302533.4156864/94‑1 (MQ=60)
                                                   GAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTT                                        <  SRR4302533.4156863/94‑1 (MQ=60)
                                                   GAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTT                                        <  SRR4302533.3191416/94‑1 (MQ=60)
                                                    AGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTC                                       <  SRR4302533.1538811/94‑1 (MQ=60)
                                                     GCAGGGGATCGGTGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACTCCAGCCCGACGGCGGCAAAAGCCGGTGGTAACTGCCTCTCTTCC                                      <  SRR4302533.1418869/94‑1 (MQ=60)
                                                     GCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCC                                      <  SRR4302533.4165895/94‑1 (MQ=60)
                                                     GCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCC                                      <  SRR4302533.420481/94‑1 (MQ=60)
                                                     GCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCC                                      <  SRR4302533.3297057/94‑1 (MQ=60)
                                                      CAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCC                                     <  SRR4302533.4384201/94‑1 (MQ=60)
                                                      CAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCC                                     <  SRR4302533.4639325/94‑1 (MQ=60)
                                                      CAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCC                                     <  SRR4302533.1615388/94‑1 (MQ=60)
                                                       AGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTGCCTCTCTTCCCT                                    <  SRR4302533.2164189/94‑1 (MQ=60)
                                                       AGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCT                                    >  SRR4302533.979212/1‑94 (MQ=60)
                                                       AGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCT                                    <  SRR4302533.2754244/94‑1 (MQ=60)
                                                        GGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTG                                   >  SRR4302533.2981725/1‑94 (MQ=60)
                                                        GGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTG                                   <  SRR4302533.5963952/94‑1 (MQ=60)
                                                         GGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGC                                  <  SRR4302533.2281988/94‑1 (MQ=60)
                                                           GATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTG                                <  SRR4302533.1053531/94‑1 (MQ=60)
                                                            ATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGG                               >  SRR4302533.5265242/1‑94 (MQ=60)
                                                             TCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTCCCTCTCTTCCCTGCTGGT                              <  SRR4302533.2451651/94‑1 (MQ=60)
                                                             TCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGT                              <  SRR4302533.4099682/94‑1 (MQ=60)
                                                             TCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGT                              <  SRR4302533.4098506/94‑1 (MQ=60)
                                                             TCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGT                              <  SRR4302533.1829990/94‑1 (MQ=60)
                                                             TCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGT                              <  SRR4302533.5871712/94‑1 (MQ=60)
                                                              CGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAGCCCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTCCCTCTCTTCCCTGCTGGTG                             <  SRR4302533.2806249/94‑1 (MQ=60)
                                                              CGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTG                             <  SRR4302533.3047746/94‑1 (MQ=60)
                                                              CGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTG                             <  SRR4302533.3182256/94‑1 (MQ=60)
                                                               GATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGG                            <  SRR4302533.4756247/94‑1 (MQ=60)
                                                                ATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGG                           >  SRR4302533.2019306/1‑94 (MQ=60)
                                                                ATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGG                           >  SRR4302533.2905514/1‑94 (MQ=60)
                                                                ATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGG                           <  SRR4302533.1478557/94‑1 (MQ=60)
                                                                 TGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGTTGGTAACTACCTCTCTTCCCTGCTGGTGGGT                          <  SRR4302533.5968363/94‑1 (MQ=60)
                                                                   CGATGGTTTCCTCCTGGAACCGNGCNGCACCAAACACNATCCCGANGGCGGNGANAGCCGGTGNTAACTANNNCNCTTCCNTNCTNGNNNGNAG                        >  SRR4302533.1222835/1‑94 (MQ=42)
                                                                   CGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAG                        <  SRR4302533.674190/94‑1 (MQ=60)
                                                                   CGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAG                        >  SRR4302533.640627/1‑94 (MQ=60)
                                                                   CGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAG                        >  SRR4302533.1446945/1‑94 (MQ=60)
                                                                   CGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAG                        >  SRR4302533.5746337/1‑94 (MQ=60)
                                                                    GATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGC                       >  SRR4302533.4764586/1‑94 (MQ=60)
                                                                     ATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCG                      >  SRR4302533.4844707/1‑94 (MQ=60)
                                                                      TGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGA                     >  SRR4302533.3619258/1‑94 (MQ=60)
                                                                        GTTTCCTCCTGGACCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   <  SRR4302533.17690/94‑1 (MQ=60)
                                                                        GTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   >  SRR4302533.4282738/1‑94 (MQ=60)
                                                                        GTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   <  SRR4302533.6129156/94‑1 (MQ=60)
                                                                        GTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAG                   >  SRR4302533.5791805/1‑94 (MQ=60)
                                                                         TTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGC                  <  SRR4302533.2630476/94‑1 (MQ=60)
                                                                          TTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCG                 >  SRR4302533.2368178/1‑94 (MQ=60)
                                                                          TTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCG                 <  SRR4302533.6304823/94‑1 (MQ=60)
                                                                          TTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCG                 <  SRR4302533.2709389/94‑1 (MQ=60)
                                                                          TGCCTCCTGGAACCGCGCAGCACCAACCACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCG                 <  SRR4302533.2338414/94‑1 (MQ=60)
                                                                            CCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCG               >  SRR4302533.2439655/1‑94 (MQ=60)
                                                                            CCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCG               <  SRR4302533.2459589/94‑1 (MQ=60)
                                                                            CCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCG               >  SRR4302533.4927590/1‑94 (MQ=60)
                                                                              TCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCC             >  SRR4302533.1215277/1‑94 (MQ=60)
                                                                               CCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCG            >  SRR4302533.2048784/1‑94 (MQ=60)
                                                                                CTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGC           >  SRR4302533.3965330/1‑94 (MQ=60)
                                                                                 TGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCC          >  SRR4302533.5336164/1‑94 (MQ=60)
                                                                                 TGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCC          >  SRR4302533.5722301/1‑94 (MQ=60)
                                                                                 TGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCC          <  SRR4302533.4946832/94‑1 (MQ=60)
                                                                                    AACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACG       >  SRR4302533.4299254/1‑94 (MQ=60)
                                                                                      CCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTCCCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT     <  SRR4302533.5647874/94‑1 (MQ=60)
                                                                                      CCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT     <  SRR4302533.5795513/94‑1 (MQ=60)
                                                                                      CCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT     <  SRR4302533.4771620/94‑1 (MQ=60)
                                                                                       CGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT    <  SRR4302533.1032225/94‑1 (MQ=60)
                                                                                       CGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT    >  SRR4302533.245275/1‑94 (MQ=60)
                                                                                       CGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT    <  SRR4302533.5541142/94‑1 (MQ=60)
                                                                                       CGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT    >  SRR4302533.4728990/1‑94 (MQ=60)
                                                                                        GCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTA   >  SRR4302533.3418416/1‑94 (MQ=60)
                                                                                         CGCAGCACCAAACCCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  <  SRR4302533.4654289/94‑1 (MQ=60)
                                                                                             |                                                                                         
GTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAT  >  NC_000913/3950464‑3950646

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