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: