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
A2 F36 I0 R1 80 95.4 3476014 97.7% 3396065 133.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,182,929 A→T 100% E562V (GAA→GTA)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,9290AT95.1% 337.7 / 5.2 101E562V (GAA→GTA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (51/45);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 103 total aligned reads displayed.

TTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGC  >  NC_000913/4182784‑4183067
                                                                                                                                                 |                                                                                                                                          
tttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa                                                                                                                                          >  3:16225/1‑148 (MQ=255)
 ttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGtaa                                                                                                                                          <  6:211378/147‑1 (MQ=255)
     ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                                                                     >  1:152005/1‑148 (MQ=255)
      gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCt                                                                                                                                     >  3:268838/1‑147 (MQ=255)
      gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCt                                                                                                                                     >  4:263986/1‑147 (MQ=255)
          aGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAgg                                                                                                                                >  7:47121/1‑148 (MQ=255)
             acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGt                                                                                                                               >  8:189122/1‑146 (MQ=255)
                acCCGCTGTCTGATATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   <  5:94590/139‑1 (MQ=255)
                acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGaa                                                                                                                          >  3:385768/1‑148 (MQ=255)
                 cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGATACGTACACACGACTCAATACGGTCGCGTATTTCAAATCGTAACCCCTGAAgagatatg                                                                                                                          >  6:94590/1‑140 (MQ=255)
                 cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGaa                                                                                                                          >  1:323886/1‑147 (MQ=255)
                        tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGACCAGGCGGTCTGACCCGGGAACTTTCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCgg                                                                                                                    <  5:430556/146‑1 (MQ=255)
                         cTGAGATTACGCACAAACGTCGTATCTCCGCACACGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGCCGTACAACAGACTCACTACCGTCGCGTATGTCCAATCGTAACCCCTGAAGTTCCGAACATCGGTCt                                                                                                                 >  5:223790/1‑148 (MQ=255)
                          tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCTCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCt                                                                                                                 >  8:142406/1‑147 (MQ=255)
                            agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGACTGACCCGTGAAAGAGCATGCTTCGAAGTTCGACACATACACCCGACTCACTACGGTCGAGTATGTCCAATCGTAACCCATGAAGGTCCGAACATCGGtatta                                                                                                               >  2:66355/1‑143 (MQ=255)
                              aTTACGCACAACCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGATCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAGGTCCAATCGTAACCCCGGAAGGTCCGACCATCGGTCGGagc                                                                                                             >  8:161331/1‑145 (MQ=255)
                               ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGCACGTGCAGGCTTCGAAGTTCGAGCCGGACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCaa                                                                                                           <  2:131233/148‑1 (MQ=255)
                                tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCTCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   >  8:160475/1‑123 (MQ=255)
                                tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   <  7:160475/123‑1 (MQ=255)
                                tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAAc                                                                                                          >  1:59015/1‑148 (MQ=255)
                                    cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAATATCgg                                                                                                                    >  3:188084/1‑134 (MQ=255)
                                    cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAATATCgg                                                                                                                    <  4:188084/134‑1 (MQ=255)
                                      caAACGTCGTATCTCCGCACTCGGCACATGCGGTCTGACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCCTATGTCCAATCGAAACCCCATAAGCTCCGAACAACTGTCTGATCAACTCTCTg                                                                                                    >  5:325660/1‑148 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   >  7:234352/1‑115 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   <  8:234352/115‑1 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                                                                  <  4:180691/148‑1 (MQ=255)
                                          cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCg                                                                                                <  8:175955/148‑1 (MQ=255)
                                               tATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                                                                           <  8:47121/148‑1 (MQ=255)
                                                aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCAGTGTACg                                                                                          >  4:70305/1‑148 (MQ=255)
                                                  ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                                                                         <  2:323886/147‑1 (MQ=255)
                                                    ccGCACTCGGCCCAGGCGGTCTTACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGAACACCCGACTCACTACGGTCGCGTATGTCAAATCGTAACCCCTGAAGGTCCGAACATCTGTCTGCACAACTCTCTGTCCGTGTACGcaca                                                                                      >  7:261680/1‑148 (MQ=255)
                                                    ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                                                        >  3:275212/1‑146 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  1:243684/146‑1 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  2:59015/146‑1 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  5:336745/146‑1 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACACGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                                                                         <  3:263986/144‑1 (MQ=255)
                                                       cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                                                                                   >  7:13418/1‑148 (MQ=255)
                                                           cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                               >  7:385781/1‑148 (MQ=255)
                                                            ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                                                                <  1:384706/146‑1 (MQ=255)
                                                            ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                                                                <  7:204339/146‑1 (MQ=255)
                                                            ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                                                              <  5:176547/148‑1 (MQ=255)
                                                             gctcAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                                                             <  7:142406/145‑1 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                                                            >  8:272703/1‑147 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                                                           <  4:385768/148‑1 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                                                           <  4:16225/148‑1 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCTTATGTCCAATCGTAACCCCTTAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                                                           <  7:384037/148‑1 (MQ=255)
                                                                cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                                                          <  7:189122/148‑1 (MQ=255)
                                                                cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                                                          >  4:342652/1‑148 (MQ=255)
                                                                  ggCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                                                        >  8:34779/1‑148 (MQ=255)
                                                                   gCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTACCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                                                       >  7:432086/1‑148 (MQ=255)
                                                                      gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                                                    >  6:54693/1‑148 (MQ=255)
                                                                       tCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            <  3:116766/107‑1 (MQ=255)
                                                                       tCTGACCCGAGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            >  4:116766/1‑107 (MQ=255)
                                                                       gcTGACCCGTGAACGTGCAGGCTTCGAAGTCCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                                                                   <  8:385781/147‑1 (MQ=255)
                                                                        cTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                                                                  >  1:419717/1‑148 (MQ=255)
                                                                          gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  3:314385/1‑148 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                                                                                 <  3:165095/130‑1 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                                                                                 >  4:165095/1‑130 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                                                               <  2:17625/148‑1 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                                                               <  4:268838/148‑1 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            <  5:41380/102‑1 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            >  6:41380/1‑102 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  3:23911/1‑146 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  2:68618/1‑146 (MQ=255)
                                                                                    cGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                                                       >  5:283253/1‑147 (MQ=255)
                                                                                       gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                                                             <  3:142656/106‑1 (MQ=255)
                                                                                       gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                                                             >  4:142656/1‑106 (MQ=255)
                                                                                           gCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa                                               <  2:152005/148‑1 (MQ=255)
                                                                                           gCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACCATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa                                               <  3:338975/148‑1 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  4:129635/1‑130 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                <  3:129635/130‑1 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              >  4:186218/1‑148 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              <  4:275212/148‑1 (MQ=255)
                                                                                              tCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt                                            >  8:121316/1‑148 (MQ=255)
                                                                                              tCGAAGTTAGATACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCTTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCATGAGACTCCGTATCGTaaa                                              >  1:311373/1‑146 (MQ=255)
                                                                                                gAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa                                          <  4:314385/148‑1 (MQ=255)
                                                                                                    ttCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGa                                      >  2:121037/1‑148 (MQ=255)
                                                                                                       gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACgg                                   <  3:70305/148‑1 (MQ=255)
                                                                                                            gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGtgttg                              <  3:342652/148‑1 (MQ=255)
                                                                                                             tACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              >  2:329676/1‑131 (MQ=255)
                                                                                                             tACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              <  1:329676/131‑1 (MQ=255)
                                                                                                               cacCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTa                            <  6:283253/147‑1 (MQ=255)
                                                                                                                  ccGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGGACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATAGTAAAGTGACCGACCGTGTTGTAACt                         >  5:210887/1‑147 (MQ=255)
                                                                                                                   cGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              <  7:325637/125‑1 (MQ=255)
                                                                                                                   cGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              >  8:325637/1‑125 (MQ=255)
                                                                                                                   cGACTCACAACGGTCGCGAATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCACCTCTCTCTCCGTGTACGCACACACGAACGAATACGGGTTCCTTGCGACTCCGTATCGTAAATTGCCCCACGGTGTTGTAACTGa                       >  6:181742/1‑148 (MQ=255)
                                                                                                                         aCTACGGTCGCGTATGTACAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAt                 <  7:392348/148‑1 (MQ=255)
                                                                                                                          cTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                >  8:106435/1‑148 (MQ=255)
                                                                                                                          cTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                <  3:186218/148‑1 (MQ=255)
                                                                                                                           tACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTaa                           <  7:6287/136‑1 (MQ=255)
                                                                                                                           tACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                >  4:299667/1‑147 (MQ=255)
                                                                                                                             cGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTc               <  5:54693/146‑1 (MQ=255)
                                                                                                                             cGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTc               <  4:178048/146‑1 (MQ=255)
                                                                                                                                 cgcgTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTACCGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAcc         >  5:112802/1‑148 (MQ=255)
                                                                                                                                      aTGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       >  6:422406/1‑65 (MQ=255)
                                                                                                                                      aTGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  5:422406/65‑1 (MQ=255)
                                                                                                                                       tGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTg   <  1:68618/148‑1 (MQ=255)
                                                                                                                                         tCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGc  >  2:195561/1‑147 (MQ=255)
                                                                                                                                          ccAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc          >  5:361221/1‑138 (MQ=255)
                                                                                                                                          ccAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc          <  6:361221/138‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
TTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGC  >  NC_000913/4182784‑4183067

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAG  >  NC_000913/4182779‑4183078
                                                                                                                                                      |                                                                                                                                                     
CTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGT                                                                                                                                                       >  NB501016:230:HFVJLAFXY:4:11604:3706:11008/1‑151 (MQ=60)
   aaccacgcagcgttacctgatgtgacgccgacaattctcatcatcgctacATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAAC                                                                                                                                                    <  NB501016:230:HFVJLAFXY:1:21306:10901:9035/101‑1 (MQ=60)
   AGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGATACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAC                                                                                                                                                    >  NB501016:230:HFVJLAFXY:2:11103:24151:16826/1‑151 (MQ=60)
   AGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAAC                                                                                                                                                    <  NB501016:230:HFVJLAFXY:3:21509:21032:2542/151‑1 (MQ=60)
    GTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACC                                                                                                                                                   <  NB501016:230:HFVJLAFXY:3:11401:10564:2086/151‑1 (MQ=60)
        ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTG                                                                                                                                               >  NB501016:230:HFVJLAFXY:1:11301:2439:7471/1‑151 (MQ=60)
         TGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTG                                                                                                                                               >  NB501016:230:HFVJLAFXY:2:21111:2616:18081/1‑150 (MQ=60)
             CCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGT                                                                                                                                          >  NB501016:230:HFVJLAFXY:4:11408:5805:6289/1‑151 (MQ=60)
              CAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTC                                                                                                                                         >  NB501016:230:HFVJLAFXY:2:21309:5063:12028/1‑151 (MQ=60)
                   CAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAAC                                                                                                                                    >  NB501016:230:HFVJLAFXY:2:21308:2115:8581/1‑151 (MQ=60)
                    AACCCGCTGTCTGATATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAG                                                                                                                                            <  NB501016:230:HFVJLAFXY:3:11505:14390:1032/142‑1 (MQ=60)
                    AACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACC                                                                                                                                    >  NB501016:230:HFVJLAFXY:1:21204:5597:15891/1‑150 (MQ=60)
                            GTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGACCAGGCGGTCTGACCCGGGAACTTTCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTC                                                                                                                             <  NB501016:230:HFVJLAFXY:3:21612:11772:14394/149‑1 (MQ=60)
                            GTCTGAGATTACGCACAAACGTCGTATCTCCGCACACGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGCCGTACAACAGACTCACTACCGTCGCGTATGTCCAATCGTAACCCCTGAAGTTCCGAACATCGGTCTG                                                                                                                           >  NB501016:230:HFVJLAFXY:3:21402:10958:2779/1‑151 (MQ=60)
                                   ATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACT                                                                                                                    >  NB501016:230:HFVJLAFXY:1:11109:9171:16408/1‑151 (MQ=60)
                                   ATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAG                                                                                                                                            <  NB501016:230:HFVJLAFXY:4:11601:4514:19398/127‑1 (MQ=60)
                                       CGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAATATCGGTC                                                                                                                             >  NB501016:230:HFVJLAFXY:2:11309:15910:17374/1‑138 (MQ=60)
                                         CACAAACGTCGTATCTCCGCACTCGGCACATGCGGTCTGACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCCTATGTCCAATCGAAACCCCATAAGCTCCGAACAACTGTCTGATCAACTCTCTGT                                                                                                              >  NB501016:230:HFVJLAFXY:3:21507:11445:19306/1‑151 (MQ=60)
                                           CAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAG                                                                                                                                            >  NB501016:230:HFVJLAFXY:4:21402:24629:5016/1‑119 (MQ=60)
                                                       CTCCGCACTCGGCCCAGGCGGTCTTACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGAACACCCGACTCACTACGGTCGCGTATGTCAAATCGTAACCCCTGAAGGTCCGAACATCTGTCTGCACAACTCTCTGTCCGTGTACGCACAC                                                                                                >  NB501016:230:HFVJLAFXY:4:21406:7504:11741/1‑151 (MQ=60)
                                                       CTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACCTCGGTCTGATCAACTCTCTGTCCGTGTACGCAC                                                                                                  >  NB501016:230:HFVJLAFXY:2:21201:8218:2891/1‑149 (MQ=60)
                                                         CCGCACTCGGCCCAGGCTGTCTGACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACACGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCAC                                                                                                  <  NB501016:230:HFVJLAFXY:2:21110:5260:20209/147‑1 (MQ=60)
                                                         CCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAG                                                                                                <  NB501016:230:HFVJLAFXY:1:21103:18265:19782/149‑1 (MQ=60)
                                                         CCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAG                                                                                                <  NB501016:230:HFVJLAFXY:3:21509:18031:16747/149‑1 (MQ=60)
                                                          CGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACT                                                                                             >  NB501016:230:HFVJLAFXY:4:11403:14162:3959/1‑151 (MQ=60)
                                                              CTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACG                                                                                         >  NB501016:230:HFVJLAFXY:4:21602:10296:19969/1‑151 (MQ=60)
                                                                CGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAA                                                                                       <  NB501016:230:HFVJLAFXY:3:11606:6883:9783/151‑1 (MQ=60)
                                                                CGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACG                                                                                         <  NB501016:230:HFVJLAFXY:4:11609:1787:5300/149‑1 (MQ=60)
                                                                CGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACG                                                                                         <  NB501016:230:HFVJLAFXY:1:21303:10502:1662/149‑1 (MQ=60)
                                                                 GGCTCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAT                                                                                      <  NB501016:230:HFVJLAFXY:4:11511:1712:2833/151‑1 (MQ=60)
                                                                   CCCAGGCGGTCTGACCCGTGAACGTGCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCTTATGTCCAATCGTAACCCCTTAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAC                                                                                    <  NB501016:230:HFVJLAFXY:4:21602:1500:14602/151‑1 (MQ=60)
                                                                    CCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACG                                                                                   <  NB501016:230:HFVJLAFXY:4:11606:2259:14043/151‑1 (MQ=60)
                                                                      AGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTACCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCCCAGACTAACGAATACGGC                                                                                 >  NB501016:230:HFVJLAFXY:4:21610:13426:16456/1‑151 (MQ=60)
                                                                          GGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAAC                                                                                                                     <  NB501016:230:HFVJLAFXY:2:11209:4341:19141/111‑1 (MQ=60)
                                                                           GTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCT                                                                            >  NB501016:230:HFVJLAFXY:1:21308:23194:16269/1‑151 (MQ=60)
                                                                             CTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTG                                                                          >  NB501016:230:HFVJLAFXY:2:21208:7031:4766/1‑151 (MQ=60)
                                                                              acACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCACTGTCCGTGTACGCACAGACTAACGAAAACGGATTCCTTG                                                                          >  NB501016:230:HFVJLAFXY:3:21510:4078:6546/3‑150 (MQ=60)
                                                                              TGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAC                                                                                          <  NB501016:230:HFVJLAFXY:2:11306:11275:3535/134‑1 (MQ=60)
                                                                               GACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTG                                                                          >  NB501016:230:HFVJLAFXY:2:11105:23176:5339/1‑149 (MQ=60)
                                                                               GACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAAC                                                                                                                     <  NB501016:230:HFVJLAFXY:3:11408:5040:3132/106‑1 (MQ=60)
                                                                                       AACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGT                                                                 >  NB501016:230:HFVJLAFXY:3:21411:14505:18013/1‑150 (MQ=60)
                                                                                          GTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAC                                                                                                      <  NB501016:230:HFVJLAFXY:2:11302:23191:8985/110‑1 (MQ=60)
                                                                                               GGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGA                                                                         <  NB501016:230:HFVJLAFXY:2:11212:18047:5443/134‑1 (MQ=60)
                                                                                               GGCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACCATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAG                                                        <  NB501016:230:HFVJLAFXY:2:21212:14138:11834/151‑1 (MQ=60)
                                                                                                 CTTCGAAGTTAGATACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCTTAACCCCTGAAGGTCCGAACATAGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCATGAGACTCCGTATCGTAAAG                                                        >  NB501016:230:HFVJLAFXY:1:21202:8714:14384/1‑149 (MQ=60)
                                                                                                           CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTG                                            <  NB501016:230:HFVJLAFXY:2:11201:6433:12499/151‑1 (MQ=60)
                                                                                                                CGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTA                                       <  NB501016:230:HFVJLAFXY:2:21301:13778:5033/151‑1 (MQ=60)
                                                                                                                CGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGT                                                       <  NB501016:230:HFVJLAFXY:1:21205:7725:15382/135‑1 (MQ=60)
                                                                                                                     ACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGGACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATAGTAAAGTGACCGACCGTGTTGTAACTG                                   >  NB501016:230:HFVJLAFXY:3:11611:15779:17361/1‑150 (MQ=60)
                                                                                                                      CCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGT                                                       <  NB501016:230:HFVJLAFXY:4:21505:25499:3057/129‑1 (MQ=60)
                                                                                                                             CACTACGGTCGCGTATGTACAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTC                          <  NB501016:230:HFVJLAFXY:4:21604:4271:2333/151‑1 (MQ=60)
                                                                                                                              ACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCA                         <  NB501016:230:HFVJLAFXY:2:11309:1956:11655/151‑1 (MQ=60)
                                                                                                                              ACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACT                                    <  NB501016:230:HFVJLAFXY:4:11402:17702:1402/140‑1 (MQ=60)
                                                                                                                                 ACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCAC                        <  NB501016:230:HFVJLAFXY:3:11410:8529:8060/149‑1 (MQ=60)
                                                                                                                                    GTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTACCGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCT                   >  NB501016:230:HFVJLAFXY:3:11508:5822:3019/1‑151 (MQ=60)
                                                                                                                                         GTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAG                                                                                                <  NB501016:230:HFVJLAFXY:3:21611:1467:9433/69‑1 (MQ=60)
                                                                                                                                           ATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCT            <  NB501016:230:HFVJLAFXY:1:11111:8698:7657/151‑1 (MQ=60)
                                                                                                                                             GTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCT                   >  NB501016:230:HFVJLAFXY:3:21601:14277:18019/1‑142 (MQ=60)
                                                                                                                                                    CGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAA   <  NB501016:230:HFVJLAFXY:1:11208:10425:3690/151‑1 (MQ=60)
                                                                                                                                                     GTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAG  <  NB501016:230:HFVJLAFXY:4:11406:15467:9250/151‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                                     
CTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAG  >  NC_000913/4182779‑4183078

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