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.5% 675.0 / NA 209intergenic (+127/‑430)ygdH/sdaChypothetical protein/putative serine transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (115/93);  minor base A (1/0);  total (116/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

CTCGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAATATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATGAT  >  NC_000913/2925731‑2925910
                                                                                          |                                                                                         
ctcGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAacat                                                                                         <  2:904916/93‑1 (MQ=255)
ctcGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAacat                                                                                         <  2:164707/93‑1 (MQ=255)
  cGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGt                                                                                       <  1:2452257/93‑1 (MQ=255)
   gTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTc                                                                                      >  2:1894167/1‑93 (MQ=255)
     ttttGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCAt                                                                                    >  2:1592623/1‑93 (MQ=255)
     ttttGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCAt                                                                                    <  1:1992950/93‑1 (MQ=255)
      tttGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATa                                                                                   <  1:1601094/93‑1 (MQ=255)
             aCTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTa                                                                            >  2:2136541/1‑93 (MQ=255)
             aCTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTa                                                                            >  2:4197372/1‑93 (MQ=255)
              cTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTAt                                                                           >  2:785182/1‑93 (MQ=255)
               ttACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATg                                                                          >  1:4210545/1‑93 (MQ=255)
                tACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGt                                                                         <  2:3233599/93‑1 (MQ=255)
                  cTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTg                                                                       >  1:3620205/1‑93 (MQ=255)
                   tcatcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGc                                                                      <  2:3341917/93‑1 (MQ=255)
                   tcatcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGc                                                                      <  2:3783842/93‑1 (MQ=255)
                   tcatcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGc                                                                      <  1:571620/93‑1 (MQ=255)
                   tcatcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGc                                                                      >  1:4100605/1‑93 (MQ=255)
                    catcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCa                                                                     >  1:2636376/1‑93 (MQ=255)
                     atcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCaa                                                                    >  1:1467777/1‑93 (MQ=255)
                      tcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAAc                                                                   >  2:275926/1‑93 (MQ=255)
                      tcaACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAAc                                                                   >  1:2563399/1‑93 (MQ=255)
                        aaCTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGc                                                                 >  1:12825/1‑93 (MQ=255)
                        aaCTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGc                                                                 >  2:360426/1‑93 (MQ=255)
                         aCTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCa                                                                >  2:456944/1‑93 (MQ=255)
                           tcatttcattTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaaa                                                              >  1:2856184/1‑93 (MQ=255)
                           tcatttcattTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCaaa                                                              >  2:4000365/1‑93 (MQ=255)
                             atttcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACg                                                            >  1:2759180/1‑93 (MQ=255)
                              tttcatttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGt                                                           >  2:4317523/1‑93 (MQ=255)
                                 catttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTc                                                        <  1:1337894/93‑1 (MQ=255)
                                  atttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTcc                                                       >  1:2931774/1‑93 (MQ=255)
                                  atttGTTATATGAATGTTTCTTACCAACCTCACGCGACAAATATCATCACAGTTAACATGGCCTACCATTTATGTTGTAACGCAAACGTTTcc                                                       >  2:803698/1‑93 (MQ=255)
                                   tttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTccc                                                      <  1:4614793/93‑1 (MQ=255)
                                    ttGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCt                                                     <  1:166137/93‑1 (MQ=255)
                                      gTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAt                                                   >  1:3893425/1‑93 (MQ=255)
                                       ttATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAtt                                                  <  2:2551846/93‑1 (MQ=255)
                                        tataTGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTAttt                                                 <  2:1786767/93‑1 (MQ=255)
                                          tatGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTc                                               >  1:251098/1‑93 (MQ=255)
                                           atGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCa                                              >  2:624163/1‑93 (MQ=255)
                                           atGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCa                                              >  2:3532059/1‑93 (MQ=255)
                                             gAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATa                                            >  2:1077873/1‑93 (MQ=255)
                                              aaTGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATaa                                           >  2:1954934/1‑93 (MQ=255)
                                               aTGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATaaa                                          <  1:849796/93‑1 (MQ=255)
                                                 gTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAAcc                                        >  2:3669291/1‑93 (MQ=255)
                                                 gTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAAcc                                        >  1:3242938/1‑93 (MQ=255)
                                                  tttCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGt                                                                         <  1:330003/59‑1 (MQ=255)
                                                   ttCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGt                                      >  1:93287/1‑93 (MQ=255)
                                                    tCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTa                                     >  2:3876264/1‑93 (MQ=255)
                                                    tCTTACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTa                                     >  2:1374795/1‑93 (MQ=255)
                                                      ttACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAtt                                   <  1:1762655/93‑1 (MQ=255)
                                                      ttACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAtt                                   >  2:2533055/1‑93 (MQ=255)
                                                      ttACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAtt                                   >  2:2920740/1‑93 (MQ=255)
                                                      ttACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAtt                                   <  2:2951551/93‑1 (MQ=255)
                                                       tACCACCCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTAttt                                  <  1:128551/93‑1 (MQ=255)
                                                         ccaccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTa                                <  2:4624533/93‑1 (MQ=255)
                                                          caccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTAt                               >  2:2157583/1‑93 (MQ=255)
                                                          caccCTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTAt                               <  1:1007362/93‑1 (MQ=255)
                                                           ccccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATc                              <  1:2230170/92‑1 (MQ=255)
                                                            cgcTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCg                             >  2:4488128/3‑93 (MQ=255)
                                                            cccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCg                             <  2:2081109/93‑1 (MQ=255)
                                                             ccTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGc                            >  1:2103505/1‑93 (MQ=255)
                                                              cTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCt                           >  1:502712/1‑93 (MQ=255)
                                                              cTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCt                           >  2:2589897/1‑93 (MQ=255)
                                                              cTCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTAATAAACCGAATTTTATCGCt                           >  2:2661618/1‑93 (MQ=255)
                                                               tCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTa                          <  1:497706/93‑1 (MQ=255)
                                                               tCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTa                          >  2:217840/1‑93 (MQ=255)
                                                               tCACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTa                          >  2:2313353/1‑93 (MQ=255)
                                                                cACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCtgt                         <  2:2880682/93‑3 (MQ=255)
                                                                cACGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTAt                         >  2:1042643/1‑93 (MQ=255)
                                                                 aCGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATg                        <  1:3661245/93‑1 (MQ=255)
                                                                 aCGCGACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATg                        >  2:31308/1‑93 (MQ=255)
                                                                  cgcgACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGa                       <  2:1975119/93‑1 (MQ=255)
                                                                   gcgACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGaa                      <  1:3391387/93‑1 (MQ=255)
                                                                   gcgACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGaa                      <  2:3005226/93‑1 (MQ=255)
                                                                     gACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGaaaa                    <  2:2703588/93‑1 (MQ=255)
                                                                     gACAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGaaaa                    >  1:3017156/1‑93 (MQ=255)
                                                                      aCAAATATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCGTAAACCGTATTTTATCGCTATGaaaaa                   >  1:1356120/1‑93 (MQ=255)
                                                                        aaaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGa                 <  1:1574404/93‑1 (MQ=255)
                                                                        aaaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGa                 >  1:629282/1‑93 (MQ=255)
                                                                        aaaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGa                 <  1:189028/93‑1 (MQ=255)
                                                                         aaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaa                <  1:4039115/93‑1 (MQ=255)
                                                                         aaTATCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaa                >  1:3477289/1‑93 (MQ=255)
                                                                          atatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaaa               <  2:330060/93‑1 (MQ=255)
                                                                          atatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaaa               <  2:3453420/93‑1 (MQ=255)
                                                                          atatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaaa               <  2:1266018/93‑1 (MQ=255)
                                                                          atatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGaaa               <  2:505550/93‑1 (MQ=255)
                                                                           tatCATCACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAAt              >  2:3847849/1‑93 (MQ=255)
                                                                            atcatcACAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAAtt             <  1:2068667/93‑1 (MQ=255)
                                                                             tcatcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCTCTATGAAAAAGAAATTa            <  1:3647482/93‑1 (MQ=255)
                                                                             tcatcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTa            <  1:3910056/93‑1 (MQ=255)
                                                                              catcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTAt           <  1:3235138/93‑1 (MQ=255)
                                                                               atcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATc          <  2:2439274/93‑1 (MQ=255)
                                                                                tcaCAGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCg         >  2:440655/1‑93 (MQ=255)
                                                                                  acaGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGcc       >  1:1937265/1‑93 (MQ=255)
                                                                                  acaGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGcc       >  1:2545761/1‑93 (MQ=255)
                                                                                    aGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     >  1:2357679/1‑93 (MQ=255)
                                                                                    aGTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     >  2:348883/1‑93 (MQ=255)
                                                                                    aGTTAACATGTCATACAAATTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCAt     >  1:4473873/1‑93 (MQ=255)
                                                                                     gTTAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATg    >  1:1325294/1‑93 (MQ=255)
                                                                                      ttAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATGa   <  2:845405/93‑1 (MQ=255)
                                                                                       tAACATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATGAt  >  1:3102934/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTCGTTTTTGGTCACTTACTCATCAACTCATTTCATTTGTTATATGAATGTTTCTTACCACCCTCACGCGACAAATATCATCACAGTTAATATGTCATACAATTTATGTTGCAACGCAAACGTTTCCCTATTTTCATAAACCGTATTTTATCGCTATGAAAAAGAAATTATCGCCATGAT  >  NC_000913/2925731‑2925910

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