Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 744 338.8 17650688 97.8% 17262372 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,235,692 T→C Q258R (CAG→CGG)  ycgB ← putative sporulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,235,6920TC99.6% 858.8 / ‑5.1 252Q258R (CAG→CGG) ycgBputative sporulation protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (127/124);  total (128/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 266 total aligned reads displayed.

AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACC  >  NC_000913/1235602‑1235781
                                                                                          |                                                                                         
aaGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCAcgccgc                                                                                         <  1:1915938/93‑1 (MQ=255)
 aGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTAGCTCACCTTACGCACAATACGCAGGCTTTCACACCGcc                                                                                        >  1:2417115/1‑93 (MQ=255)
  ggTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCa                                                                                       <  1:1026491/93‑1 (MQ=255)
  ggTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCa                                                                                       <  1:3855570/93‑1 (MQ=255)
  ggTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCa                                                                                       >  2:8567952/1‑93 (MQ=255)
   gTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATCCGCAGGATTTCACGCCGCCAt                                                                                      <  1:3594550/93‑1 (MQ=255)
     cGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGa                                                                                    >  1:1208441/1‑93 (MQ=255)
     cGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGa                                                                                    <  1:5469565/93‑1 (MQ=255)
     cGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGa                                                                                    <  2:442212/93‑1 (MQ=255)
      gCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAt                                                                                   <  2:7683784/93‑1 (MQ=255)
        ccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTc                                                                                 >  2:2272901/1‑93 (MQ=255)
         cAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTcc                                                                                <  1:2166681/93‑1 (MQ=255)
         cAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTcc                                                                                >  1:1539388/1‑93 (MQ=255)
          aGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCca                                                                               >  2:3362780/1‑93 (MQ=255)
          aGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCca                                                                               >  1:1381447/1‑93 (MQ=255)
          aGCCTTCGTTCATCACCTGAGGCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCca                                                                               >  2:5028685/1‑93 (MQ=255)
            ccTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagc                                                                             >  1:1672868/1‑93 (MQ=255)
             cTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCaagta                                                                            >  2:6822290/1‑91 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                           >  1:5917709/1‑93 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                           <  2:326639/93‑1 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                           >  1:4354408/1‑93 (MQ=255)
                tgTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTg                                                                         <  2:6311615/92‑1 (MQ=255)
                cGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTg                                                                         <  1:3315023/93‑1 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                        >  1:2158775/1‑93 (MQ=255)
                  ttCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgc                                                                       >  2:8640623/1‑93 (MQ=255)
                  ttCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgc                                                                       <  2:5480644/93‑1 (MQ=255)
                   tcatcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgcg                                                                      >  2:450376/1‑93 (MQ=255)
                     atcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCa                                                                    <  1:1758619/93‑1 (MQ=255)
                     atcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGACGCCATGATTCCAGCAGTGGCGCa                                                                    >  1:6766317/1‑93 (MQ=255)
                      tcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAt                                                                   <  2:540623/93‑1 (MQ=255)
                      tcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAt                                                                   >  2:1882875/1‑93 (MQ=255)
                      tcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAt                                                                   <  2:5724998/93‑1 (MQ=255)
                        aCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAttt                                                                 >  1:6156078/1‑93 (MQ=255)
                         ccTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtttt                                                                >  1:2864093/1‑93 (MQ=255)
                          cTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAttttt                                                               >  1:975599/1‑93 (MQ=255)
                             aGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtctc                                                            <  1:7114390/93‑1 (MQ=255)
                             aGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtctc                                                            >  2:3765027/1‑93 (MQ=255)
                             aGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtctc                                                            <  1:297999/93‑1 (MQ=255)
                              gTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTcc                                                           <  1:6609044/93‑1 (MQ=255)
                              gTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTcc                                                           <  2:5650622/93‑1 (MQ=255)
                              gTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTcc                                                           <  2:1220520/93‑1 (MQ=255)
                               tCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCa                                                          >  2:6711721/1‑93 (MQ=255)
                                cTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCAt                                                         >  1:8311750/1‑93 (MQ=255)
                                 tGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATa                                                        <  2:353531/93‑1 (MQ=255)
                                 tGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATa                                                        >  1:5395797/1‑93 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                       >  2:720621/1‑93 (MQ=255)
                                   tttttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaaa                                                      >  2:8809882/1‑93 (MQ=255)
                                     tttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAt                                                    >  2:6455143/1‑93 (MQ=255)
                                      ttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATa                                                   >  1:7736423/1‑93 (MQ=255)
                                      ttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATa                                                   <  2:571062/93‑1 (MQ=255)
                                      ttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATa                                                   <  2:6024071/93‑1 (MQ=255)
                                      ttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATa                                                   >  1:5598379/1‑93 (MQ=255)
                                          ggATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGc                                               >  1:8086105/1‑93 (MQ=255)
                                          ggATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATATCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGc                                               <  1:8515591/93‑1 (MQ=255)
                                            aTAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAg                                             <  2:1346537/93‑1 (MQ=255)
                                             tAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAgg                                            <  1:4610519/93‑1 (MQ=255)
                                              aaaaaTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGt                                           <  2:7331770/93‑1 (MQ=255)
                                               aaaaTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGtt                                          <  1:4870798/93‑1 (MQ=255)
                                                aaaTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGttt                                         >  2:2759563/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       >  1:378547/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       <  2:3106477/93‑1 (MQ=255)
                                                   tACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCt                                      >  2:6852281/1‑93 (MQ=255)
                                                     cTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtg                                    >  1:1997100/1‑93 (MQ=255)
                                                     cTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtg                                    <  2:6239362/93‑1 (MQ=255)
                                                     cTGGCTCAACTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtg                                    <  1:3961225/93‑1 (MQ=255)
                                                      tGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtgt                                   >  1:3232989/1‑93 (MQ=255)
                                                        gCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTgg                                 <  2:2514468/93‑1 (MQ=255)
                                                         cTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGt                                <  1:1433261/93‑1 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                               <  2:6571827/93‑1 (MQ=255)
                                                           cACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTc                              <  2:3851342/93‑1 (MQ=255)
                                                            cccTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcg                             <  1:1658751/92‑1 (MQ=255)
                                                            aCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcg                             >  1:2661487/1‑93 (MQ=255)
                                                             ccTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgg                            <  1:721604/93‑1 (MQ=255)
                                                             ccTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgg                            >  2:7758837/1‑93 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           <  1:7861961/93‑1 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           >  1:1380440/1‑93 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           <  1:8436450/93‑1 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           <  1:6095865/93‑1 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           >  2:164986/1‑93 (MQ=255)
                                                               ttACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                          <  1:4795042/93‑1 (MQ=255)
                                                               ttACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                          >  1:3678078/1‑93 (MQ=255)
                                                               ttACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                          >  1:5313747/1‑93 (MQ=255)
                                                                 aCGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacgg                        >  2:4963814/1‑93 (MQ=255)
                                                                     ccaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAg                    <  1:5145319/92‑1 (MQ=255)
                                                                     acaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAg                    <  2:4869565/93‑1 (MQ=255)
                                                                      caATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgc                   >  1:6721279/1‑93 (MQ=255)
                                                                      caATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgc                   >  2:6953496/1‑93 (MQ=255)
                                                                       aaTACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcg                  <  2:108340/93‑1 (MQ=255)
                                                                        aTACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcgg                 >  2:5858758/1‑93 (MQ=255)
                                                                         tACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggg                >  2:3966518/1‑92 (MQ=255)
                                                                           cGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGgcgc              >  2:518282/1‑93 (MQ=255)
                                                                                 aTTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCa        <  1:6530167/93‑1 (MQ=255)
                                                                                   ttCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGc      <  1:3054475/93‑1 (MQ=255)
                                                                                    tCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCa     >  1:7405085/1‑93 (MQ=255)
                                                                                    tCACGCCGCCAGGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGGTGTCTTGTGGTTCGGACGGATAGCGGCGTGCTgcagca     >  2:7794196/1‑93 (MQ=255)
                                                                                     cACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCaa    <  2:2992469/93‑1 (MQ=255)
                                                                                     cACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCaa    >  1:5042890/1‑93 (MQ=255)
                                                                                      aCGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAc   >  1:7336788/1‑93 (MQ=255)
                                                                                      aCGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAc   <  2:4100551/93‑1 (MQ=255)
                                                                                       cgccgcCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAcc  >  2:8270884/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACC  >  NC_000913/1235602‑1235781

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGAAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTT  >  NC_000913/1235599‑1235785
                                                                                             |                                                                                             
CAGAAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCC                                                                                               >  SRR4302541.3924677/1‑94 (MQ=60)
CAGAAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCC                                                                                               <  SRR4302541.5533350/94‑1 (MQ=60)
CAGAAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCC                                                                                               >  SRR4302541.2307865/1‑94 (MQ=60)
CAGAAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCC                                                                                               >  SRR4302541.638274/1‑94 (MQ=60)
   AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCC                                                                                            <  SRR4302541.1934368/94‑1 (MQ=60)
   AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTAGCTCACCTTACGCACAATACGCAGGCTTTCACACCGCC                                                                                            >  SRR4302541.2440932/1‑94 (MQ=60)
     GGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCAT                                                                                          <  SRR4302541.1035976/94‑1 (MQ=60)
     GGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCAT                                                                                          <  SRR4302541.3894362/94‑1 (MQ=60)
     GGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCAT                                                                                          <  SRR4302541.5413481/94‑1 (MQ=60)
     GGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCAT                                                                                          >  SRR4302541.6715133/1‑94 (MQ=60)
      TTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTTCGGATAAAAATACTGGCTCACCTTACTCACAATACGCAGGATTTCACTCCGCCATG                                                                                         >  SRR4302541.7871051/1‑94 (MQ=60)
      GTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATCCGCAGGATTTCACGCCGCCATG                                                                                         <  SRR4302541.3630994/94‑1 (MQ=60)
      GTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATG                                                                                         >  SRR4302541.7871052/1‑94 (MQ=60)
       TCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGA                                                                                        >  SRR4302541.1220206/1‑94 (MQ=60)
        CGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAT                                                                                       <  SRR4302541.8854441/94‑1 (MQ=60)
        CGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAT                                                                                       <  SRR4302541.5525666/94‑1 (MQ=60)
          CCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTC                                                                                     >  SRR4302541.1529177/1‑94 (MQ=60)
          CCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTC                                                                                     >  SRR4302541.5729564/1‑94 (MQ=60)
           CCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCC                                                                                    >  SRR4302541.1554854/1‑94 (MQ=60)
           CCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCC                                                                                    <  SRR4302541.8501476/94‑1 (MQ=60)
           CCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCC                                                                                    >  SRR4302541.6713592/1‑94 (MQ=60)
            CAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCA                                                                                   >  SRR4302541.1395010/1‑94 (MQ=60)
            CAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCA                                                                                   <  SRR4302541.7437474/94‑1 (MQ=60)
            CAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCA                                                                                   <  SRR4302541.2187634/94‑1 (MQ=60)
            CAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCA                                                                                   <  SRR4302541.2187633/94‑1 (MQ=60)
            CAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCA                                                                                   >  SRR4302541.7340712/1‑94 (MQ=60)
            CAGCATTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCA                                                                                   >  SRR4302541.3827447/1‑94 (MQ=60)
             AGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAG                                                                                  >  SRR4302541.4048674/1‑94 (MQ=60)
              GCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGC                                                                                 >  SRR4302541.1689346/1‑94 (MQ=60)
                CTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAG                                                                               >  SRR4302541.1893073/1‑94 (MQ=60)
                CTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAG                                                                               >  SRR4302541.4397899/1‑94 (MQ=60)
                CTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAG                                                                               >  SRR4302541.5977958/1‑94 (MQ=60)
                CTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAG                                                                               >  SRR4302541.5730869/1‑94 (MQ=60)
                  TCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTG                                                                             <  SRR4302541.5603486/94‑1 (MQ=60)
                  TCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTAACGCCGCCATGATTCCAGCAGTG                                                                             <  SRR4302541.1543952/94‑1 (MQ=60)
                   CGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGG                                                                            >  SRR4302541.2179647/1‑94 (MQ=60)
                   CGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGG                                                                            <  SRR4302541.3348292/94‑1 (MQ=60)
                   CGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGATCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGG                                                                            >  SRR4302541.7287690/1‑94 (MQ=60)
                    GTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGC                                                                           <  SRR4302541.8826808/94‑1 (MQ=60)
                     TTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCG                                                                          >  SRR4302541.6822004/1‑94 (MQ=60)
                      TCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGC                                                                         <  SRR4302541.2671465/94‑1 (MQ=60)
                      TCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGC                                                                         >  SRR4302541.196416/1‑94 (MQ=60)
                       CATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCA                                                                        >  SRR4302541.1982114/1‑94 (MQ=60)
                       CATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGACGCCATGATTCCAGCAGTGGCGCA                                                                        >  SRR4302541.6838674/1‑94 (MQ=60)
                        ATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAT                                                                       <  SRR4302541.1688631/94‑1 (MQ=60)
                        ATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAT                                                                       >  SRR4302541.3854635/1‑94 (MQ=60)
                        ATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAT                                                                       >  SRR4302541.1688047/1‑94 (MQ=60)
                        ATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAT                                                                       <  SRR4302541.1775750/94‑1 (MQ=60)
                         TCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATT                                                                      >  SRR4302541.3758241/1‑94 (MQ=60)
                         TCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATT                                                                      >  SRR4302541.3718540/1‑94 (MQ=60)
                          CACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTT                                                                     >  SRR4302541.6218983/1‑94 (MQ=60)
                          CACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTT                                                                     <  SRR4302541.6820266/94‑1 (MQ=60)
                          CACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTT                                                                     <  SRR4302541.1610642/94‑1 (MQ=60)
                           ACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTT                                                                    <  SRR4302541.5269803/94‑1 (MQ=60)
                           ACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTT                                                                    >  SRR4302541.2892872/1‑94 (MQ=60)
                            CCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTT                                                                   <  SRR4302541.5203942/94‑1 (MQ=60)
                            CCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTT                                                                   >  SRR4302541.984747/1‑94 (MQ=60)
                               GAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTCGCGCATTTTTCTC                                                                >  SRR4302541.319029/1‑94 (MQ=60)
                                AGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCC                                                               <  SRR4302541.7190164/94‑1 (MQ=60)
                                AGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCC                                                               <  SRR4302541.300607/94‑1 (MQ=60)
                                AGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCC                                                               >  SRR4302541.5388365/1‑94 (MQ=60)
                                 GTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCA                                                              <  SRR4302541.6676768/94‑1 (MQ=60)
                                  TCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCAT                                                             >  SRR4302541.8402088/1‑94 (MQ=60)
                                  TCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCAT                                                             >  SRR4302541.6630732/1‑94 (MQ=60)
                                   CTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATA                                                            >  SRR4302541.5518041/1‑94 (MQ=60)
                                   CTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATA                                                            >  SRR4302541.5451168/1‑94 (MQ=60)
                                   CTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATA                                                            >  SRR4302541.1394100/1‑94 (MQ=60)
                                    TGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAA                                                           <  SRR4302541.6245370/94‑1 (MQ=60)
                                     GTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAA                                                          >  SRR4302541.1736560/1‑94 (MQ=60)
                                     GTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAA                                                          >  SRR4302541.3274557/1‑94 (MQ=60)
                                      TTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAA                                                         <  SRR4302541.6184869/94‑1 (MQ=60)
                                      TTTTTGCGCATAAAAAGACTGGCTAACATTACGCCCAATACGCAGGATTTCACGCCGCCATGATCCCAGCAGTGGCGCATTTTTCTCCATAAAA                                                         <  SRR4302541.4899117/94‑1 (MQ=60)
                                       TTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAT                                                        <  SRR4302541.2743124/94‑1 (MQ=60)
                                       TTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAT                                                        <  SRR4302541.2492359/94‑1 (MQ=60)
                                       TTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAT                                                        <  SRR4302541.1174517/94‑1 (MQ=60)
                                        TTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATA                                                       >  SRR4302541.5655687/1‑94 (MQ=60)
                                        TTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATA                                                       >  SRR4302541.7818727/1‑94 (MQ=60)
                                          TGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGA                                                     >  SRR4302541.2159065/1‑94 (MQ=60)
                                           GCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAG                                                    <  SRR4302541.5632598/94‑1 (MQ=60)
                                            CGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGC                                                   >  SRR4302541.8174760/1‑94 (MQ=60)
                                             GGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCA                                                  <  SRR4302541.2225973/94‑1 (MQ=60)
                                             GGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATATCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCA                                                  <  SRR4302541.8607253/94‑1 (MQ=60)
                                               ATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGG                                                <  SRR4302541.7289750/94‑1 (MQ=60)
                                                TAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGT                                               <  SRR4302541.6125812/94‑1 (MQ=60)
                                                TAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGT                                               <  SRR4302541.4657346/94‑1 (MQ=60)
                                                TAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGT                                               >  SRR4302541.1224004/1‑94 (MQ=60)
                                                 AAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTT                                              <  SRR4302541.7143198/94‑1 (MQ=60)
                                                  AAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTT                                             <  SRR4302541.4920483/94‑1 (MQ=60)
                                                  AAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTT                                             <  SRR4302541.4920482/94‑1 (MQ=60)
                                                    AATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTC                                           >  SRR4302541.382033/1‑94 (MQ=60)
                                                      TACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTT                                         <  SRR4302541.4489183/94‑1 (MQ=60)
                                                      TACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTT                                         >  SRR4302541.8208638/1‑94 (MQ=60)
                                                       ACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTG                                        >  SRR4302541.2016430/1‑94 (MQ=60)
                                                       ACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTG                                        <  SRR4302541.2897265/94‑1 (MQ=60)
                                                        CTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGT                                       <  SRR4302541.1446265/94‑1 (MQ=60)
                                                        CTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGT                                       >  SRR4302541.3265633/1‑94 (MQ=60)
                                                        CTGGCTCAACTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGT                                       <  SRR4302541.4000570/94‑1 (MQ=60)
                                                         TGGCTCCCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTG                                      <  SRR4302541.1006633/94‑1 (MQ=60)
                                                         TGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTG                                      >  SRR4302541.3020488/1‑94 (MQ=60)
                                                         TGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTG                                      <  SRR4302541.3281690/94‑1 (MQ=60)
                                                          GGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGG                                     >  SRR4302541.733402/1‑94 (MQ=60)
                                                          GGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGG                                     <  SRR4302541.8322486/94‑1 (MQ=60)
                                                           GCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGT                                    <  SRR4302541.1298573/94‑1 (MQ=60)
                                                            CTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTT                                   <  SRR4302541.1447194/94‑1 (MQ=60)
                                                             TCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTC                                  >  SRR4302541.6583054/1‑94 (MQ=60)
                                                              CACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCG                                 >  SRR4302541.2687953/1‑94 (MQ=60)
                                                              CACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCG                                 >  SRR4302541.5021595/1‑94 (MQ=60)
                                                               CCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGG                                <  SRR4302541.1675117/94‑1 (MQ=60)
                                                               ACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGG                                <  SRR4302541.8535709/94‑1 (MQ=60)
                                                                CCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGA                               >  SRR4302541.1394000/1‑94 (MQ=60)
                                                                CCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGA                               <  SRR4302541.728426/94‑1 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              <  SRR4302541.8527585/94‑1 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              >  SRR4302541.7735099/1‑94 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              >  SRR4302541.5368447/1‑94 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              <  SRR4302541.7948317/94‑1 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              <  SRR4302541.6158095/94‑1 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              >  SRR4302541.5132715/1‑94 (MQ=60)
                                                                 CTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGAC                              >  SRR4302541.3715338/1‑94 (MQ=60)
                                                                  TTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACG                             <  SRR4302541.4843960/94‑1 (MQ=60)
                                                                     CGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGAT                          <  SRR4302541.2103772/94‑1 (MQ=60)
                                                                      GCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATA                         >  SRR4302541.8420076/1‑94 (MQ=60)
                                                                       CACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAG                        >  SRR4302541.3708022/1‑94 (MQ=60)
                                                                        CCAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGC                       <  SRR4302541.5198243/94‑1 (MQ=60)
                                                                        ACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGC                       >  SRR4302541.6792999/1‑94 (MQ=60)
                                                                        ACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGC                       >  SRR4302541.6206515/1‑94 (MQ=60)
                                                                         CAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCG                      <  SRR4302541.2249042/94‑1 (MQ=60)
                                                                          AATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGG                     >  SRR4302541.3602669/1‑94 (MQ=60)
                                                                           ATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGC                    <  SRR4302541.7780224/94‑1 (MQ=60)
                                                                           ATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGC                    <  SRR4302541.8101775/94‑1 (MQ=60)
                                                                             TCGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGC                  <  SRR4302541.475066/94‑1 (MQ=60)
                                                                             ACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGC                  >  SRR4302541.1410719/1‑94 (MQ=60)
                                                                             ACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGC                  <  SRR4302541.393302/94‑1 (MQ=60)
                                                                                CAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCT               >  SRR4302541.4551838/1‑94 (MQ=60)
                                                                                 AGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTT              >  SRR4302541.4301654/1‑94 (MQ=60)
                                                                                  GGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTC             >  SRR4302541.1773860/1‑94 (MQ=60)
                                                                                  GGATTCCCCGCCGCCATGATTCGAGCAGTGGCGCATTTTTCTGCATAAAATAGAGCAGTTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTC             <  SRR4302541.3028851/94‑1 (MQ=60)
                                                                                    ATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAG           >  SRR4302541.2978541/1‑94 (MQ=60)
                                                                                    ATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAG           <  SRR4302541.6596826/94‑1 (MQ=60)
                                                                                      TTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCA         <  SRR4302541.3085403/94‑1 (MQ=60)
                                                                                      TTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCA         >  SRR4302541.7484288/1‑94 (MQ=60)
                                                                                       TCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAA        >  SRR4302541.5094782/1‑94 (MQ=60)
                                                                                        CACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAC       >  SRR4302541.7415131/1‑94 (MQ=60)
                                                                                         ACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACC      >  SRR4302541.4946215/1‑94 (MQ=60)
                                                                                          CGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCG     <  SRR4302541.3601253/94‑1 (MQ=60)
                                                                                           GCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGT    >  SRR4302541.663336/1‑94 (MQ=60)
                                                                                           GCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGT    <  SRR4302541.8868500/94‑1 (MQ=60)
                                                                                             CGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTT  <  SRR4302541.1239668/94‑1 (MQ=60)
                                                                                             CGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTT  >  SRR4302541.6553081/1‑94 (MQ=60)
                                                                                             CGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTT  <  SRR4302541.2021594/94‑1 (MQ=60)
                                                                                             CGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTT  <  SRR4302541.1036328/94‑1 (MQ=60)
                                                                                             |                                                                                             
CAGAAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTT  >  NC_000913/1235599‑1235785

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