Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,375 +A coding (458/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3721.A99.0% 346.5 / ‑3.8 96coding (455/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (32/63);  total (33/63)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449
                                                                                              |                                                                                               
gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                                  <  2:341848/93‑1 (MQ=255)
gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTc                                                                                                  <  2:4402295/93‑1 (MQ=255)
 cTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCa                                                                                                >  2:2326490/1‑93 (MQ=255)
 cTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCa                                                                                                >  1:5244410/1‑93 (MQ=255)
  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaa                                                                                               >  2:2027251/1‑93 (MQ=255)
  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaa                                                                                               <  1:472365/93‑1 (MQ=255)
  tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaa                                                                                               <  1:3256673/93‑1 (MQ=255)
   agAGCAGGCGATTAACGCTCGTTGGTGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaaa                                                                                              <  2:1074239/92‑1 (MQ=255)
    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaaaa                                                                                             <  1:2285683/93‑2 (MQ=255)
    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaaaa                                                                                             <  2:4163007/93‑2 (MQ=255)
    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCaaaa                                                                                             >  2:2686103/1‑92 (MQ=255)
    gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCaaaa                                                                                             >  2:4295461/1‑92 (MQ=255)
     aGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTTaaaag                                                                                            >  1:3192920/1‑91 (MQ=255)
       cAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaaaaggc                                                                                          >  2:1634344/1‑89 (MQ=255)
       cAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCaaaaggc                                                                                          <  2:5700231/93‑5 (MQ=255)
         ggCGATTAACGCTCGTTGGGGCTATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                         <  1:6112872/93‑1 (MQ=255)
         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                         <  1:1506873/93‑1 (MQ=255)
         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                         <  2:4951341/93‑1 (MQ=255)
         ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCa                                                                                         <  1:5722761/93‑1 (MQ=255)
          gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAt                                                                                        <  2:1599366/93‑1 (MQ=255)
          gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAt                                                                                        >  2:4819971/1‑93 (MQ=255)
           cGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATg                                                                                       <  2:2169652/93‑1 (MQ=255)
           cGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATg                                                                                       <  1:4269990/93‑1 (MQ=255)
            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                      <  1:2157699/93‑1 (MQ=255)
            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                      <  2:5492824/93‑1 (MQ=255)
            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                      >  1:5172308/1‑93 (MQ=255)
            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                      <  1:3756419/93‑1 (MQ=255)
            gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgc                                                                                      <  1:2478711/93‑1 (MQ=255)
             aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAcgcg                                                                                     >  2:1385819/1‑93 (MQ=255)
             aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcg                                                                                     >  1:1442960/1‑93 (MQ=255)
              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATgcgc                                                                                    >  1:5228684/1‑93 (MQ=255)
              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATgcgc                                                                                    >  1:5071661/1‑93 (MQ=255)
              ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATgcgc                                                                                    >  2:5994215/1‑93 (MQ=255)
                aaCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTg                                                                                  <  2:2092422/93‑1 (MQ=255)
                 aCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTgg                                                                                 <  2:2390309/93‑1 (MQ=255)
                  cGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTTCAAAAGACATGCGCTTGa                                                                                >  2:557163/1‑93 (MQ=255)
                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGaa                                                                               <  2:1607153/93‑1 (MQ=255)
                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGaa                                                                               >  1:2885126/1‑93 (MQ=255)
                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGaa                                                                               <  1:2474865/93‑1 (MQ=255)
                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGaa                                                                               <  1:868213/93‑1 (MQ=255)
                   gCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGaa                                                                               <  1:926211/93‑1 (MQ=255)
                    cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                             <  1:5707178/93‑1 (MQ=255)
                      cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                           <  2:4574996/93‑3 (MQ=255)
                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGGTTCAAAAGGCCATGCGCTGGAAgcgg                                                                          >  1:4813315/1‑90 (MQ=255)
                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                          <  1:4623304/93‑4 (MQ=255)
                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                          <  2:3083575/93‑4 (MQ=255)
                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                          <  2:3097571/93‑4 (MQ=255)
                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                          <  2:3282794/93‑4 (MQ=255)
                       gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCATGCGCTGGAAgcgg                                                                          >  1:4041485/1‑90 (MQ=255)
                        ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                         <  2:3738063/93‑5 (MQ=255)
                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                >  1:882407/1‑93 (MQ=255)
                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                <  2:722109/93‑1 (MQ=255)
                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                >  2:2003960/1‑93 (MQ=255)
                             gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGGCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                >  2:4689360/1‑93 (MQ=255)
                              cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                               <  2:2103880/93‑1 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                             <  2:5143571/93‑3 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                             <  2:219070/93‑3 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                             <  2:4328059/93‑3 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                             <  2:3925526/93‑3 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                             >  1:4242722/1‑91 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                             >  1:3268770/1‑91 (MQ=255)
                                aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCACTGGCTGGAAGGGGTGGCG‑‑‑‑aaaa                                                             >  1:1410130/1‑91 (MQ=255)
                                 tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                            >  1:2474435/1‑90 (MQ=255)
                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                           <  1:3617388/93‑5 (MQ=255)
                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                           <  2:2759649/93‑5 (MQ=255)
                                                                  tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                        <  1:3067422/93‑1 (MQ=255)
                                                                  tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                        >  2:5367497/1‑93 (MQ=255)
                                                                   cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                       <  1:3423220/93‑2 (MQ=255)
                                                                   cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                       >  2:2690426/1‑92 (MQ=255)
                                                                   cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                       >  2:2607369/1‑92 (MQ=255)
                                                                    tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  2:4708075/1‑91 (MQ=255)
                                                                    tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      <  1:1051706/93‑3 (MQ=255)
                                                                    tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  1:5951006/1‑91 (MQ=255)
                                                                    tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  1:5995372/1‑91 (MQ=255)
                                                                    tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                      >  2:1156237/1‑91 (MQ=255)
                                                                     ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                     >  2:630739/1‑90 (MQ=255)
                                                                             tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa             <  1:100493/93‑1 (MQ=255)
                                                                             tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa             <  2:349093/93‑1 (MQ=255)
                                                                             tGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa             >  2:6045798/1‑93 (MQ=255)
                                                                              ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt            <  2:2148939/93‑1 (MQ=255)
                                                                              ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAAACCGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt            <  2:1686766/93‑1 (MQ=255)
                                                                               gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg           <  2:2236559/93‑1 (MQ=255)
                                                                               gcgggcggCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg           <  2:790947/93‑1 (MQ=255)
                                                                                cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGa          <  1:5478423/93‑1 (MQ=255)
                                                                                 gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa         <  1:2635131/93‑1 (MQ=255)
                                                                                 gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa         <  2:3445771/93‑1 (MQ=255)
                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        <  2:2772791/93‑1 (MQ=255)
                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        <  1:4957056/93‑1 (MQ=255)
                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        >  1:3993000/1‑93 (MQ=255)
                                                                                  ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac        <  2:331868/93‑1 (MQ=255)
                                                                                   gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg       >  1:4814006/1‑93 (MQ=255)
                                                                                   gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg       <  1:2957114/93‑1 (MQ=255)
                                                                                     ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAacgac     >  2:5773184/1‑93 (MQ=255)
                                                                                     ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     <  2:4498313/93‑1 (MQ=255)
                                                                                     ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     >  2:5281397/1‑93 (MQ=255)
                                                                                     ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     <  2:4279944/93‑1 (MQ=255)
                                                                                     ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac     >  2:3119511/1‑93 (MQ=255)
                                                                                      gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg    >  2:4300379/1‑93 (MQ=255)
                                                                                       cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc   >  2:5938791/1‑93 (MQ=255)
                                                                                         gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc  <  1:4418546/93‑1 (MQ=255)
                                                                                              |                                                                                               
GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008282‑4008458
                                                                                               |                                                                                             
TGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAA                                                                                               <  SRR4302533.475190/94‑1 (MQ=60)
TGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAA                                                                                               <  SRR4302533.3276758/94‑1 (MQ=60)
  GAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTTAAAAG                                                                                             >  SRR4302533.3212560/1‑94 (MQ=60)
  GAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAG                                                                                             <  SRR4302533.2299831/94‑1 (MQ=60)
  GAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAG                                                                                             <  SRR4302533.5679763/94‑1 (MQ=60)
   AGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGTGTTTCATAAGG                                                                                             >  SRR4302533.781467/1‑94 (MQ=60)
     CAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAaagcc                                                                                          >  SRR4302533.4589789/1‑89 (MQ=60)
      AGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAaagcca                                                                                         >  SRR4302533.38905/1‑88 (MQ=60)
       GGCGATTAACGCTCGTTGGGGCTATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAT                                                                                         <  SRR4302533.6149765/94‑1 (MQ=60)
       GGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAT                                                                                         <  SRR4302533.1515956/94‑1 (MQ=60)
       GGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAT                                                                                         <  SRR4302533.5757314/94‑1 (MQ=60)
       GGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAT                                                                                         <  SRR4302533.4482448/94‑1 (MQ=60)
        GCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAaatcaagg                                                                                       >  SRR4302533.1434716/1‑86 (MQ=60)
        GCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAaggccctg                                                                                       >  SRR4302533.1714856/1‑86 (MQ=60)
         CGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGC                                                                                       <  SRR4302533.4295693/94‑1 (MQ=60)
         CGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGC                                                                                       >  SRR4302533.5203517/1‑94 (MQ=60)
          GATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCG                                                                                      >  SRR4302533.394157/1‑94 (MQ=60)
          GATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCG                                                                                      <  SRR4302533.2170981/94‑1 (MQ=60)
          GATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCG                                                                                      <  SRR4302533.2493997/94‑1 (MQ=60)
          GATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCG                                                                                      <  SRR4302533.3779444/94‑1 (MQ=60)
          GATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCG                                                                                      >  SRR4302533.1451666/1‑94 (MQ=60)
           ATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGC                                                                                     >  SRR4302533.5260139/1‑94 (MQ=60)
           ATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGC                                                                                     >  SRR4302533.5102315/1‑94 (MQ=60)
              AACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAaagccacgagctgg                                                                                 >  SRR4302533.4806168/1‑80 (MQ=60)
                CGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA                                                                                <  SRR4302533.5075433/94‑1 (MQ=60)
                CGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAA                                                                                >  SRR4302533.2902864/1‑94 (MQ=60)
                 GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAG                                                                               <  SRR4302533.931646/94‑1 (MQ=60)
                 GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAG                                                                               <  SRR4302533.2490125/94‑1 (MQ=60)
                 GCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAG                                                                               <  SRR4302533.873300/94‑1 (MQ=60)
                  CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGC                                                                             <  SRR4302533.5741651/94‑1 (MQ=60)
                  ATCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAaagccctgcgctggaagc                                                                             >  SRR4302533.777670/1‑76 (MQ=60)
                    CGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGGTTCAAAAGGCCATGCGCTGGAAGcgg                                                                            >  SRR4302533.4842487/1‑91 (MQ=60)
                    CGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgg                                                                            <  SRR4302533.3215368/94‑4 (MQ=60)
                    CGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCATGCGCTGGAAGcgg                                                                            >  SRR4302533.4066083/1‑91 (MQ=60)
                     GTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggt                                                                           <  SRR4302533.4651328/94‑5 (MQ=60)
                     GTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggt                                                                           <  SRR4302533.2032402/94‑5 (MQ=60)
                      TTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtg                                                                          >  SRR4302533.3893437/1‑89 (MQ=60)
                       TGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtgg                                                                         <  SRR4302533.4178070/94‑7 (MQ=60)
                        GGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggc                                                                        >  SRR4302533.2294013/1‑87 (MQ=60)
                        GGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggc                                                                        <  SRR4302533.2958628/94‑8 (MQ=60)
                         GGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcg                                                                       <  SRR4302533.3340360/94‑9 (MQ=60)
                          GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑A                                                                  >  SRR4302533.887590/1‑94 (MQ=60)
                          GGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑A                                                                  <  SRR4302533.5337889/94‑1 (MQ=60)
                            CGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggggaag                                                                    >  SRR4302533.1620270/1‑83 (MQ=60)
                             GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaaga                                                                   >  SRR4302533.4268339/1‑82 (MQ=60)
                             GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaaga                                                                   >  SRR4302533.3288937/1‑82 (MQ=60)
                             GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAaggccactggctggaaggggtggcgaaaa                                                                  >  SRR4302533.1418657/1‑65 (MQ=60)
                             GATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCGAAAGGCCATGCGCTGGAAGcgggggcggaga                                                                   >  SRR4302533.5631381/1‑82 (MQ=60)
                              ATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcggagaa                                                                  >  SRR4302533.1150890/1‑81 (MQ=60)
                              ATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaa                                                                  >  SRR4302533.2489694/1‑81 (MQ=60)
                                CGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGGTTCAAAAAGCCATGCGCTGGAAGcggtggcgaagaaac                                                                >  SRR4302533.4661568/1‑79 (MQ=60)
                                CGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaac                                                                >  SRR4302533.4570836/1‑79 (MQ=60)
                                CGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaac                                                                <  SRR4302533.3639640/94‑16 (MQ=60)
                                 GCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggcgaagaaact                                                               >  SRR4302533.1959643/1‑78 (MQ=60)
                                  CGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGGTTCAAAAGGCCATGCGCTGGAAGGgggggggaagcaaccg                                                             >  SRR4302533.4074910/1‑78 (MQ=60)
                                   GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                             <  SRR4302533.2181474/94‑19 (MQ=60)
                                   GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                             >  SRR4302533.5196848/1‑76 (MQ=60)
                                   GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                             <  SRR4302533.247755/94‑19 (MQ=60)
                                   GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                             <  SRR4302533.1419427/94‑19 (MQ=60)
                                   GTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgg                                                             <  SRR4302533.4252206/94‑19 (MQ=60)
                                    TCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactggg                                                            <  SRR4302533.1253243/94‑20 (MQ=60)
                                     CAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggc                                                           >  SRR4302533.2600926/1‑74 (MQ=60)
                                     CAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcgggggcgaagaaactgggc                                                           >  SRR4302533.2040920/1‑74 (MQ=60)
                                      AACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggct                                                          >  SRR4302533.4416745/1‑73 (MQ=60)
                                      AACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggct                                                          >  SRR4302533.4062766/1‑73 (MQ=60)
                                           CAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagc                                                     >  SRR4302533.2545756/1‑68 (MQ=60)
                                            AGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcc                                                    >  SRR4302533.4012320/1‑67 (MQ=60)
                                            AGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcc                                                    <  SRR4302533.1574414/94‑28 (MQ=60)
                                            AGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctaaagcc                                                    >  SRR4302533.214174/1‑67 (MQ=60)
                                             GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcct                                                   >  SRR4302533.1757581/1‑66 (MQ=60)
                                             GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcct                                                   <  SRR4302533.4606144/94‑29 (MQ=60)
                                             GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcct                                                   <  SRR4302533.6253719/94‑29 (MQ=60)
                                             GTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTGCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcct                                                   <  SRR4302533.2381622/94‑29 (MQ=60)
                                              TTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctg                                                  >  SRR4302533.47652/1‑65 (MQ=60)
                                               TTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctga                                                 >  SRR4302533.3195947/1‑64 (MQ=60)
                                               TTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctga                                                 >  SRR4302533.2601426/1‑64 (MQ=60)
                                               TTCTCTACCTTAACCTGTCCGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctga                                                 <  SRR4302533.2760742/94‑31 (MQ=60)
                                                TCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaa                                                <  SRR4302533.5359189/94‑32 (MQ=60)
                                                 CTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaag                                               >  SRR4302533.5244329/1‑62 (MQ=60)
                                                 CTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaaatgggctacagcctgaag                                               >  SRR4302533.2848517/1‑62 (MQ=60)
                                                  TCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaagg                                              <  SRR4302533.4555648/94‑34 (MQ=60)
                                                     ACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggatt                                           >  SRR4302533.4602454/1‑58 (MQ=60)
                                                      CCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGaaagcggtggcgaagaaactgggctacagcctgaaggattg                                          <  SRR4302533.396174/94‑42 (MQ=60)
                                                      CCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattg                                          <  SRR4302533.272655/94‑38 (MQ=60)
                                                       CTTAGCCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgt                                         <  SRR4302533.674543/94‑39 (MQ=60)
                                                       CTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCTCTGGAAGcggtggcgaagaaactggtctacagcctgaaggattgt                                         >  SRR4302533.3921045/1‑56 (MQ=60)
                                                       CTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgt                                         <  SRR4302533.1374313/94‑39 (MQ=60)
                                                        TTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgta                                        <  SRR4302533.4840798/94‑40 (MQ=60)
                                                           ACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagactgaaggattgtattg                                     >  SRR4302533.7657/1‑52 (MQ=60)
                                                             CTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcg                                   <  SRR4302533.2486450/94‑45 (MQ=60)
                                                             CTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcg                                   <  SRR4302533.4946983/94‑45 (MQ=60)
                                                             CTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcg                                   <  SRR4302533.3164330/94‑45 (MQ=60)
                                                              TGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgt                                  >  SRR4302533.3979667/1‑49 (MQ=60)
                                                                TCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgttt                                <  SRR4302533.3086262/94‑48 (MQ=60)
                                                                 CTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttg                               >  SRR4302533.6031684/1‑46 (MQ=60)
                                                                 CTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttg                               >  SRR4302533.5987001/1‑46 (MQ=60)
                                                                 CTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttg                               <  SRR4302533.3444365/94‑49 (MQ=60)
                                                                  TGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactggggtacagcctgaaggattgtattgcgtttgg                              >  SRR4302533.1478632/1‑45 (MQ=60)
                                                                  TGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttgg                              <  SRR4302533.1057898/94‑50 (MQ=60)
                                                                      AGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgac                          <  SRR4302533.714606/94‑54 (MQ=60)
                                                                       GTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGcggtggcgaagaaactgggctacagcctgaaggattgtattgcgtttggtgacg                         <  SRR4302533.2403475/94‑55 (MQ=60)
                                                                                         tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA       >  SRR4302533.4854775/4‑94 (MQ=16)
                                                                                         tcaAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAA       >  SRR4302533.4854776/4‑94 (MQ=16)
                                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT      <  SRR4302533.5801257/92‑1 (MQ=18)
                                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT      >  SRR4302533.5860484/3‑94 (MQ=18)
                                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT      <  SRR4302533.3455385/92‑1 (MQ=18)
                                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT      <  SRR4302533.1001159/92‑1 (MQ=18)
                                                                                          caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT      <  SRR4302533.4581160/92‑1 (MQ=18)
                                                                                          AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG     >  SRR4302533.292135/1‑94 (MQ=19)
                                                                                              AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT   <  SRR4302533.6229189/94‑1 (MQ=21)
                                                                                              AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT   <  SRR4302533.2953665/94‑1 (MQ=21)
                                                                                               AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG  >  SRR4302533.2757655/1‑94 (MQ=23)
                                                                                               AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG  <  SRR4302533.3758008/94‑1 (MQ=23)
                                                                                               AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG  <  SRR4302533.1204246/94‑1 (MQ=23)
                                                                                               AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG  >  SRR4302533.49854/1‑94 (MQ=23)
                                                                                               |                                                                                             
TGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008282‑4008458

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