Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,577,071 A→G *114Q (TAA→CAA)  yjiW ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,577,0710AG99.5% 737.6 / ‑5.2 209*114Q (TAA→CAA) yjiWhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (94/114);  total (95/114)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 225 total aligned reads displayed.

AAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG  >  NC_000913/4576981‑4577160
                                                                                          |                                                                                         
aaaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGcg                                                                                         <  2:1670625/93‑1 (MQ=255)
aaaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACGGGAGTAAACGCTGTTGcg                                                                                         >  1:1087940/1‑93 (MQ=255)
 aaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTgcgc                                                                                        >  1:5393023/1‑93 (MQ=255)
   aTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGa                                                                                      <  1:3421068/93‑1 (MQ=255)
    acAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGAc                                                                                     <  2:4866083/92‑1 (MQ=255)
      aGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGGGCGACtt                                                                                   >  1:4675376/1‑93 (MQ=255)
      aGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACtt                                                                                   <  1:5140967/93‑1 (MQ=255)
        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTTTTGCGCGACTTTc                                                                                 <  1:4843478/93‑1 (MQ=255)
        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACtttt                                                                                 >  2:5374535/1‑92 (MQ=255)
        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTc                                                                                 <  1:3084810/93‑1 (MQ=255)
        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGCGGTAAAAGCTGTTGCGCGACATTc                                                                                 >  2:1815115/1‑93 (MQ=255)
         cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCt                                                                                >  1:1956165/1‑93 (MQ=255)
         cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCt                                                                                >  2:4156267/1‑93 (MQ=255)
          tGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTg                                                                               <  1:209597/93‑1 (MQ=255)
           ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGt                                                                              >  2:5238530/1‑93 (MQ=255)
           ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGGAAAAGCTGTTGCGCGACTTTCTGt                                                                              >  2:323894/1‑93 (MQ=255)
            gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGAGAGTAAACGCTGTTGCGCGACTTTCTGtt                                                                             <  2:4273442/93‑1 (MQ=255)
               aaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGGTCCCGGTGGTTAACCGGGAGGAAACGCCGTTGCGCGACTTTCTGTTTAc                                                                          >  2:3521822/1‑93 (MQ=255)
                atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         >  1:4058235/1‑93 (MQ=255)
                atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         <  1:5268423/93‑1 (MQ=255)
                atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         >  2:5711354/1‑93 (MQ=255)
                 tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTGACCGGGAGTAAAAGCCGTTTCGCGACTTTTTGTTTAACg                                                                        >  1:3514737/1‑93 (MQ=255)
                 tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCg                                                                        >  2:5849284/1‑93 (MQ=255)
                  atTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAACgg                                                                       >  2:3871074/1‑93 (MQ=255)
                  atTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGGTAACCGGGAGGAAACGCTGGTGCGCGACGTTCTGTTTACCgg                                                                       >  2:237075/1‑93 (MQ=255)
                   tttACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAAAGCTGTTGCGCGACTTTCTGTTTACCGGc                                                                      >  2:3018250/1‑93 (MQ=255)
                     tACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCaa                                                                    >  1:5002514/1‑93 (MQ=255)
                       cTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATc                                                                  <  2:885080/93‑1 (MQ=255)
                       cTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATc                                                                  <  1:6138310/93‑1 (MQ=255)
                        tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAAttc                                                                 >  1:3376789/1‑91 (MQ=255)
                        tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCa                                                                 <  2:2914341/93‑1 (MQ=255)
                        tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCa                                                                 <  1:5727249/93‑1 (MQ=255)
                            ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTcc                                                             <  1:2939248/93‑1 (MQ=255)
                             gggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGCCTTTCTGTTTACCGGCAATCACTCCa                                                            >  2:902807/1‑93 (MQ=255)
                               gAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAAt                                                          <  1:2371269/93‑1 (MQ=255)
                                aaaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATa                                                         <  1:3599628/93‑1 (MQ=255)
                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGCCTTTCTGTTTACCGGCAATCACTCCAATaa                                                        >  1:298483/1‑93 (MQ=255)
                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCAATCAATCCAATaa                                                        >  1:1468030/1‑93 (MQ=255)
                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                        >  2:2465032/1‑93 (MQ=255)
                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                        <  1:4901117/93‑1 (MQ=255)
                                  aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       <  2:2258821/93‑1 (MQ=255)
                                  aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       <  1:2564316/93‑1 (MQ=255)
                                     ttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGc                                                    <  1:1242288/93‑1 (MQ=255)
                                      tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCCATAAACGcc                                                   >  1:5840047/1‑93 (MQ=255)
                                      tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGcc                                                   >  2:3586747/1‑93 (MQ=255)
                                       cTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCt                                                  >  1:2176767/1‑93 (MQ=255)
                                        tATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTg                                                 >  1:1881239/1‑93 (MQ=255)
                                         aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                <  2:4363637/93‑1 (MQ=255)
                                         aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                <  2:2597133/93‑1 (MQ=255)
                                         aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                >  1:4243976/1‑93 (MQ=255)
                                         aTTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                >  1:1516558/1‑93 (MQ=255)
                                          ttGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCa                                               >  2:402240/1‑93 (MQ=255)
                                          ttGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCa                                               >  1:942427/1‑93 (MQ=255)
                                           tGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAc                                              <  1:4951362/93‑1 (MQ=255)
                                           tGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAc                                              <  2:5077812/93‑1 (MQ=255)
                                            gAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAcc                                             <  1:1430193/93‑1 (MQ=255)
                                                 aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAACCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCtt                                        <  2:269865/93‑1 (MQ=255)
                                                 aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTGCCGGCAATCACTCCAATAAACGCCTGCACCTGCtt                                        <  1:1615996/93‑1 (MQ=255)
                                                  gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCCCCTGCttt                                       >  2:5363572/1‑93 (MQ=255)
                                                   tGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCtttt                                      >  1:3870787/1‑93 (MQ=255)
                                                      aaGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGtt                                   <  1:4371810/93‑1 (MQ=255)
                                                        ggCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 <  2:4988357/93‑1 (MQ=255)
                                                        ggCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 >  1:3619010/1‑93 (MQ=255)
                                                        ggCTCCCGGTGGTTAACCGGGAGTAAAAGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 >  2:6042114/1‑93 (MQ=255)
                                                          cTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcg                               <  2:4687551/93‑1 (MQ=255)
                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              <  1:5316389/93‑1 (MQ=255)
                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              >  2:2945488/1‑93 (MQ=255)
                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              >  1:1746164/1‑93 (MQ=255)
                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              <  2:771252/93‑1 (MQ=255)
                                                            cccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgcg                             >  2:3817112/1‑93 (MQ=255)
                                                            cccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgcg                             >  2:4507367/1‑93 (MQ=255)
                                                             ccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACgcgc                            >  2:5596274/1‑93 (MQ=255)
                                                              cGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGcc                           <  2:2794024/93‑1 (MQ=255)
                                                               ggtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                          >  1:4030039/1‑93 (MQ=255)
                                                                gtggtTAAGCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGa                         <  1:2791732/93‑1 (MQ=255)
                                                                 tggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGAc                        <  2:3154001/93‑1 (MQ=255)
                                                                  ggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACa                       >  1:1834482/1‑93 (MQ=255)
                                                                    ttAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGt                     >  1:5546365/1‑93 (MQ=255)
                                                                    ttAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGt                     <  1:4505726/93‑1 (MQ=255)
                                                                     tAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                    <  1:75708/93‑1 (MQ=255)
                                                                      aaCCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGttt                   <  2:4403941/93‑1 (MQ=255)
                                                                         cGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGca                <  2:193093/93‑1 (MQ=255)
                                                                          gggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcac               >  2:1895305/1‑93 (MQ=255)
                                                                           ggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca              <  1:479262/93‑1 (MQ=255)
                                                                           ggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca              <  2:1471362/93‑1 (MQ=255)
                                                                             aGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACAcc            <  1:4861099/93‑1 (MQ=255)
                                                                             aGTAAACCCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAATCGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCATAcc            <  2:1333969/93‑1 (MQ=255)
                                                                              gTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCt           <  2:5262248/93‑1 (MQ=255)
                                                                              gTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCt           <  1:4919243/93‑1 (MQ=255)
                                                                               tAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTg          <  1:5959537/93‑1 (MQ=255)
                                                                                  aCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCg       >  1:3835744/1‑93 (MQ=255)
                                                                                  aCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCg       >  2:1752071/1‑93 (MQ=255)
                                                                                   cGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGt      >  2:272219/1‑93 (MQ=255)
                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     <  1:2156440/93‑1 (MQ=255)
                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     <  2:3662910/93‑1 (MQ=255)
                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     <  1:187432/93‑1 (MQ=255)
                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     >  2:118375/1‑93 (MQ=255)
                                                                                     cTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAg    >  1:1858099/1‑93 (MQ=255)
                                                                                      tGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAg  >  2:3076167/1‑92 (MQ=255)
                                                                                      tGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGc   >  2:5806786/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG  >  NC_000913/4576981‑4577160

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGA  >  NC_000913/4576978‑4577161
                                                                                             |                                                                                          
TGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTG                                                                                            <  SRR4302532.3918475/94‑1 (MQ=60)
TGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTG                                                                                            <  SRR4302532.4507632/94‑1 (MQ=60)
 GCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGC                                                                                           <  SRR4302532.599550/94‑1 (MQ=60)
  CAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCG                                                                                          <  SRR4302532.2165548/94‑1 (MQ=60)
  CAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACGGGAGTAAACGCTGTTGCG                                                                                          >  SRR4302532.1124373/1‑94 (MQ=60)
   AAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTGCGC                                                                                         >  SRR4302532.5579868/1‑94 (MQ=60)
      ATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGAC                                                                                      <  SRR4302532.3535669/94‑1 (MQ=60)
        CAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGGGCGACTT                                                                                    >  SRR4302532.4834791/1‑94 (MQ=60)
        CAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTT                                                                                    <  SRR4302532.3092010/94‑1 (MQ=60)
        CAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTT                                                                                    >  SRR4302532.2432878/1‑94 (MQ=60)
         AGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTT                                                                                   <  SRR4302532.5317824/94‑1 (MQ=60)
           ACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTTTTGCGCGACTTTCT                                                                                 <  SRR4302532.5010204/94‑1 (MQ=60)
           ACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCT                                                                                 <  SRR4302532.1724325/94‑1 (MQ=60)
           ACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCT                                                                                 >  SRR4302532.2021145/1‑94 (MQ=60)
           ACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCT                                                                                 <  SRR4302532.3186871/94‑1 (MQ=60)
            CTGGTAAATATTTACTGAGGGGGAAATTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTG                                                                                <  SRR4302532.2748407/94‑1 (MQ=60)
            CTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTGCGCGACTTTCTG                                                                                >  SRR4302532.4625080/1‑94 (MQ=60)
             TGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGT                                                                               <  SRR4302532.216753/94‑1 (MQ=60)
             TGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCCGTGGTTAACCGGGAGGAAACGCTGTTGCGCGACTTTTTGT                                                                               >  SRR4302532.4610346/1‑94 (MQ=60)
               GTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTTCCGGTGGTTAACCGGGAGGAAACGCTGTATAGCGACTTTCTGTTT                                                                             >  SRR4302532.1917233/1‑94 (MQ=60)
                 AAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAC                                                                           <  SRR4302532.2521207/94‑1 (MQ=60)
                  AATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACC                                                                          >  SRR4302532.4194671/1‑94 (MQ=60)
                   ATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTGACCGGGAGTAAAAGCCGTTTCGCGACTTTTTGTTTAACG                                                                         >  SRR4302532.3632414/1‑94 (MQ=60)
                   ATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCG                                                                         >  SRR4302532.1191159/1‑94 (MQ=60)
                   ATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCG                                                                         <  SRR4302532.5450096/94‑1 (MQ=60)
                       TTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCAA                                                                     >  SRR4302532.5173989/1‑94 (MQ=60)
                         ACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATC                                                                   >  SRR4302532.3532682/1‑94 (MQ=60)
                         ACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATC                                                                   >  SRR4302532.3993781/1‑94 (MQ=60)
                         ACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTACGCGCCTTTCTGTTTACCGGCAATC                                                                   >  SRR4302532.5566505/1‑94 (MQ=60)
                          CTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATtc                                                                  >  SRR4302532.3489720/1‑92 (MQ=60)
                          CTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCA                                                                  <  SRR4302532.711389/94‑1 (MQ=60)
                          CTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCA                                                                  <  SRR4302532.6350066/94‑1 (MQ=60)
                          CTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGCTGCGCGACTTTCTGTTTACAGGCCATCA                                                                  >  SRR4302532.516065/1‑94 (MQ=60)
                          CTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGGTAACCGGGAGGAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCA                                                                  >  SRR4302532.5952958/1‑94 (MQ=60)
                           TGAGGGTGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCAC                                                                 <  SRR4302532.1491275/94‑1 (MQ=60)
                           TGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCAC                                                                 <  SRR4302532.5926716/94‑1 (MQ=60)
                              GGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACGTTCTGTTTACCGGCAATCACTCC                                                              >  SRR4302532.3668768/1‑94 (MQ=60)
                               GGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCA                                                             <  SRR4302532.3036866/94‑1 (MQ=60)
                                  GAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATA                                                          <  SRR4302532.4173200/94‑1 (MQ=60)
                                  GAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATA                                                          <  SRR4302532.4518902/94‑1 (MQ=60)
                                  GAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATA                                                          >  SRR4302532.1943198/1‑94 (MQ=60)
                                  GAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATA                                                          <  SRR4302532.2450325/94‑1 (MQ=60)
                                   AAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGCCTTTCTGTTTACCGGCAATCACTCCAATAA                                                         >  SRR4302532.308788/1‑94 (MQ=60)
                                   AAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCAATCAATCCAATAA                                                         >  SRR4302532.1517483/1‑94 (MQ=60)
                                   AAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAA                                                         <  SRR4302532.3719968/94‑1 (MQ=60)
                                   AAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAA                                                         >  SRR4302532.575866/1‑94 (MQ=60)
                                    AAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAA                                                        <  SRR4302532.5069533/94‑1 (MQ=60)
                                     AGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAAC                                                       <  SRR4302532.2649895/94‑1 (MQ=60)
                                     AGTTTCTATTGAGTCAGTGGAAGGCTACCGGTGGTTAACCGGGAGAAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAAC                                                       >  SRR4302532.4320366/1‑94 (MQ=60)
                                      GTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACG                                                      <  SRR4302532.1565/94‑1 (MQ=60)
                                      GTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACG                                                      <  SRR4302532.5086192/94‑1 (MQ=60)
                                      GTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTGCGCGACTTTCTGTTTACCCGCAATCACTCCAATAAACG                                                      >  SRR4302532.2624136/1‑94 (MQ=60)
                                        TTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCCATAAACGCC                                                    >  SRR4302532.6042735/1‑94 (MQ=60)
                                        TTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCC                                                    <  SRR4302532.1283850/94‑1 (MQ=60)
                                         TCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCT                                                   >  SRR4302532.2249467/1‑94 (MQ=60)
                                          CTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTG                                                  >  SRR4302532.1943600/1‑94 (MQ=60)
                                          CTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATAACTCCAATAAACGCCTG                                                  >  SRR4302532.1943601/1‑94 (MQ=60)
                                           TATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGC                                                 >  SRR4302532.3062449/1‑94 (MQ=60)
                                           TATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGC                                                 >  SRR4302532.4386566/1‑94 (MQ=60)
                                           TATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGC                                                 <  SRR4302532.2215091/94‑1 (MQ=60)
                                           TATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGC                                                 <  SRR4302532.2229887/94‑1 (MQ=60)
                                           TATTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGC                                                 >  SRR4302532.1567782/1‑94 (MQ=60)
                                            ATTGGGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCCGTTGCGCGACTTTCTGTTTACCCGCAATCACTCCAATACACGCCTGCC                                                >  SRR4302532.3653534/1‑94 (MQ=60)
                                            ATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCA                                                >  SRR4302532.974708/1‑94 (MQ=60)
                                              TGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTcccac                                              >  SRR4302532.4602140/1‑89 (MQ=60)
                                              TGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACC                                              <  SRR4302532.5121263/94‑1 (MQ=60)
                                              TGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACC                                              <  SRR4302532.4902799/94‑1 (MQ=60)
                                              TGAGTCAGTGGAAGGCCCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACC                                              <  SRR4302532.5121262/94‑1 (MQ=60)
                                               GAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCT                                             <  SRR4302532.1478290/94‑1 (MQ=60)
                                                 GTCAGTGGAAGGGTCCCGGTGGTTAACCGGGAGTAAACGCTTTTTTGAGAGTTTCTGTTTACCGGCAATCACcccactaaaagcctctacccgc                                           >  SRR4302532.5254005/1‑72 (MQ=48)
                                                 GTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAAAGCCTGCACCTGC                                           >  SRR4302532.5574029/1‑94 (MQ=60)
                                                 GTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGAAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCCGC                                           >  SRR4302532.5367049/1‑94 (MQ=60)
                                                  TCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCT                                          <  SRR4302532.1816601/94‑1 (MQ=60)
                                                  TCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCT                                          <  SRR4302532.3805452/94‑1 (MQ=60)
                                                    AGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTGCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTT                                        <  SRR4302532.1670682/94‑1 (MQ=60)
                                                    AGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAAAGCCTGCACCTGCCTT                                        >  SRR4302532.4036765/1‑94 (MQ=60)
                                                     GTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTT                                       >  SRR4302532.4001405/1‑94 (MQ=60)
                                                       GGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGT                                     <  SRR4302532.1413605/94‑1 (MQ=60)
                                                         AAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTT                                   <  SRR4302532.4520140/94‑1 (MQ=60)
                                                          AGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCAGTCCAATAAACGCCTGCACCTGCTTTTGTTTA                                  <  SRR4302532.5584145/94‑1 (MQ=60)
                                                          AGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTA                                  >  SRR4302532.3740141/1‑94 (MQ=60)
                                                            GTTCCTGGTGGTTAACCGGGATTAATCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG                                <  SRR4302532.6134633/94‑1 (MQ=60)
                                                            GCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG                                >  SRR4302532.473229/1‑94 (MQ=60)
                                                             CTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGC                               >  SRR4302532.1804889/1‑94 (MQ=60)
                                                              TCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCG                              <  SRR4302532.5500326/94‑1 (MQ=60)
                                                               CCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGC                             <  SRR4302532.4035789/94‑1 (MQ=60)
                                                                CCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCC                            >  SRR4302532.5019441/1‑94 (MQ=60)
                                                                 CGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCG                           >  SRR4302532.4165825/1‑94 (MQ=60)
                                                                   GTGGTTAAGCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGAC                         <  SRR4302532.2885331/94‑1 (MQ=60)
                                                                    TGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACA                        >  SRR4302532.1895475/1‑94 (MQ=60)
                                                                     GGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAG                       <  SRR4302532.1926056/94‑1 (MQ=60)
                                                                     GGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAG                       <  SRR4302532.153997/94‑1 (MQ=60)
                                                                      GTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGT                      >  SRR4302532.5738865/1‑94 (MQ=60)
                                                                       TTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTT                     <  SRR4302532.3724749/94‑1 (MQ=60)
                                                                       TTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTT                     <  SRR4302532.4659567/94‑1 (MQ=60)
                                                                        TAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTT                    <  SRR4302532.78417/94‑1 (MQ=60)
                                                                        TAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTT                    >  SRR4302532.3347389/1‑94 (MQ=60)
                                                                           CCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCA                 <  SRR4302532.5231566/94‑1 (MQ=60)
                                                                            CGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCAC                <  SRR4302532.3185547/94‑1 (MQ=60)
                                                                            CGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCAC                <  SRR4302532.5497828/94‑1 (MQ=60)
                                                                             GGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACA               >  SRR4302532.2058791/1‑94 (MQ=60)
                                                                             GGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACA               <  SRR4302532.6238940/94‑1 (MQ=60)
                                                                              GGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACAC              <  SRR4302532.495970/94‑1 (MQ=60)
                                                                              GGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACAC              <  SRR4302532.5035055/94‑1 (MQ=60)
                                                                              GGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACAC              <  SRR4302532.495971/94‑1 (MQ=60)
                                                                                AGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCT            <  SRR4302532.5028329/94‑1 (MQ=60)
                                                                                AGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCT            <  SRR4302532.5028330/94‑1 (MQ=60)
                                                                                 GTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTG           <  SRR4302532.5088217/94‑1 (MQ=60)
                                                                                 ATAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTG           <  SRR4302532.191651/94‑1 (MQ=60)
                                                                                  TAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGG          <  SRR4302532.6165651/94‑1 (MQ=60)
                                                                                   AAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGC         <  SRR4302532.2114726/94‑1 (MQ=60)
                                                                                   AAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGC         >  SRR4302532.823786/1‑94 (MQ=60)
                                                                                    AACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCG        >  SRR4302532.3965244/1‑94 (MQ=60)
                                                                                    AACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCG        >  SRR4302532.5166197/1‑94 (MQ=60)
                                                                                      CGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTA      >  SRR4302532.3005941/1‑94 (MQ=60)
                                                                                      CGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTA      <  SRR4302532.1654147/94‑1 (MQ=60)
                                                                                      CGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTA      <  SRR4302532.3703995/94‑1 (MQ=60)
                                                                                       GCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAG     <  SRR4302532.6089303/94‑1 (MQ=60)
                                                                                       GCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAG     <  SRR4302532.193743/94‑1 (MQ=60)
                                                                                       GCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAG     >  SRR4302532.1919696/1‑94 (MQ=60)
                                                                                       GCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAG     <  SRR4302532.2228524/94‑1 (MQ=60)
                                                                                         TGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG   >  SRR4302532.3571582/1‑94 (MQ=60)
                                                                                          GTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGA  >  SRR4302532.711427/1‑94 (MQ=60)
                                                                                          GTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGA  >  SRR4302532.4666881/1‑94 (MQ=60)
                                                                                             |                                                                                          
TGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGA  >  NC_000913/4576978‑4577161

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