Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,925,821 T→C intergenic (+127/‑430) ygdH → / → sdaC hypothetical protein/putative serine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,925,8210TC99.4% 1113.2 / ‑3.8 349intergenic (+127/‑430)ygdH/sdaChypothetical protein/putative serine transporter
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (186/161);  total (186/163)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 361 total aligned reads displayed.

CTCGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAATATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATGAT  >  NC_000913/2925731‑2925910
                                                                                          |                                                                                         
ctcGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAacat                                                                                         >  1:141777/1‑93 (MQ=255)
ctcGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAacat                                                                                         <  2:4027881/93‑1 (MQ=255)
 tcGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATg                                                                                        >  1:1178047/1‑93 (MQ=255)
   gTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTc                                                                                      >  1:7771111/1‑93 (MQ=255)
    tttttGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCa                                                                                     >  1:3718945/1‑93 (MQ=255)
       ttGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATAc                                                                                  <  1:7458940/93‑1 (MQ=255)
          gTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAAt                                                                               >  1:3049453/1‑93 (MQ=255)
           tCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAAtt                                                                              <  1:974149/93‑1 (MQ=255)
           tCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAAtt                                                                              >  2:1287858/1‑93 (MQ=255)
             aCTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTa                                                                            <  2:265864/93‑1 (MQ=255)
             aCTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTa                                                                            <  1:6820720/93‑1 (MQ=255)
             aCTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTa                                                                            <  1:2768747/93‑1 (MQ=255)
             aCTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTa                                                                            <  2:5699881/93‑1 (MQ=255)
              cTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTAt                                                                           <  2:1963538/93‑1 (MQ=255)
              cTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTAt                                                                           >  2:6282323/1‑93 (MQ=255)
              cTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTAt                                                                           <  1:4337120/93‑1 (MQ=255)
              cTTACTCATCAACTCATTTCATGTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTAt                                                                           <  1:4942857/93‑1 (MQ=255)
                  cTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTg                                                                       >  2:3277716/1‑93 (MQ=255)
                    catcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCa                                                                     <  2:2113478/93‑1 (MQ=255)
                    catcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCa                                                                     >  2:1863540/1‑93 (MQ=255)
                    catcaACTCATTTCATTTGTTATATGAATGCTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCa                                                                     <  2:7730167/93‑1 (MQ=255)
                     atcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCaa                                                                    <  1:4246469/93‑1 (MQ=255)
                      tcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAAc                                                                   <  1:7129280/93‑1 (MQ=255)
                      tcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAAc                                                                   <  2:5385931/93‑1 (MQ=255)
                      tcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAAc                                                                   >  1:7826755/1‑93 (MQ=255)
                        aaCTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGc                                                                 <  2:7511689/93‑1 (MQ=255)
                         aCTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCa                                                                >  2:5206110/1‑93 (MQ=255)
                         aCTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCa                                                                <  1:1270190/93‑1 (MQ=255)
                          cTCATTTCATTTGTTATATGAATGTTTCTTCCCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaa                                                               <  1:2047601/93‑1 (MQ=255)
                          cTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaa                                                               >  2:5259547/1‑93 (MQ=255)
                           tcatttcattTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaaa                                                              <  2:1648624/93‑1 (MQ=255)
                           tcatttcattTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaaa                                                              <  1:824552/93‑1 (MQ=255)
                           tcatttcattTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaaa                                                              <  1:6186187/93‑1 (MQ=255)
                            catttcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAAc                                                             <  1:178799/93‑1 (MQ=255)
                            catttcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAAc                                                             <  1:2324590/93‑1 (MQ=255)
                            catttcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAAc                                                             <  1:7550792/93‑1 (MQ=255)
                              tttcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGt                                                           >  1:5021576/1‑93 (MQ=255)
                                tcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGttt                                                         >  1:6591672/1‑93 (MQ=255)
                                tcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGttt                                                         >  1:6730832/1‑93 (MQ=255)
                                 catttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTc                                                        <  1:5115651/93‑1 (MQ=255)
                                  atttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTcc                                                       <  1:1563873/93‑1 (MQ=255)
                                   tttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTccc                                                      <  1:5731926/93‑1 (MQ=255)
                                   tttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTccc                                                      >  1:4506755/1‑93 (MQ=255)
                                   tttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTccc                                                      <  1:6322794/93‑1 (MQ=255)
                                   tttGTTAAATGAAGTTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTccc                                                      <  2:5824067/93‑1 (MQ=255)
                                      gTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAt                                                   <  2:2286213/93‑1 (MQ=255)
                                      gTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAt                                                   >  1:4686265/1‑93 (MQ=255)
                                       ttATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAtt                                                  <  1:6910586/93‑1 (MQ=255)
                                       ttATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAtt                                                  >  2:4467574/1‑93 (MQ=255)
                                        tataTGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAttt                                                 <  2:291406/93‑1 (MQ=255)
                                         atatGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAtttt                                                >  2:2742282/1‑93 (MQ=255)
                                          tatGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTc                                               >  1:5984060/1‑93 (MQ=255)
                                           atGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCa                                              >  2:6857065/1‑93 (MQ=255)
                                            tGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCAt                                             <  2:2421423/93‑1 (MQ=255)
                                              aaTGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATaa                                           <  1:5170812/93‑1 (MQ=255)
                                               aTGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATaaa                                          <  1:2449012/93‑1 (MQ=255)
                                                 gTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAAcc                                        >  2:4206362/1‑93 (MQ=255)
                                                 gTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAAcc                                        >  2:6388049/1‑93 (MQ=255)
                                                 gTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAAcc                                        >  2:6129227/1‑93 (MQ=255)
                                                 gTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAAcc                                        <  2:4403215/93‑1 (MQ=255)
                                                  tttCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCg                                       >  1:1020068/1‑93 (MQ=255)
                                                     cTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAt                                    >  1:6083036/1‑93 (MQ=255)
                                                      ttACCACCCTCACGCGACAATTATCATCACAGTTAACATGTCATACAATTTATGTTACAACGCAAACGTTTCCCTATTTTCATAAACCGTAtt                                   >  2:2217249/1‑93 (MQ=255)
                                                      ttACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAtt                                   >  2:7182465/1‑93 (MQ=255)
                                                       tACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAttt                                  <  2:7852178/93‑1 (MQ=255)
                                                        accaccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAtttt                                 <  1:2730946/93‑1 (MQ=255)
                                                         ccaccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTa                                <  2:6591386/93‑1 (MQ=255)
                                                         ccaccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTa                                <  2:1819937/93‑1 (MQ=255)
                                                          caccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTAt                               >  1:3786329/1‑93 (MQ=255)
                                                           ccccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATc                              <  2:6726669/92‑1 (MQ=255)
                                                           ccccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATc                              <  2:471274/92‑1 (MQ=255)
                                                           accCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATc                              >  1:4095382/1‑93 (MQ=255)
                                                            cccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCg                             >  2:3537172/1‑93 (MQ=255)
                                                            cccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCg                             >  1:6468844/1‑93 (MQ=255)
                                                               tCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTa                          <  2:5483287/93‑1 (MQ=255)
                                                                cACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTAt                         <  2:5927678/93‑1 (MQ=255)
                                                                cACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTAt                         >  2:2605951/1‑93 (MQ=255)
                                                                    cgACAAATATCATCACAGTTAACATGTCATACAGTTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGaaa                     <  1:7630007/93‑1 (MQ=255)
                                                                       cAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAg                  <  1:6607383/93‑1 (MQ=255)
                                                                        aaaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGa                 >  1:3891651/1‑93 (MQ=255)
                                                                        aaaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGa                 >  2:1062290/1‑93 (MQ=255)
                                                                        aaaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGa                 >  2:1977837/1‑93 (MQ=255)
                                                                         aaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaa                >  2:3815169/1‑93 (MQ=255)
                                                                          atatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaaa               <  1:2908031/93‑1 (MQ=255)
                                                                           tatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAAt              <  2:4528115/93‑1 (MQ=255)
                                                                           tatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAAt              <  1:5388770/93‑1 (MQ=255)
                                                                            atcatcACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAAtt             <  1:7341622/93‑1 (MQ=255)
                                                                             tcatcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTa            <  2:5160128/93‑1 (MQ=255)
                                                                             tcatcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTa            <  2:5020550/93‑1 (MQ=255)
                                                                             tcatcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTa            >  2:1203621/1‑93 (MQ=255)
                                                                             tcatcaCAGTTAACATGTCATACAATTTATGATGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTa            >  1:3510688/1‑93 (MQ=255)
                                                                               atcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATc          >  1:2547967/1‑93 (MQ=255)
                                                                               atcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATc          <  2:4345777/93‑1 (MQ=255)
                                                                               atcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATc          <  2:4054125/93‑1 (MQ=255)
                                                                                 cacaGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGc        >  1:7227664/1‑93 (MQ=255)
                                                                                   caGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCa      >  1:5669631/1‑93 (MQ=255)
                                                                                    aGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     <  2:4836540/93‑1 (MQ=255)
                                                                                    aGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     <  1:3257752/93‑1 (MQ=255)
                                                                                    aGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     >  1:3579236/1‑93 (MQ=255)
                                                                                    aGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     <  2:828213/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTCGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAATATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATGAT  >  NC_000913/2925731‑2925910

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