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,183,028 T→C V70A (GTC→GCC)  b3042 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,0280TC99.3% 443.6 / ‑4.7 144V70A (GTC→GCC) b3042hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (71/72);  total (71/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

Only 100 of 150 total aligned reads displayed.

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG  >  NC_000913/3182938‑3183117
                                                                                          |                                                                                         
aTGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGccc                                                                                         <  1:4460297/93‑1 (MQ=255)
  gtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTg                                                                                       >  2:4134545/1‑93 (MQ=255)
   tgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCTAGCCCTGt                                                                                      <  2:4571626/93‑1 (MQ=255)
   tgGAGAAAAAAACCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTGGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGt                                                                                      <  2:2997399/93‑1 (MQ=255)
     aagaAAAAAACCCGCCAACCCCTACAAGCGCAGCGGCCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTa                                                                                    <  2:2654008/92‑1 (MQ=255)
      agaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTAc                                                                                   >  2:3274848/1‑93 (MQ=255)
      agaAAAAAATCCGCCAAACCCGACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTAc                                                                                   <  2:387371/93‑1 (MQ=255)
       gaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                  >  1:2498179/1‑93 (MQ=255)
       aaaaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                  <  2:4316901/92‑1 (MQ=255)
        aaaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACgg                                                                                 >  1:1993557/1‑93 (MQ=255)
          aaaaaCCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGAc                                                                               <  1:226949/93‑1 (MQ=255)
           aaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACt                                                                              >  2:3719788/1‑93 (MQ=255)
            aaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTc                                                                             >  1:4350651/1‑93 (MQ=255)
             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                            >  2:4195733/1‑93 (MQ=255)
             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                            >  2:1864467/1‑93 (MQ=255)
             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                            >  2:4002684/1‑93 (MQ=255)
              acccGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGt                                                                           <  2:946955/91‑1 (MQ=255)
              aTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGt                                                                           <  1:3027177/93‑1 (MQ=255)
               tCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTg                                                                          >  1:1553804/1‑93 (MQ=255)
                  gCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaa                                                                       >  2:3814706/1‑93 (MQ=255)
                   ccAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaa                                                                      >  1:1889257/1‑93 (MQ=255)
                       ccccTGCAGGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTg                                                                  <  2:3823100/92‑1 (MQ=255)
                       aCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTg                                                                  >  1:1232838/1‑93 (MQ=255)
                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                 >  2:3069471/1‑93 (MQ=255)
                        cccAACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                 <  1:56177/93‑1 (MQ=255)
                          cTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcg                                                               >  1:667321/1‑93 (MQ=255)
                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  2:400458/1‑93 (MQ=255)
                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              <  2:3687163/93‑1 (MQ=255)
                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  2:3606176/1‑93 (MQ=255)
                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  1:4137486/1‑93 (MQ=255)
                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              <  1:3747860/93‑1 (MQ=255)
                            aCAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgct                                                             >  2:1606505/1‑93 (MQ=255)
                            aCAAGCCCAGCGGCCGCGCGTCGACCTTGTAAGCCGGGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgct                                                             <  1:816972/93‑1 (MQ=255)
                              aaGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgc                                                           <  1:4230609/93‑1 (MQ=255)
                               aGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgct                                                          >  2:1967684/1‑93 (MQ=255)
                                gcgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                         <  1:465613/93‑1 (MQ=255)
                                gcgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                         <  1:3812273/93‑1 (MQ=255)
                                gcgcAGCTGACGCGCCTCGACCTTGTAACCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                         <  1:1122040/93‑1 (MQ=255)
                                  gcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa                                                       >  2:4501349/1‑93 (MQ=255)
                                  gcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa                                                       >  2:2721080/1‑93 (MQ=255)
                                   cAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGaa                                                      <  1:2389363/93‑1 (MQ=255)
                                    cgCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  1:4179204/92‑1 (MQ=255)
                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  1:1609513/93‑1 (MQ=255)
                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  1:2623214/93‑1 (MQ=255)
                                      cTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa                                                   <  2:1709391/93‑1 (MQ=255)
                                       tGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa                                                  >  1:1017413/1‑93 (MQ=255)
                                        gtcgcgCCTCGATCGCTTAAGCCGTGAAGATTACGACGTCGAACCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAg                                                 <  2:2389788/91‑1 (MQ=255)
                                        cgcgcgCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAg                                                 <  2:2516453/91‑1 (MQ=255)
                                         aCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa                                                <  2:450218/93‑1 (MQ=255)
                                         aCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa                                                <  2:2295232/93‑1 (MQ=255)
                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  1:2058589/93‑1 (MQ=255)
                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  2:4304002/93‑1 (MQ=255)
                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  1:3961765/93‑1 (MQ=255)
                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCCACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  2:787931/93‑1 (MQ=255)
                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAGAACTGGCGCTGCTTGAACAAAGAATCa                                            <  1:1827335/93‑1 (MQ=255)
                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                            <  2:3040364/93‑1 (MQ=255)
                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                            <  1:1556472/93‑1 (MQ=255)
                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                            <  2:736203/93‑1 (MQ=255)
                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                           >  2:1799235/1‑93 (MQ=255)
                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                           <  1:3940785/93‑1 (MQ=255)
                                               cTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGc                                          >  1:1719821/1‑93 (MQ=255)
                                                tCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCg                                         <  1:4458578/93‑1 (MQ=255)
                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      <  1:1800969/93‑1 (MQ=255)
                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      >  2:3257795/1‑93 (MQ=255)
                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      >  2:4044592/1‑93 (MQ=255)
                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  2:1929021/1‑93 (MQ=255)
                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  1:2911093/1‑93 (MQ=255)
                                                        gTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaa                                 >  2:2973065/1‑93 (MQ=255)
                                                        gTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaa                                 >  1:1399459/1‑93 (MQ=255)
                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                <  1:104166/93‑1 (MQ=255)
                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                <  2:3398766/93‑1 (MQ=255)
                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                >  1:4556417/1‑93 (MQ=255)
                                                           aGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAAt                              <  2:93006/93‑1 (MQ=255)
                                                            gCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                             <  1:36355/93‑1 (MQ=255)
                                                            gCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                             <  2:1908811/93‑1 (MQ=255)
                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:3493684/93‑1 (MQ=255)
                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  1:2609747/93‑1 (MQ=255)
                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTAGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:774556/93‑1 (MQ=255)
                                                                 gAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGt                        >  1:1254590/1‑93 (MQ=255)
                                                                  aaGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTa                       >  1:3120133/1‑93 (MQ=255)
                                                                  aaGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTa                       <  1:3274848/93‑1 (MQ=255)
                                                                  aaGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTa                       >  2:243142/1‑93 (MQ=255)
                                                                   agagTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTAc                      >  1:2764545/1‑93 (MQ=255)
                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                    <  2:1170332/93‑1 (MQ=255)
                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                    >  2:472870/1‑93 (MQ=255)
                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                    <  2:2160598/93‑1 (MQ=255)
                                                                      gTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                   <  2:1358399/93‑1 (MQ=255)
                                                                      attCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                   <  2:4178880/92‑1 (MQ=255)
                                                                        tCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgaga                 >  2:1714800/1‑93 (MQ=255)
                                                                         cGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGAt                <  1:3363374/93‑1 (MQ=255)
                                                                          gACGTCCAAACGCAAGCCCTGTTACGGACGCGTGAAAAACTGGCGCTGCGTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATc               <  1:1610645/93‑1 (MQ=255)
                                                                              tCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaa           <  2:150349/93‑1 (MQ=255)
                                                                                 aaaCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagc        <  2:2712406/93‑1 (MQ=255)
                                                                                  aaCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagca       >  1:316068/1‑93 (MQ=255)
                                                                                   aCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagcag      <  2:3745157/93‑1 (MQ=255)
                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc    >  2:4251741/1‑93 (MQ=255)
                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc    >  2:198624/1‑93 (MQ=255)
                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc    <  2:693449/93‑1 (MQ=255)
                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCa   <  1:2323720/93‑1 (MQ=255)
                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCa   >  2:1394591/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG  >  NC_000913/3182938‑3183117

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGATC  >  NC_000913/3182936‑3183120
                                                                                            |                                                                                            
AGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCC                                                                                             >  SRR4302535.2552316/1‑94 (MQ=60)
  ATGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCT                                                                                           <  SRR4302535.4490889/94‑1 (MQ=60)
    GTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGT                                                                                         <  SRR4302535.3007287/94‑1 (MQ=60)
     GGGAGAAAAAAACCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTT                                                                                        <  SRR4302535.4244568/94‑1 (MQ=60)
        AGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACG                                                                                     >  SRR4302535.2516202/1‑94 (MQ=60)
        AGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACG                                                                                     <  SRR4302535.3360117/94‑1 (MQ=60)
         GAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGG                                                                                    >  SRR4302535.2007912/1‑94 (MQ=60)
         GAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGCCCAAACGCAAGCCCTGTTACGG                                                                                    >  SRR4302535.1584903/1‑94 (MQ=60)
           AAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGAC                                                                                  >  SRR4302535.2383460/1‑94 (MQ=60)
            AAAAACCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACT                                                                                 <  SRR4302535.228511/94‑1 (MQ=60)
             AAAATCCGCCAAACCCTACAAGCGCCGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTG                                                                                >  SRR4302535.722858/1‑94 (MQ=60)
             AAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTC                                                                                >  SRR4302535.4380418/1‑94 (MQ=60)
                ATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTG                                                                             >  SRR4302535.1565000/1‑94 (MQ=60)
                ATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTG                                                                             <  SRR4302535.3048370/94‑1 (MQ=60)
                   CGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAA                                                                          >  SRR4302535.1072763/1‑94 (MQ=60)
                    GCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAA                                                                         >  SRR4302535.1902865/1‑94 (MQ=60)
                     CCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAA                                                                        <  SRR4302535.1815080/94‑1 (MQ=60)
                       AAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACT                                                                      >  SRR4302535.2755224/1‑94 (MQ=60)
                        AACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTG                                                                     >  SRR4302535.4105116/1‑94 (MQ=60)
                        AACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTG                                                                     >  SRR4302535.1241655/1‑94 (MQ=60)
                          CCCAACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGC                                                                   <  SRR4302535.56584/94‑1 (MQ=60)
                           CCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCG                                                                  >  SRR4302535.4550143/1‑94 (MQ=60)
                           CCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCG                                                                  >  SRR4302535.671972/1‑94 (MQ=60)
                            CTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGC                                                                 >  SRR4302535.4165818/1‑94 (MQ=60)
                             TACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCT                                                                <  SRR4302535.3773788/94‑1 (MQ=60)
                              ACAAGCCCAGCGGCCGCGCGTCGACCTTGTAAGCCGGGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTG                                                               <  SRR4302535.822702/94‑1 (MQ=60)
                                AAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCT                                                             <  SRR4302535.4259559/94‑1 (MQ=60)
                                  GCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTG                                                           <  SRR4302535.468869/94‑1 (MQ=60)
                                  GCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTG                                                           <  SRR4302535.3838581/94‑1 (MQ=60)
                                  GCGCAGCTGACGCGCCTCGACCTTGTAACCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTG                                                           <  SRR4302535.1130046/94‑1 (MQ=60)
                                     CAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAC                                                        <  SRR4302535.2406540/94‑1 (MQ=60)
                                     CAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAC                                                        >  SRR4302535.16925/1‑94 (MQ=60)
                                      CGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACA                                                       <  SRR4302535.4207809/94‑1 (MQ=60)
                                      AGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACA                                                       <  SRR4302535.2217579/94‑1 (MQ=60)
                                      AGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACA                                                       <  SRR4302535.1621118/94‑1 (MQ=60)
                                      AGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACA                                                       <  SRR4302535.2642013/94‑1 (MQ=60)
                                        CTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAA                                                     >  SRR4302535.1024644/1‑94 (MQ=60)
                                              CGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCA                                               <  SRR4302535.3988821/94‑1 (MQ=60)
                                              CGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCA                                               <  SRR4302535.2073419/94‑1 (MQ=60)
                                               GCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAGAACTGGCGCTGCTTGAACAAAGAATCAG                                              <  SRR4302535.1840527/94‑1 (MQ=60)
                                               GCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAG                                              >  SRR4302535.3170606/1‑94 (MQ=60)
                                               GCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAG                                              <  SRR4302535.1567692/94‑1 (MQ=60)
                                                CCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGC                                             <  SRR4302535.3967687/94‑1 (MQ=60)
                                                CCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGC                                             >  SRR4302535.1732261/1‑94 (MQ=60)
                                                 CTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCG                                            >  SRR4302535.4399629/1‑94 (MQ=60)
                                                  TCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGA                                           <  SRR4302535.4489157/94‑1 (MQ=60)
                                                  TCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGA                                           <  SRR4302535.1256882/94‑1 (MQ=60)
                                                     ACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACT                                        <  SRR4302535.1813987/94‑1 (MQ=60)
                                                     ACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACT                                        >  SRR4302535.2931654/1‑94 (MQ=60)
                                                         TGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAA                                    >  SRR4302535.1409449/1‑94 (MQ=60)
                                                          GTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAA                                   >  SRR4302535.4587532/1‑94 (MQ=60)
                                                           TAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAA                                  <  SRR4302535.104909/94‑1 (MQ=60)
                                                            AAGCCGGTAAGAGTTCGACGGCCAAACCGAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAAT                                 <  SRR4302535.1528601/94‑1 (MQ=60)
                                                             AGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATC                                >  SRR4302535.525407/1‑94 (MQ=60)
                                                              GCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCG                               <  SRR4302535.36637/94‑1 (MQ=60)
                                                                 GTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAG                            <  SRR4302535.2628474/94‑1 (MQ=60)
                                                                  TGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGT                           >  SRR4302535.1263562/1‑94 (MQ=60)
                                                                  TGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGT                           >  SRR4302535.3893568/1‑94 (MQ=60)
                                                                   GAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTA                          <  SRR4302535.19038/94‑1 (MQ=60)
                                                                   GAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTA                          >  SRR4302535.3700866/1‑94 (MQ=60)
                                                                   GAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTA                          >  SRR4302535.3141857/1‑94 (MQ=60)
                                                                    AAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTAC                         >  SRR4302535.2784237/1‑94 (MQ=60)
                                                                    AAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTAC                         <  SRR4302535.3297593/94‑1 (MQ=60)
                                                                           CGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATC                  <  SRR4302535.3386777/94‑1 (MQ=60)
                                                                            GACGTCCAAACGCAAGCCCTGTTACGGACGCGTGAAAAACTGGCGCTGCGTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCA                 <  SRR4302535.1622255/94‑1 (MQ=60)
                                                                                   AAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCA          >  SRR4302535.318262/1‑94 (MQ=60)
                                                                                      CGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGACTCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCC       <  SRR4302535.2531096/94‑1 (MQ=60)
                                                                                        CAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG     <  SRR4302535.2340347/94‑1 (MQ=60)
                                                                                          AGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGAT   <  SRR4302535.3969805/94‑1 (MQ=60)
                                                                                          AGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGAT   <  SRR4302535.502540/94‑1 (MQ=60)
                                                                                          AGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGAT   <  SRR4302535.2652815/94‑1 (MQ=60)
                                                                                          AGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGAT   >  SRR4302535.17698/1‑94 (MQ=60)
                                                                                           GCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGATC  <  SRR4302535.2081347/94‑1 (MQ=60)
                                                                                            |                                                                                            
AGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAGATC  >  NC_000913/3182936‑3183120

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