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 3,333,863 A→G R91R (CGT→CGC)  murA ← UDP‑N‑glucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,8630AG99.0% 713.9 / ‑2.8 207R91R (CGT→CGC) murAUDP‑N‑glucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (97/108);  total (99/108)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 214 total aligned reads displayed.

GGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCG  >  NC_000913/3333774‑3333951
                                                                                         |                                                                                        
 gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCa                                                                                      >  1:963989/1‑93 (MQ=255)
 gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCa                                                                                      <  1:1001378/93‑1 (MQ=255)
 gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCa                                                                                      <  2:2373751/93‑1 (MQ=255)
    gcgcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATgg                                                                                   <  1:5438117/93‑1 (MQ=255)
     cgcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGt                                                                                  <  1:3715860/93‑1 (MQ=255)
      gcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGtt                                                                                 <  2:273736/93‑1 (MQ=255)
          cGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTaa                                                                             >  2:6026136/1‑93 (MQ=255)
           gATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAc                                                                            >  2:2838565/1‑93 (MQ=255)
            aTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAcc                                                                           >  1:103138/1‑93 (MQ=255)
            aTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAcc                                                                           >  2:1446763/1‑93 (MQ=255)
            aTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAcc                                                                           >  1:58156/1‑93 (MQ=255)
             tCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCa                                                                          <  1:3887034/93‑1 (MQ=255)
                tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAt                                                                       >  2:5817912/1‑93 (MQ=255)
                tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAt                                                                       >  1:3849769/1‑93 (MQ=255)
                tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAt                                                                       >  2:2437608/1‑93 (MQ=255)
                 aCAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATc                                                                      <  2:1036004/93‑1 (MQ=255)
                   aaCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCgc                                                                    >  2:3787695/1‑92 (MQ=255)
                   aaCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGt                                                                    >  2:503048/1‑93 (MQ=255)
                    aCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTa                                                                   <  1:1866733/93‑1 (MQ=255)
                     ccGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTaa                                                                  <  2:3497869/93‑1 (MQ=255)
                       gCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAgg                                                                >  2:2105401/1‑93 (MQ=255)
                        ccAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGt                                                               <  2:5756754/93‑1 (MQ=255)
                        ccAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGt                                                               <  1:2474320/93‑1 (MQ=255)
                         cAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTg                                                              >  2:4849654/1‑93 (MQ=255)
                         cAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTg                                                              >  2:3008396/1‑93 (MQ=255)
                          aGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTgc                                                             >  2:3607118/1‑93 (MQ=255)
                            gTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTgcgc                                                           >  2:6052039/1‑93 (MQ=255)
                            gTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTgcgc                                                           >  1:2802113/1‑93 (MQ=255)
                                tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                                       >  2:22876/1‑93 (MQ=255)
                                tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                                       >  2:1994070/1‑93 (MQ=255)
                                tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                                       <  1:6147859/93‑1 (MQ=255)
                                 gAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAAt                                                      <  2:829880/93‑1 (MQ=255)
                                  aaaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATa                                                     <  1:1233953/93‑1 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGGGCATGGTTTTAACCAGATCGTAAGGTGC‑CAGAATACa                                                   <  1:5016049/93‑1 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                                                    <  1:1289283/93‑1 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                                                    >  1:3257492/1‑93 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATAc                                                    >  2:4223538/1‑93 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                                                    >  1:4861024/1‑93 (MQ=255)
                                     cTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACAt                                                  <  2:4564033/93‑1 (MQ=255)
                                      ttGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAtt                                                 >  2:2501931/1‑93 (MQ=255)
                                      ttGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAtt                                                 >  2:69402/1‑93 (MQ=255)
                                      ttGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACAtt                                                 >  1:822672/1‑93 (MQ=255)
                                       tGCCCCTGACCAAAGCGCGCGACCAGCGGCCCCAGCGCCAAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTa                                                <  2:936742/93‑1 (MQ=255)
                                         ccccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAAc                                              <  2:4790932/93‑1 (MQ=255)
                                          cccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACg                                             >  1:1323963/1‑93 (MQ=255)
                                           ccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACgc                                            >  2:2557483/1‑92 (MQ=255)
                                            cTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTc                                           <  2:4016695/93‑1 (MQ=255)
                                             tGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcg                                          >  2:4943109/1‑93 (MQ=255)
                                              gACCAAAGCGCGCTGCCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgc                                         <  1:79037/93‑1 (MQ=255)
                                              gACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgc                                         <  2:3259081/93‑1 (MQ=255)
                                               aCCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgcg                                        <  1:4463269/93‑1 (MQ=255)
                                                ccAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCgg                                       >  1:5390394/1‑93 (MQ=255)
                                                 cAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCggg                                      <  2:1115118/93‑1 (MQ=255)
                                                  aaaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGGAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGc                                     >  1:6083474/1‑93 (MQ=255)
                                                  aaaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGc                                     >  1:6104097/1‑93 (MQ=255)
                                                   aaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCa                                    <  2:410006/93‑1 (MQ=255)
                                                   aaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCa                                    >  2:5348827/1‑93 (MQ=255)
                                                    aGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAt                                   <  2:889067/93‑1 (MQ=255)
                                                     gcgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc                                  <  1:3991161/93‑1 (MQ=255)
                                                     gcgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc                                  >  1:4072685/1‑93 (MQ=255)
                                                     gcgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc                                  <  2:3060461/93‑1 (MQ=255)
                                                       gcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCaa                                >  2:5700866/1‑93 (MQ=255)
                                                        cgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAat                               <  1:4226635/93‑1 (MQ=255)
                                                        cgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAat                               <  2:5656086/93‑1 (MQ=255)
                                                        cgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAat                               <  2:5214811/93‑1 (MQ=255)
                                                         gcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAata                              >  2:5857858/1‑93 (MQ=255)
                                                         gcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAata                              <  2:5932285/93‑1 (MQ=255)
                                                          cTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAatat                             >  2:1404557/1‑93 (MQ=255)
                                                            aCCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGc                           <  1:517543/93‑1 (MQ=255)
                                                             ccAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGca                          >  2:3446349/1‑93 (MQ=255)
                                                              cAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcac                         <  1:546763/93‑1 (MQ=255)
                                                              cAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcac                         <  2:2148568/93‑1 (MQ=255)
                                                               aGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcaca                        <  2:4474385/93‑1 (MQ=255)
                                                               aGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcaca                        >  2:1778828/1‑93 (MQ=255)
                                                                gCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAg                       >  1:2863651/1‑93 (MQ=255)
                                                                gCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAg                       >  2:1693633/1‑93 (MQ=255)
                                                                  ggCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGaa                     <  1:3275499/93‑1 (MQ=255)
                                                                   gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAc                    >  2:1984592/1‑93 (MQ=255)
                                                                   gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAc                    >  1:607413/1‑93 (MQ=255)
                                                                   gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAc                    >  1:5474168/1‑93 (MQ=255)
                                                                        aGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTa               <  1:3758420/93‑1 (MQ=255)
                                                                        aGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTa               <  2:2610200/93‑1 (MQ=255)
                                                                         gcgcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAc              >  1:2637526/1‑93 (MQ=255)
                                                                           gcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGt            <  2:1679811/93‑1 (MQ=255)
                                                                           gcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGt            <  1:3312072/93‑1 (MQ=255)
                                                                           gcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGt            >  1:4721841/1‑93 (MQ=255)
                                                                            cccAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGtt           >  1:1520656/1‑93 (MQ=255)
                                                                            cccAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACGGAACCATTACGtt           >  1:1823768/1‑93 (MQ=255)
                                                                              cAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCt         >  1:4150114/1‑93 (MQ=255)
                                                                              cAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCt         >  1:2595568/1‑93 (MQ=255)
                                                                               aGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTa        <  2:3622230/93‑1 (MQ=255)
                                                                                gATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTAc       <  2:3273865/93‑1 (MQ=255)
                                                                                gATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACAATTACGTTCTAc       >  2:1279163/1‑93 (MQ=255)
                                                                                 aTAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACt      <  2:2571072/93‑1 (MQ=255)
                                                                                  tAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     <  1:3029188/93‑1 (MQ=255)
                                                                                  tAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     >  2:4255927/1‑93 (MQ=255)
                                                                                  tAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     <  1:914719/93‑1 (MQ=255)
                                                                                  aaGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     <  1:2165060/92‑1 (MQ=255)
                                                                                   aGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACttt    <  1:4970690/93‑1 (MQ=255)
                                                                                     aaGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCg  <  2:507639/93‑1 (MQ=255)
                                                                                         |                                                                                        
GGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCG  >  NC_000913/3333774‑3333951

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ 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
AACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACC  >  NC_000913/3333770‑3333955
                                                                                             |                                                                                            
AACCGGACGCGGACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCG                                                                                              <  SRR4302532.1294558/94‑1 (MQ=60)
AACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCG                                                                                              >  SRR4302532.5874843/1‑94 (MQ=60)
  CCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCG                                                                                            <  SRR4302532.5798909/94‑1 (MQ=60)
    GGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCA                                                                                          >  SRR4302532.996882/1‑94 (MQ=60)
     GACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCAT                                                                                         <  SRR4302532.1035464/94‑1 (MQ=60)
       CGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGG                                                                                       <  SRR4302532.338777/94‑1 (MQ=60)
       CGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGG                                                                                       >  SRR4302532.1218343/1‑94 (MQ=60)
        GCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGT                                                                                      <  SRR4302532.5626581/94‑1 (MQ=60)
         CGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTT                                                                                     <  SRR4302532.3840583/94‑1 (MQ=60)
          GCACCGTTCGTACAACCGCCAGGTCGTGCACCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTT                                                                                    <  SRR4302532.5657733/94‑1 (MQ=60)
              CGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTATC                                                                                >  SRR4302532.6303808/1‑94 (MQ=60)
               GATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACC                                                                               >  SRR4302532.1514445/1‑94 (MQ=60)
               GATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACC                                                                               <  SRR4302532.1910110/94‑1 (MQ=60)
               GATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACC                                                                               >  SRR4302532.106708/1‑94 (MQ=60)
               GATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACC                                                                               >  SRR4302532.60263/1‑94 (MQ=60)
                 TCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCCTGGGTTTAACCAG                                                                             >  SRR4302532.2796065/1‑94 (MQ=60)
                 TCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAG                                                                             <  SRR4302532.4018192/94‑1 (MQ=60)
                  CGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGA                                                                            >  SRR4302532.6068265/1‑94 (MQ=60)
                   GTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAT                                                                           >  SRR4302532.3979763/1‑94 (MQ=60)
                   GTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGAGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAT                                                                           >  SRR4302532.2900855/1‑94 (MQ=60)
                     ACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCG                                                                         <  SRR4302532.1947980/94‑1 (MQ=60)
                      CAACCGCCAGGAAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGAGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGT                                                                        >  SRR4302532.437014/1‑94 (MQ=60)
                       AACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTA                                                                       >  SRR4302532.584221/1‑94 (MQ=60)
                       AACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTA                                                                       >  SRR4302532.6285752/1‑94 (MQ=60)
                       AACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTA                                                                       <  SRR4302532.5387421/94‑1 (MQ=60)
                       AACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTA                                                                       <  SRR4302532.1788072/94‑1 (MQ=60)
                        ACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAA                                                                      >  SRR4302532.4917174/1‑94 (MQ=60)
                        ACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAA                                                                      <  SRR4302532.1928576/94‑1 (MQ=60)
                         CCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAG                                                                     <  SRR4302532.4008586/94‑1 (MQ=60)
                           GCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGG                                                                   >  SRR4302532.3082639/1‑94 (MQ=60)
                            CCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTG                                                                  <  SRR4302532.2556723/94‑1 (MQ=60)
                              AGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCG                                                                >  SRR4302532.82326/1‑94 (MQ=60)
                               GGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGC                                                               >  SRR4302532.2895996/1‑94 (MQ=60)
                                GTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTGAGGTGCGCA                                                              >  SRR4302532.5313631/1‑94 (MQ=60)
                                GTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCA                                                              >  SRR4302532.6359295/1‑94 (MQ=60)
                                 TAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAG                                                             >  SRR4302532.2243204/1‑94 (MQ=60)
                                  AGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGA                                                            >  SRR4302532.5631924/1‑94 (MQ=60)
                                    TGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAAT                                                          <  SRR4302532.6359890/94‑1 (MQ=60)
                                      AAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAC                                                        >  SRR4302532.3365726/1‑94 (MQ=60)
                                      AAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAC                                                        <  SRR4302532.1275234/94‑1 (MQ=60)
                                      AAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAC                                                        >  SRR4302532.5028255/1‑94 (MQ=60)
                                       AACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGGGCATGGTTTTAACCAGATCGTAAGGTGC‑CAGAATACAT                                                      <  SRR4302532.5188017/94‑1 (MQ=60)
                                       AACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACA                                                       <  SRR4302532.1332421/94‑1 (MQ=60)
                                        ACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAT                                                      <  SRR4302532.425828/94‑1 (MQ=60)
                                         CTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATT                                                     <  SRR4302532.3344380/94‑1 (MQ=60)
                                         CTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATT                                                     >  SRR4302532.850755/1‑94 (MQ=60)
                                           TGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAA                                                   <  SRR4302532.888316/94‑1 (MQ=60)
                                             CCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACG                                                 >  SRR4302532.1368201/1‑94 (MQ=60)
                                                CTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCG                                              >  SRR4302532.1938988/1‑94 (MQ=60)
                                                CTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCG                                              <  SRR4302532.4135832/94‑1 (MQ=60)
                                                  GACCAAAGCGCGCTGCCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCG                                            <  SRR4302532.81842/94‑1 (MQ=60)
                                                   ACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGG                                           >  SRR4302532.5577131/1‑94 (MQ=60)
                                                   ACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGG                                           <  SRR4302532.4615415/94‑1 (MQ=60)
                                                     CAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGGAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGC                                         >  SRR4302532.6293543/1‑94 (MQ=60)
                                                     CAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGC                                         >  SRR4302532.6314750/1‑94 (MQ=60)
                                                     CAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGC                                         <  SRR4302532.984204/94‑1 (MQ=60)
                                                      aaagagagctccccGCGGCCCCTGCGCCTAGATAGCAGCACGCATGGTTTTAACCAGATCGTAAGcttcgaagaatcgcttaacgtcgagggca                                        >  SRR4302532.4426870/15‑65 (MQ=47)
                                                       AAGCGCGCTACCAGCGGCCCCAGCGCCGAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAT                                       <  SRR4302532.2624677/94‑1 (MQ=60)
                                                       AAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAT                                       <  SRR4302532.4613461/94‑1 (MQ=60)
                                                        AGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATC                                      >  SRR4302532.4209465/1‑94 (MQ=60)
                                                         GCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCA                                     <  SRR4302532.4125603/94‑1 (MQ=60)
                                                            CGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATA                                  <  SRR4302532.4368512/94‑1 (MQ=60)
                                                                ACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCA                              <  SRR4302532.535612/94‑1 (MQ=60)
                                                                ACCAGCGCCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCA                              <  SRR4302532.5136049/94‑1 (MQ=60)
                                                                 CCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCAC                             >  SRR4302532.5579576/1‑94 (MQ=60)
                                                                  CAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACA                            <  SRR4302532.565806/94‑1 (MQ=60)
                                                                   AGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAG                           >  SRR4302532.2959422/1‑94 (MQ=60)
                                                                    GCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGA                          <  SRR4302532.319597/94‑1 (MQ=60)
                                                                      GGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAC                        >  SRR4302532.628611/1‑94 (MQ=60)
                                                                      GGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAC                        <  SRR4302532.3384427/94‑1 (MQ=60)
                                                                      GGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAC                        >  SRR4302532.5663968/1‑94 (MQ=60)
                                                                        CCCCAGCGCCCAGATAGAAGCGCTCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCA                      <  SRR4302532.3611345/94‑1 (MQ=60)
                                                                        CCCCAGCGCACAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCGCAGAACCA                      <  SRR4302532.2717357/94‑1 (MQ=60)
                                                                         CCCGGCGCCCAGGTAGAAGTGCGCATGGTTTTAAGCAGATCGTGAGGGGGGCAGAATACATTAACGTCGCGGGCAGCAGTATGAACAGAACgct                     >  SRR4302532.3738788/1‑91 (MQ=29)
                                                                          CCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATT                    <  SRR4302532.1275542/94‑1 (MQ=60)
                                                                           CAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTA                   <  SRR4302532.1336891/94‑1 (MQ=60)
                                                                            AGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAC                  <  SRR4302532.3884852/94‑1 (MQ=60)
                                                                            AGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAC                  <  SRR4302532.3884853/94‑1 (MQ=60)
                                                                            AGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAC                  >  SRR4302532.5049048/1‑94 (MQ=60)
                                                                            AGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAC                  >  SRR4302532.2725934/1‑94 (MQ=60)
                                                                            AGCGCCCAGATAGAAGCGCGCATGGTTTTAACAAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAC                  <  SRR4302532.926568/94‑1 (MQ=60)
                                                                             cccCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTTCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACG                 <  SRR4302532.4076840/91‑1 (MQ=60)
                                                                             GCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACG                 <  SRR4302532.1052750/94‑1 (MQ=60)
                                                                              CGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGT                >  SRR4302532.4882945/1‑94 (MQ=60)
                                                                              CGCCCAGATAGAAGCGCGCATGGGTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGT                >  SRR4302532.3873113/1‑94 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTT               >  SRR4302532.4747656/1‑94 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTT               >  SRR4302532.1572009/1‑94 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTT               >  SRR4302532.646218/1‑94 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTT               <  SRR4302532.3422515/94‑1 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTT               <  SRR4302532.3422516/94‑1 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTT               >  SRR4302532.5154994/1‑94 (MQ=60)
                                                                               GCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACGGAACCATTACGTT               >  SRR4302532.1884502/1‑94 (MQ=60)
                                                                                 CCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCT             >  SRR4302532.2682329/1‑94 (MQ=60)
                                                                                 CCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCT             >  SRR4302532.4289334/1‑94 (MQ=60)
                                                                                   AGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTAC           >  SRR4302532.5767989/1‑94 (MQ=60)
                                                                                      TAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTT        <  SRR4302532.3129163/94‑1 (MQ=60)
                                                                                      TAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTT        <  SRR4302532.946003/94‑1 (MQ=60)
                                                                                      AAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTT        <  SRR4302532.2237473/94‑1 (MQ=60)
                                                                                       AGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTC       <  SRR4302532.4018441/94‑1 (MQ=60)
                                                                                       AGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTC       <  SRR4302532.5141098/94‑1 (MQ=60)
                                                                                            CGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACC  >  SRR4302532.353410/1‑94 (MQ=60)
                                                                                             |                                                                                            
AACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACC  >  NC_000913/3333770‑3333955

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