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 593 264.2 13311096 97.8% 13018251 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,117,945 G→A P451S (CCG→TCG)  yghK ← putative permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,117,9450GA86.2% 562.1 / 59.5 195P451S (CCG→TCG) yghKputative permease
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base A (115/53);  total (127/68)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 222 total aligned reads displayed.

GGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGC  >  NC_000913/3117853‑3118034
                                                                                            |                                                                                         
gtgcGATTGCAGTGAACCAAACGGGGCATAAGGGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACgg                                                                                           <  1:2384772/90‑1 (MQ=25)
ggTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACgg                                                                                           >  2:2473156/1‑93 (MQ=255)
   ggtttGCAGTAACCAAAACGGGGCGTAAGAGGAGGTGCCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGaa                                                                                        <  2:5777063/90‑1 (MQ=21)
   cGATTGCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGTTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGaa                                                                                        <  1:5343256/93‑1 (MQ=25)
    gATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAAc                                                                                       >  1:968672/1‑93 (MQ=255)
     aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACa                                                                                      >  2:2969266/1‑93 (MQ=255)
     aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACa                                                                                      >  1:1902201/1‑93 (MQ=255)
       tGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATc                                                                                    >  1:1900505/1‑93 (MQ=255)
        gCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCa                                                                                   >  1:603504/1‑93 (MQ=255)
        gCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCa                                                                                   <  1:203284/93‑1 (MQ=255)
         cAGTGACCAAAACGGGCCGTTAGAGGAGGTGCCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCAc                                                                                  <  1:6403471/93‑1 (MQ=25)
            tGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGcc                                                                               >  1:4506266/1‑93 (MQ=255)
            gaacccAAACGGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGcc                                                                               <  2:3814939/89‑1 (MQ=255)
              aaCCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTg                                                                             >  2:3443462/1‑93 (MQ=255)
                ccAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTa                                                                           >  2:124541/1‑93 (MQ=255)
                ccAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTa                                                                           >  2:4132774/1‑93 (MQ=255)
                 cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                                                          >  2:6301940/1‑93 (MQ=255)
                 cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                                                          >  1:4662038/1‑93 (MQ=255)
                 cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                                                          >  2:4318449/1‑93 (MQ=255)
                  aaaaAGGGTGTTAGGGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAcc                                                                         <  1:1958748/93‑1 (MQ=25)
                   aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCt                                                                        >  2:6010800/1‑93 (MQ=255)
                    acgggggGTAAGAGGAGGGGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTg                                                                       <  2:2557488/90‑1 (MQ=255)
                    aCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTg                                                                       >  1:4996463/1‑93 (MQ=255)
                       gggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCa                                                                    >  2:889654/1‑93 (MQ=255)
                       gggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCa                                                                    >  1:5443636/1‑93 (MQ=255)
                       gggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACAGCACGCCTGTACCTGCCa                                                                    >  1:3720094/1‑93 (MQ=255)
                        ggCGTTAGCGGAGGTGTCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAg                                                                   <  2:854557/93‑1 (MQ=37)
                        ggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAg                                                                   >  1:3105561/1‑93 (MQ=255)
                           gTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAc                                                                >  1:109029/1‑93 (MQ=255)
                           gTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAc                                                                >  2:3614032/1‑93 (MQ=255)
                             tgggggaggTTCCCGAGCGGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAcca                                                              <  2:1459488/89‑1 (MQ=255)
                              agagGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccag                                                             >  2:4980090/1‑93 (MQ=255)
                               gaggagGTGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagc                                                            <  2:5394554/93‑1 (MQ=255)
                                aggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATGACGCCTGTACCTGCCAGTAccagcg                                                           >  1:5513416/1‑93 (MQ=255)
                                aggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcg                                                           >  1:4464902/1‑93 (MQ=255)
                                 ggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgc                                                          >  1:2752082/1‑93 (MQ=255)
                                  gaggTACCCGACCCGTTAAGAAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAAAACATCACGCCTGTACCTGCCAGTAccagcgcc                                                         <  1:3803865/93‑1 (MQ=255)
                                    ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                                                       >  2:5483829/1‑93 (MQ=255)
                                    ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                                                       >  1:3610822/1‑93 (MQ=255)
                                    ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                                                       >  1:4335208/1‑93 (MQ=255)
                                     gtgtCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccagc                                                      >  1:2960456/1‑93 (MQ=255)
                                        tCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtg                                                   >  2:1511734/1‑93 (MQ=255)
                                        gcggAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtg                                                   <  2:4475702/90‑1 (MQ=255)
                                          cGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt                                                 <  1:6149973/93‑1 (MQ=255)
                                          cGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt                                                 <  1:4099929/93‑1 (MQ=255)
                                            aGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGtggtgg                                               >  1:659453/1‑93 (MQ=255)
                                              cgggTAAGGAATGCGCACGGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAc                                             <  1:743423/91‑1 (MQ=255)
                                              ccGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAc                                             >  2:4654746/1‑93 (MQ=255)
                                               cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGACCGAGAACATCACGCTTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACa                                            <  2:5914295/93‑1 (MQ=255)
                                               cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCGGCGCCAGCGTGGTGGACa                                            >  1:5908432/1‑93 (MQ=255)
                                               cGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACa                                            <  2:6466579/93‑1 (MQ=255)
                                                 gTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAcgtc                                          >  2:512525/1‑90 (MQ=255)
                                                    aGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCa                                       <  2:2680041/93‑1 (MQ=255)
                                                      gAATCCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGa                                     <  1:3205736/93‑1 (MQ=255)
                                                      gAATCCCCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCCCGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGa                                     <  1:1684850/93‑1 (MQ=38)
                                                       aaggcGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGaa                                    <  2:1130456/89‑1 (MQ=255)
                                                       aaTGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGaa                                    >  2:2883552/1‑93 (MQ=255)
                                                         gccGCCGAGCAAGCCGAGAATCGCTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGTACATGCCAGAATa                                  <  2:2171877/91‑1 (MQ=255)
                                                          gcGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAg                                 <  2:3204526/92‑1 (MQ=255)
                                                             cccAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTg                              >  1:4215726/1‑93 (MQ=255)
                                                                agccagccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTg                           >  2:4178306/1‑93 (MQ=255)
                                                                 gccagccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGa                          >  1:5392979/1‑93 (MQ=255)
                                                                   cagccGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACg                        >  2:369176/1‑93 (MQ=255)
                                                                    agccGAGAAACGGTGAGAAGAACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGa                       <  1:1606470/93‑1 (MQ=255)
                                                                    agccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGa                       >  2:6169871/1‑93 (MQ=255)
                                                                     gcggagaAAGGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      <  2:1195833/90‑1 (MQ=255)
                                                                     gccGAGAAAGGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACgca                      <  1:3465831/93‑3 (MQ=255)
                                                                     gccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      >  1:2085463/1‑93 (MQ=255)
                                                                     gccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      <  2:1860469/93‑1 (MQ=255)
                                                                     gccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      <  2:696140/93‑1 (MQ=255)
                                                                      ccGAGAAAGTGTGAGAAGAACGAGAACCTCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAg                     <  2:5592849/93‑1 (MQ=255)
                                                                      ccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAg                     <  1:6367213/93‑1 (MQ=255)
                                                                      ccGAGAAACGCTGAGAAGAGCGAGACCTTCCCGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAg                     <  1:5455319/93‑1 (MQ=255)
                                                                       cGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAgg                    <  2:5477515/93‑1 (MQ=255)
                                                                       cGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGGGACGAAgg                    >  2:656329/1‑93 (MQ=255)
                                                                        gagaAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGc                   >  1:413429/1‑93 (MQ=255)
                                                                        gagaAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGc                   <  2:4915628/93‑1 (MQ=255)
                                                                         agaAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCg                  <  1:1141775/93‑1 (MQ=255)
                                                                         agaAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCg                  >  1:4001152/1‑93 (MQ=255)
                                                                          gaAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGa                 >  2:4550675/1‑93 (MQ=255)
                                                                           aaaCGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGaa                <  1:4292591/93‑1 (MQ=255)
                                                                           aaaCGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGaa                >  1:5841928/1‑93 (MQ=255)
                                                                            aaCGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAAc               >  2:441329/1‑93 (MQ=255)
                                                                              cGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGc             <  2:5193632/93‑1 (MQ=255)
                                                                               ggTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGcc            >  2:169296/1‑93 (MQ=255)
                                                                               ggTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGcc            <  1:6357012/93‑1 (MQ=255)
                                                                                gTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGCACGCCa           >  2:1400391/1‑93 (MQ=255)
                                                                                 ggagaAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAg          <  2:2739730/92‑1 (MQ=255)
                                                                                  gagaAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGc         >  2:4378743/1‑93 (MQ=255)
                                                                                   agaagaACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCa        <  1:1158615/93‑1 (MQ=255)
                                                                                   agaagaACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCa        >  1:4969150/1‑93 (MQ=255)
                                                                                    gaagaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAc       <  2:5373781/93‑1 (MQ=255)
                                                                                    aaagaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAc       <  2:2629817/92‑1 (MQ=255)
                                                                                     aagaaCGGGAACATCCCGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAcc      <  2:4274395/93‑1 (MQ=255)
                                                                                     aagaaCGGAAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAcc      <  2:6397433/93‑1 (MQ=255)
                                                                                      agaaCGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCa     >  2:5663711/1‑93 (MQ=255)
                                                                                      agaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCa     >  2:4844919/1‑93 (MQ=255)
                                                                                      agaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAAATGCCAGAATAGTTGGTGACGAAGGCGCACGCCCGCAccc     >  1:5158560/1‑92 (MQ=255)
                                                                                         aCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGc  <  2:330031/93‑1 (MQ=255)
                                                                                         aCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGc  >  1:250222/1‑93 (MQ=255)
                                                                                            |                                                                                         
GGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGC  >  NC_000913/3117853‑3118034

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ 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
TGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAAT  >  NC_000913/3117852‑3118038
                                                                                             |                                                                                             
TGGTGGATTGCAGTGACCAAAACAGGGCGTGAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGA                                                                                               <  SRR4302547.930434/94‑1 (MQ=60)
TGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGG                                                                                               >  SRR4302547.1159076/1‑94 (MQ=60)
 gtgCGATTGCAGTGAACCAAACGGGGCATAAGGGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGG                                                                                              <  SRR4302547.2433084/91‑1 (MQ=60)
 GGTCGATTGCAGTGAACCAAACGGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGG                                                                                              <  SRR4302547.4681928/94‑1 (MQ=60)
 GGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAG                                                                                              >  SRR4302547.990614/1‑94 (MQ=60)
 GGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAG                                                                                              >  SRR4302547.3612811/1‑94 (MQ=60)
  GTGGATTGCAGAGACCAAAACAGGGCGTGAGAGGAGGGGCCGGAGCGGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGA                                                                                             <  SRR4302547.4169469/94‑1 (MQ=42)
  GTCGATTGCAGTGACCAAAACAGGGCGTTAGAGGAGGTGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGA                                                                                             <  SRR4302547.3670756/94‑1 (MQ=60)
    CGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAAC                                                                                           >  SRR4302547.987476/1‑94 (MQ=60)
    CGATTGCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGTTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAAC                                                                                           <  SRR4302547.5453413/94‑1 (MQ=60)
     GATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACA                                                                                          >  SRR4302547.1705176/1‑94 (MQ=60)
     GATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACA                                                                                          >  SRR4302547.4283700/1‑94 (MQ=60)
     GATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACA                                                                                          >  SRR4302547.1939831/1‑94 (MQ=60)
      ATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACAT                                                                                         >  SRR4302547.5300524/1‑94 (MQ=60)
      ATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACAT                                                                                         >  SRR4302547.2428527/1‑94 (MQ=60)
       ttgcaggaacccaccCAGGGCGTAAGAGGAGGTGCCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATC                                                                                        <  SRR4302547.6072014/79‑1 (MQ=60)
       TTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATC                                                                                        >  SRR4302547.1938108/1‑94 (MQ=60)
       TTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATC                                                                                        >  SRR4302547.3160204/1‑94 (MQ=60)
        TGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCA                                                                                       >  SRR4302547.614701/1‑94 (MQ=60)
         GCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCAC                                                                                      <  SRR4302547.207086/94‑1 (MQ=60)
          CAGTGACCAAAACGGGCCGTTAGAGGAGGTGCCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACG                                                                                     <  SRR4302547.6537033/94‑1 (MQ=60)
          CAGTGACCAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACG                                                                                     <  SRR4302547.5567313/94‑1 (MQ=60)
           AGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGC                                                                                    >  SRR4302547.2562352/1‑94 (MQ=60)
            GTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCC                                                                                   >  SRR4302547.4599321/1‑94 (MQ=60)
                acaaaagggggcgttggtgggggggtcgGAGCCGGTAGGGAAGGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTA                                                                               <  SRR4302547.1997232/66‑1 (MQ=60)
                ACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGCAAGGAATGCGCCCAGCCAGCCGAGACACGGTGAGAAGAACGAGCACATCACGCCTGTA                                                                               >  SRR4302547.1845931/1‑94 (MQ=60)
                 CCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAC                                                                              >  SRR4302547.4758590/1‑94 (MQ=60)
                   aagcgggggttaagagggggtggcgGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCT                                                                            <  SRR4302547.5616246/69‑1 (MQ=60)
                   aaacggggtttaagaggggttgccgGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCT                                                                            <  SRR4302547.1528975/69‑1 (MQ=60)
                   AAAAAGGGTGTTAGGGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCT                                                                            <  SRR4302547.1997633/94‑1 (MQ=60)
                    acaAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTG                                                                           <  SRR4302547.3393957/91‑1 (MQ=60)
                    AACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTG                                                                           >  SRR4302547.5099782/1‑94 (MQ=60)
                    AACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTG                                                                           >  SRR4302547.5265925/1‑94 (MQ=60)
                     ACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGC                                                                          >  SRR4302547.5970686/1‑94 (MQ=60)
                       AGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCA                                                                        >  SRR4302547.5556118/1‑94 (MQ=60)
                       AGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACAGCACGCCTGTACCTGCCA                                                                        >  SRR4302547.3796032/1‑94 (MQ=60)
                        GGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAG                                                                       >  SRR4302547.3167110/1‑94 (MQ=60)
                           CGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAC                                                                    >  SRR4302547.111017/1‑94 (MQ=60)
                              tagaggggcgcccgGAGCCGGTAAGGAACGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAG                                                                 <  SRR4302547.5144454/80‑1 (MQ=60)
                                GAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATGACGCCTGTACCTGCCAGTACCAGCG                                                               >  SRR4302547.5627236/1‑94 (MQ=60)
                                GAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCG                                                               >  SRR4302547.1470319/1‑94 (MQ=60)
                                GAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCG                                                               >  SRR4302547.4557154/1‑94 (MQ=60)
                                 ggggggtggCCGAGCGGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGC                                                              <  SRR4302547.1930185/85‑1 (MQ=60)
                                 AGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGC                                                              >  SRR4302547.2807263/1‑94 (MQ=60)
                                   gggtggccgGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCA                                                            <  SRR4302547.3050353/85‑1 (MQ=60)
                                   gaggtacccgacccgttaagaaatgCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAAAACATCACGCCTGTACCTGCCAGTACCAGCGCCA                                                            <  SRR4302547.3881408/69‑1 (MQ=60)
                                    AGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAG                                                           >  SRR4302547.4424387/1‑94 (MQ=60)
                                    AGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAG                                                           >  SRR4302547.3684475/1‑94 (MQ=60)
                                    AGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAG                                                           >  SRR4302547.6496437/1‑94 (MQ=60)
                                     GGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGC                                                          >  SRR4302547.3019379/1‑94 (MQ=60)
                                        GTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTG                                                       >  SRR4302547.1986071/1‑94 (MQ=60)
                                           CGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTG                                                    <  SRR4302547.6277567/94‑1 (MQ=60)
                                           CGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTG                                                    <  SRR4302547.4591648/94‑1 (MQ=60)
                                           CGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTG                                                    <  SRR4302547.4183310/94‑1 (MQ=60)
                                            GAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGG                                                   >  SRR4302547.671599/1‑94 (MQ=60)
                                            GAGCCGGTAAGAAATCCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGG                                                   <  SRR4302547.3652098/94‑1 (MQ=60)
                                               CGGGTAAGGAATGCGCACGGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACA                                                <  SRR4302547.757415/94‑1 (MQ=60)
                                               CCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCGGCGCCAGCGTGGTGGACA                                                >  SRR4302547.6030789/1‑94 (MQ=60)
                                                ggtaaagaaagcCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACAT                                               <  SRR4302547.4628759/82‑1 (MQ=60)
                                                gggaaggaagggcccccacCCAGCAGAGAAAGGGAGAGAAGAACGAAAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACAT                                               <  SRR4302547.46713/75‑1 (MQ=60)
                                                CGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACAT                                               <  SRR4302547.3980547/94‑1 (MQ=60)
                                                    AAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCA                                           >  SRR4302547.4314169/1‑94 (MQ=60)
                                                    AAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCA                                           >  SRR4302547.2847657/1‑94 (MQ=60)
                                                     AGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAG                                          >  SRR4302547.6377772/1‑94 (MQ=60)
                                                      gaaacccccccgcCCAGCGGAGAACCGGTGAAAAGACCGGGAACATCACGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGA                                         <  SRR4302547.3472928/81‑1 (MQ=60)
                                                      GGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGA                                         >  SRR4302547.5521009/1‑94 (MQ=60)
                                                       GAATCCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAA                                        <  SRR4302547.3269670/94‑1 (MQ=60)
                                                       GAATCCCCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCCCGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAA                                        <  SRR4302547.1718009/94‑1 (MQ=60)
                                                         ATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATA                                      >  SRR4302547.2282089/1‑94 (MQ=60)
                                                         ATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCACCGTGGTGGACATGCCAGAATA                                      <  SRR4302547.6532006/94‑1 (MQ=60)
                                                           GCGCCCAGCCAGCCGAGAAACGGTGAGAAGGACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGT                                    >  SRR4302547.932828/1‑94 (MQ=60)
                                                            CGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTT                                   >  SRR4302547.569322/1‑94 (MQ=60)
                                                             GCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTG                                  >  SRR4302547.4301783/1‑94 (MQ=60)
                                                             GCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTG                                  >  SRR4302547.5585630/1‑94 (MQ=60)
                                                              CCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGG                                 >  SRR4302547.2435912/1‑94 (MQ=60)
                                                              CCCAGCCAGCCGAGAAACGGGGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGCGGACATGCCAGAATAGTGGG                                 >  SRR4302547.2058154/1‑94 (MQ=60)
                                                                CAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTG                               >  SRR4302547.5501492/1‑94 (MQ=60)
                                                                CAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTG                               >  SRR4302547.960537/1‑94 (MQ=60)
                                                                CAGCCAGCCGAGAAACGGTGAGAAGAAATAGAACAACACCCGTACACCTGCCAGTACCAGCGGCAGCGTGGTGGACATGCCAGAATAGTTGGTG                               >  SRR4302547.827634/1‑94 (MQ=53)
                                                                 AGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGA                              >  SRR4302547.5504384/1‑94 (MQ=60)
                                                                 AGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGA                              >  SRR4302547.4546338/1‑94 (MQ=60)
                                                                   CCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACG                            <  SRR4302547.5685219/94‑1 (MQ=60)
                                                                   CCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACG                            >  SRR4302547.6374777/1‑94 (MQ=60)
                                                                     AGCCGAGAAACGGTGAGAAGAACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAA                          <  SRR4302547.1638315/94‑1 (MQ=60)
                                                                     AGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAA                          >  SRR4302547.3770802/1‑94 (MQ=60)
                                                                     AGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAA                          >  SRR4302547.2127051/1‑94 (MQ=60)
                                                                     AGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAA                          >  SRR4302547.4247028/1‑94 (MQ=60)
                                                                      GCCGAGAAAGGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGCAG                         <  SRR4302547.3536513/94‑1 (MQ=60)
                                                                      GCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAG                         >  SRR4302547.1938898/1‑94 (MQ=60)
                                                                       CCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGG                        <  SRR4302547.6500049/94‑1 (MQ=60)
                                                                       CCGAGAAACGCTGAGAAGAGCGAGACCTTCCCGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGG                        <  SRR4302547.5568038/94‑1 (MQ=60)
                                                                        CGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGC                       >  SRR4302547.6778019/1‑94 (MQ=60)
                                                                        CGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGC                       <  SRR4302547.4188920/94‑1 (MQ=60)
                                                                        CGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGC                       >  SRR4302547.421101/1‑94 (MQ=60)
                                                                        CGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGC                       >  SRR4302547.439318/1‑94 (MQ=60)
                                                                        CGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGC                       <  SRR4302547.6665959/94‑1 (MQ=60)
                                                                         GAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCG                      >  SRR4302547.4082482/1‑94 (MQ=60)
                                                                          AGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGA                     <  SRR4302547.1164839/94‑1 (MQ=60)
                                                                          AGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGA                     >  SRR4302547.3192500/1‑94 (MQ=60)
                                                                           GAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAA                    >  SRR4302547.5962790/1‑94 (MQ=60)
                                                                            AAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAAC                   <  SRR4302547.4380823/94‑1 (MQ=60)
                                                                             AACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACG                  >  SRR4302547.6215154/1‑94 (MQ=60)
                                                                             AACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACG                  >  SRR4302547.3847030/1‑94 (MQ=60)
                                                                             AACGGTGAGAAGAACGAAAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACG                  <  SRR4302547.3818722/94‑1 (MQ=60)
                                                                                GGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCA               >  SRR4302547.4111825/1‑94 (MQ=60)
                                                                                GGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCA               <  SRR4302547.6489560/94‑1 (MQ=60)
                                                                                  TGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGC             >  SRR4302547.6130263/1‑94 (MQ=60)
                                                                                  TGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGC             >  SRR4302547.1178861/1‑94 (MQ=60)
                                                                                   GAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCCGCA            >  SRR4302547.2114232/1‑94 (MQ=60)
                                                                                   GAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCA            >  SRR4302547.5071974/1‑94 (MQ=60)
                                                                                   GAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCA            <  SRR4302547.3097084/94‑1 (MQ=60)
                                                                                    AGAAGAACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAC           <  SRR4302547.1182051/94‑1 (MQ=60)
                                                                                    AGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAC           <  SRR4302547.4738154/94‑1 (MQ=60)
                                                                                     GAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACC          >  SRR4302547.509881/1‑94 (MQ=60)
                                                                                      AAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAAATGCCAGAATAGTTGGTGACGAAGGCGCACGCCCGCACCC         >  SRR4302547.5264782/1‑94 (MQ=60)
                                                                                        AAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATG       <  SRR4302547.5067211/94‑1 (MQ=60)
                                                                                         AACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGC      >  SRR4302547.254906/1‑94 (MQ=60)
                                                                                         AACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGC      >  SRR4302547.5093381/1‑94 (MQ=60)
                                                                                           CGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCA    >  SRR4302547.1517175/1‑94 (MQ=60)
                                                                                             AGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAAT  >  SRR4302547.6254653/1‑94 (MQ=60)
                                                                                             |                                                                                             
TGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAAT  >  NC_000913/3117852‑3118038

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