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.5% 631.1 / ‑4.4 183Q258R (CAG→CGG) ycgBputative sporulation protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (90/92);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 191 total aligned reads displayed.

AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACC  >  NC_000913/1235602‑1235781
                                                                                          |                                                                                         
 aGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACgccgcc                                                                                        >  2:5081828/1‑93 (MQ=255)
 aGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACgccgcc                                                                                        <  2:3695866/93‑1 (MQ=255)
    tCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTGCGCACAATACGCAGGATTTCACGCCGCCATg                                                                                     <  2:1148379/93‑1 (MQ=255)
    tCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATg                                                                                     >  2:1075449/1‑93 (MQ=255)
      gCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAt                                                                                   <  2:2565673/93‑1 (MQ=255)
      gCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAt                                                                                   <  1:2404336/93‑1 (MQ=255)
       cccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAtt                                                                                  <  2:6385050/93‑1 (MQ=255)
       cccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAtt                                                                                  <  2:1061518/93‑1 (MQ=255)
        ccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTc                                                                                 <  2:1771800/93‑1 (MQ=255)
         cAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTcc                                                                                <  1:3552123/93‑1 (MQ=255)
         cAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTcc                                                                                >  1:3369144/1‑93 (MQ=255)
          aGCCTTCGTTCATCACCTGAGTCTGTATTTGCGGATAAAAATACTGGCTCTCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCca                                                                               <  2:5244831/93‑1 (MQ=255)
           gCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcag                                                                              <  2:3557130/93‑1 (MQ=255)
            ccTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagc                                                                             <  1:5248747/93‑1 (MQ=255)
             cTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagca                                                                            <  1:1149865/93‑1 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                           >  2:5197513/1‑93 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                           >  1:3175138/1‑93 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                           >  1:1515462/1‑93 (MQ=255)
               tCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGt                                                                          >  1:5769770/1‑93 (MQ=255)
               tCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGt                                                                          <  2:368604/93‑1 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                        >  2:5330730/1‑93 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                        <  2:4441734/93‑1 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                        >  2:4377106/1‑93 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                        >  1:198455/1‑93 (MQ=255)
                  ttCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgc                                                                       <  1:2448880/93‑1 (MQ=255)
                  ttCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgc                                                                       >  2:3429266/1‑93 (MQ=255)
                    catcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgcgc                                                                     <  1:5093614/93‑1 (MQ=255)
                    catcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgcgc                                                                     <  2:1345265/93‑1 (MQ=255)
                     atcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCa                                                                    >  1:6451739/1‑93 (MQ=255)
                       caCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtt                                                                  <  1:5731560/93‑1 (MQ=255)
                       caCCGGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtt                                                                  <  2:1261462/93‑1 (MQ=255)
                        aCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAttt                                                                 >  1:2169706/1‑93 (MQ=255)
                         ccTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtttt                                                                >  2:1516182/1‑93 (MQ=255)
                         ccTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtttt                                                                >  1:3892617/1‑93 (MQ=255)
                         ccTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtttt                                                                >  2:2396576/1‑93 (MQ=255)
                         ccTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtttt                                                                <  1:4973317/93‑1 (MQ=255)
                           tGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtc                                                              >  1:5375858/1‑93 (MQ=255)
                            gAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtct                                                             >  1:4387947/1‑93 (MQ=255)
                            gAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtct                                                             >  2:5338549/1‑93 (MQ=255)
                             aGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtctc                                                            <  2:1286072/93‑1 (MQ=255)
                             aGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtctc                                                            >  2:881934/1‑93 (MQ=255)
                              gTCTGTTTTTGCTGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTcc                                                           <  1:4931244/93‑1 (MQ=255)
                              gTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTcc                                                           <  2:4698281/93‑1 (MQ=255)
                               tCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCa                                                          <  2:3932967/93‑1 (MQ=255)
                                cTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCAt                                                         <  2:6368788/93‑1 (MQ=255)
                                 tGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATa                                                        >  1:376990/1‑93 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                       >  2:2148973/1‑93 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                       >  1:474229/1‑93 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                       >  2:321248/1‑93 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                       >  1:3148519/1‑93 (MQ=255)
                                    ttttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaaaa                                                     >  1:2138507/1‑93 (MQ=255)
                                     tttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAt                                                    <  1:1357437/93‑1 (MQ=255)
                                        gCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATACAATaga                                                 >  2:2294866/1‑93 (MQ=255)
                                        gCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATaga                                                 <  1:343881/93‑1 (MQ=255)
                                          ggATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGc                                               >  1:596292/1‑93 (MQ=255)
                                               aaaaTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGtt                                          >  1:4265599/1‑93 (MQ=255)
                                                aaaTACTGGCTCACCTTACGCACAATACGCAGGCTTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGttt                                         >  1:5857946/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       <  2:2795305/93‑1 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       >  2:547668/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       >  1:5526073/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       >  1:880821/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                       <  2:1446186/93‑1 (MQ=255)
                                                   tACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCt                                      <  1:2331093/93‑1 (MQ=255)
                                                   tACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCt                                      <  1:5509318/93‑1 (MQ=255)
                                                      tGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtgt                                   >  2:4432773/1‑93 (MQ=255)
                                                        gCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTgg                                 <  2:4217823/93‑1 (MQ=255)
                                                         cTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGt                                <  1:19133/93‑1 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                               <  1:5412740/93‑1 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                               <  1:253087/93‑1 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                               >  2:2926655/1‑93 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                               >  1:3592624/1‑93 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                               <  1:1439696/93‑1 (MQ=255)
                                                           cACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTc                              <  1:3942078/93‑1 (MQ=255)
                                                           cACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTc                              >  1:2979224/1‑93 (MQ=255)
                                                            aCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcg                             >  2:4246925/1‑93 (MQ=255)
                                                            aCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcg                             >  2:3950168/1‑93 (MQ=255)
                                                             ccTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgg                            >  1:1667328/1‑93 (MQ=255)
                                                             ccTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgg                            >  2:2894191/1‑93 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           <  2:4249718/93‑1 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           >  2:4956837/1‑93 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                           >  1:2546473/1‑93 (MQ=255)
                                                               ttACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGTTTcggag                          >  2:123719/1‑92 (MQ=255)
                                                               ttACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                          >  1:1013564/1‑93 (MQ=255)
                                                                 aCGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacgg                        >  1:195672/1‑93 (MQ=255)
                                                                  cGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacgga                       <  2:5305058/93‑1 (MQ=255)
                                                                   gCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGAt                      >  2:3086111/1‑93 (MQ=255)
                                                                    cacaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATa                     <  2:5426723/93‑1 (MQ=255)
                                                                       aaTACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcg                  >  1:4238969/1‑93 (MQ=255)
                                                                       aaTACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcg                  >  1:1245889/1‑93 (MQ=255)
                                                                        aTACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcgg                 <  1:1604108/93‑1 (MQ=255)
                                                                         tACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggc                >  2:2009865/1‑93 (MQ=255)
                                                                         tACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggc                <  1:4791219/93‑1 (MQ=255)
                                                                           cGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGgcgc              >  1:2216966/1‑93 (MQ=255)
                                                                            gCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGgcgcg             >  2:6064209/1‑93 (MQ=255)
                                                                                gATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTc         <  2:591670/93‑1 (MQ=255)
                                                                                gATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTc         >  2:3128778/1‑93 (MQ=255)
                                                                                 aTTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCa        <  1:6273720/93‑1 (MQ=255)
                                                                                   ttCACGCCGCCATGAGTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGc      <  1:2583760/93‑1 (MQ=255)
                                                                                     cACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCaa    >  2:823351/1‑93 (MQ=255)
                                                                                       cgccgcCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAcc  <  2:4584747/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACC  >  NC_000913/1235602‑1235781

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