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,6920TC84.9% 561.8 / 96.4 245Q258R (CAG→CGG) ycgBputative sporulation protein
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base C (100/108);  total (117/128)
Fisher's exact test for biased strand distribution p-value = 8.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

Only 100 of 247 total aligned reads displayed.

AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAC  >  NC_000913/1235602‑1235780
                                                                                          |                                                                                        
 aGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACgccgcc                                                                                       >  2:4949395/1‑93 (MQ=255)
 aGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACgccgcc                                                                                       >  1:6711426/1‑93 (MQ=255)
   gTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCTCCAt                                                                                     <  2:7369224/93‑1 (MQ=255)
    tCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATg                                                                                    <  1:1094328/93‑1 (MQ=255)
      gCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGAt                                                                                  <  2:7284498/93‑1 (MQ=255)
        ccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTc                                                                                >  1:3343104/1‑93 (MQ=255)
        ccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTc                                                                                <  1:26638/93‑1 (MQ=255)
        ccAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTc                                                                                <  1:725462/93‑1 (MQ=255)
         cAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTcc                                                                               >  2:3993874/1‑93 (MQ=255)
         cAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTcc                                                                               >  1:3030547/1‑93 (MQ=255)
          aGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATTATTCca                                                                              <  2:4721899/93‑1 (MQ=255)
          aGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTccc                                                                              <  1:4871105/93‑2 (MQ=255)
           gCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcag                                                                             >  1:7376521/1‑93 (MQ=255)
           gCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACACCGCCATGATTCcag                                                                             >  2:4842709/1‑93 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                          <  2:1319505/93‑1 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                          >  2:6343947/1‑93 (MQ=255)
              ttCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCcagcag                                                                          >  2:3021566/1‑93 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTGCGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTgg                                                                       <  2:527337/93‑1 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                       <  2:1460733/93‑1 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                       <  1:3786438/93‑1 (MQ=255)
                 gTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTgg                                                                       >  2:3402498/1‑93 (MQ=255)
                  ttCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgc                                                                      >  2:3986065/1‑93 (MQ=255)
                  ttCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGgc                                                                      <  1:4094654/93‑1 (MQ=255)
                      tcaCCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAt                                                                  >  2:1109766/1‑93 (MQ=255)
                         ccTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCAtttt                                                               >  2:2300088/1‑93 (MQ=255)
                           tGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTtc                                                             >  2:5593399/1‑93 (MQ=255)
                           tGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtc                                                             <  2:1889170/93‑1 (MQ=255)
                           tGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtc                                                             <  1:6874513/93‑1 (MQ=255)
                           tGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTtc                                                             <  1:5619181/93‑1 (MQ=255)
                              gTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTcc                                                          <  1:6392474/93‑1 (MQ=255)
                               tCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCa                                                         <  2:6547750/93‑1 (MQ=255)
                               tCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACCCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCa                                                         >  1:4646176/1‑93 (MQ=255)
                                 tGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATa                                                       >  1:6159496/1‑93 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                      <  1:6634475/93‑1 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                      <  1:3185572/93‑1 (MQ=255)
                                  gTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaa                                                      >  1:3883881/1‑93 (MQ=255)
                                   tttttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaaa                                                     >  1:5953535/1‑93 (MQ=255)
                                    ttttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATaaaa                                                    <  1:239128/93‑1 (MQ=255)
                                     tttGCGGATAAAAATACTGGCTCACCTTGCGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAt                                                   >  1:6234479/1‑93 (MQ=255)
                                     tttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAAt                                                   <  2:4451844/93‑1 (MQ=255)
                                      ttGCGGATAAAAATACTGGCTCACCTTGCGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATa                                                  >  2:3299656/1‑93 (MQ=255)
                                      ttGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATa                                                  >  2:6505469/1‑93 (MQ=255)
                                       tGCGGATAAAAATACTGGCTCACCTTGCGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATag                                                 >  1:4493371/1‑93 (MQ=255)
                                       tGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATag                                                 >  2:2768908/1‑93 (MQ=255)
                                         cGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATagag                                               <  2:1401363/93‑1 (MQ=255)
                                         cGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATagag                                               <  1:5181113/93‑1 (MQ=255)
                                          ggATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGc                                              <  1:2015705/93‑1 (MQ=255)
                                           gATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCa                                             <  2:2837519/93‑1 (MQ=255)
                                            aTAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAg                                            <  1:5501968/93‑1 (MQ=255)
                                             tAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAgg                                           >  1:2326012/1‑93 (MQ=255)
                                             tAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAgg                                           >  1:7103395/1‑93 (MQ=255)
                                               aaaaTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGtt                                         >  2:7518408/1‑93 (MQ=255)
                                                 aaTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGtttt                                       <  1:5885925/93‑1 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                      >  1:3438441/1‑93 (MQ=255)
                                                  aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTc                                      <  1:2474104/93‑1 (MQ=255)
                                                   tACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCt                                     >  1:1709439/1‑93 (MQ=255)
                                                   tACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCt                                     >  2:4113934/1‑93 (MQ=255)
                                                    aCTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCtt                                    <  1:883430/93‑1 (MQ=255)
                                                    aCTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCtt                                    <  1:7009170/93‑1 (MQ=255)
                                                    aCTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCtt                                    <  2:6216298/93‑1 (MQ=255)
                                                     cTGGCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtg                                   <  2:2199731/93‑1 (MQ=255)
                                                       ggCTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTtgtg                                 >  1:363481/1‑93 (MQ=255)
                                                        nntCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTgg                                <  1:2006673/91‑1 (MQ=255)
                                                         cTCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGt                               <  1:4809240/93‑1 (MQ=255)
                                                          tCACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGtt                              >  2:5488213/1‑93 (MQ=255)
                                                           cACCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTc                             >  1:2218375/1‑93 (MQ=255)
                                                            aCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcg                            >  1:7339463/1‑93 (MQ=255)
                                                            aCCTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcg                            <  2:3535501/93‑1 (MQ=255)
                                                              cTTACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcgga                          >  2:7011051/1‑93 (MQ=255)
                                                               ttgcGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                         >  2:82104/4‑93 (MQ=255)
                                                               ttgcGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                         >  1:656422/4‑93 (MQ=255)
                                                               ttACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggac                         >  1:1958867/1‑93 (MQ=255)
                                                                tACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacg                        <  1:1929312/93‑1 (MQ=255)
                                                                tACGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacg                        >  1:2532997/1‑93 (MQ=255)
                                                                 aCGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacgg                       <  2:4692498/93‑1 (MQ=255)
                                                                  cGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacgga                      >  2:3465300/1‑93 (MQ=255)
                                                                  cGCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTcggacgga                      <  1:2202383/93‑1 (MQ=255)
                                                                   gCACAATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGAt                     <  2:3732862/93‑1 (MQ=255)
                                                                    cacaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATa                    <  1:1521703/93‑1 (MQ=255)
                                                                     acaATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAg                   >  2:4151081/1‑93 (MQ=255)
                                                                     acaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAg                   >  2:5873344/1‑93 (MQ=255)
                                                                     acaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAg                   >  2:604684/1‑93 (MQ=255)
                                                                     acaATACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAg                   <  1:5429857/93‑1 (MQ=255)
                                                                      caATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgc                  >  1:812086/1‑93 (MQ=255)
                                                                       aaTACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcg                 >  1:6769036/1‑93 (MQ=255)
                                                                         tACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggc               >  2:1756302/1‑93 (MQ=255)
                                                                         tACGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggc               >  2:3144180/1‑93 (MQ=255)
                                                                          aCGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggcg              >  2:7120632/1‑93 (MQ=255)
                                                                          aCGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAgcggcg              >  2:5329471/1‑93 (MQ=255)
                                                                           cGCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGgcgc             <  1:1670156/93‑1 (MQ=255)
                                                                            gCAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGgcgcg            >  2:215127/1‑93 (MQ=255)
                                                                             cAGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGgcgcgc           <  1:1338074/93‑1 (MQ=255)
                                                                              aGGATTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCt          >  1:661301/1‑93 (MQ=255)
                                                                                 aTTTCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCa       <  1:5295570/93‑1 (MQ=255)
                                                                                   ttCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGc     >  1:1237239/1‑93 (MQ=255)
                                                                                    tCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCa    >  2:244176/1‑93 (MQ=255)
                                                                                    tCACGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCa    <  1:1260313/93‑1 (MQ=255)
                                                                                    tCACGCCGCCATGAGTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCa    >  2:5183696/1‑93 (MQ=255)
                                                                                      aCGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAc  <  2:147885/93‑1 (MQ=255)
                                                                                      aCGCCGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAc  <  2:4332964/93‑1 (MQ=255)
                                                                                          |                                                                                        
AAGGTCGCCCAGCCTTCGTTCATCACCTGAGTCTGTTTTTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAAC  >  NC_000913/1235602‑1235780

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