Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,527 Δ1 bp intergenic (+36/+23) yjgN → / ← valS hypothetical protein/valine tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,5270G.99.7% 1548.1 / ‑3.9 341intergenic (+36/+23)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (186/154);  total (186/155)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 352 total aligned reads displayed.

AGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCA  >  NC_000913/4478437‑4478614
                                                                                          |                                                                                       
 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGtttt                                                                                   <  1:7709245/93‑1 (MQ=255)
   cttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                                                 <  2:3456654/93‑1 (MQ=255)
    ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                                                >  1:2045704/1‑93 (MQ=255)
     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                                               >  2:4034892/1‑93 (MQ=255)
     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                                               >  1:4527125/1‑93 (MQ=255)
      tttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                                                              <  2:1552809/92‑1 (MQ=255)
      cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                                                              <  2:4921245/93‑1 (MQ=255)
       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGGATGCTCCGGCCTTCGTTTTCATCAc                                                                             >  1:5945968/1‑93 (MQ=255)
        tAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACt                                                                            <  2:830646/93‑1 (MQ=255)
          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                                          >  2:2288292/1‑93 (MQ=255)
          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                                          >  2:5738535/1‑93 (MQ=255)
          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                                          >  2:6100880/1‑93 (MQ=255)
          gCCAGTAGTTCAGGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCTTTTTCATCACtgt                                                                          <  2:6879567/93‑1 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         >  2:4800918/1‑93 (MQ=255)
            cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTgtgt                                                                        <  2:4082296/93‑1 (MQ=255)
             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                                       >  1:253456/1‑93 (MQ=255)
             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTCTTCATCACTGTGtt                                                                       >  2:4264482/1‑93 (MQ=255)
             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCAAGCTTCGGGCTTCGTTTTCATCACTGTGtt                                                                       >  1:7195954/1‑93 (MQ=255)
              gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGttt                                                                      <  1:2641793/93‑1 (MQ=255)
               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                                     <  1:4306363/93‑1 (MQ=255)
               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                                     >  2:5456863/1‑93 (MQ=255)
                  ttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAt                                                                  <  2:6500967/93‑1 (MQ=255)
                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                                                 <  1:2939816/93‑1 (MQ=255)
                    taCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTa                                                                <  1:2316353/92‑1 (MQ=255)
                     aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTAc                                                               <  1:1885969/93‑1 (MQ=255)
                      cGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                                                              <  1:1440265/93‑1 (MQ=255)
                      cGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                                                              <  2:3697356/93‑1 (MQ=255)
                       ggtgAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                                             <  1:2806743/90‑1 (MQ=255)
                       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                                             <  2:451254/93‑1 (MQ=255)
                         ggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcg                                                           >  2:1179158/1‑93 (MQ=255)
                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                                          <  1:371754/93‑1 (MQ=255)
                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                                          >  1:3476126/1‑93 (MQ=255)
                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                                          >  1:7540724/1‑93 (MQ=255)
                           aGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGcgcg                                                         >  2:7514708/1‑93 (MQ=255)
                            gTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcgg                                                        >  1:1034462/1‑93 (MQ=255)
                              tATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggcg                                                      >  1:3401876/1‑93 (MQ=255)
                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGAGGCGa                                                     >  2:1769733/1‑93 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  <  2:4594834/93‑1 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  <  1:5491911/93‑1 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  >  1:617494/1‑93 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  <  2:3969143/93‑1 (MQ=255)
                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                                 <  2:2059244/93‑1 (MQ=255)
                                   cttcttTACCATTACTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                                 >  2:325716/1‑93 (MQ=255)
                                    ttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                                                >  2:1282619/1‑93 (MQ=255)
                                    ttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                                                >  1:7491038/1‑93 (MQ=255)
                                    ttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                                                >  2:5109837/1‑93 (MQ=255)
                                     tcttAACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAg                                               >  2:5048891/1‑93 (MQ=255)
                                      cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTTATCACTGTGGTTTGATTTCAGCGCGGCGATAACAGc                                              >  1:3156054/1‑93 (MQ=255)
                                      cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACaac                                              >  2:6765082/1‑91 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             <  1:7075218/93‑1 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             >  2:7647588/1‑93 (MQ=255)
                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                                            <  1:5550696/93‑1 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                                         <  2:7239996/93‑1 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                                         >  1:4185122/1‑93 (MQ=255)
                                           ccATTTCTGTAAGAGAT‑‑AAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGGTTACAGCGCGGCGATAACAGCctgctg                                        >  1:524622/1‑93 (MQ=255)
                                            cATTTCTGTAAGAGATAAAAAAGGCCCGG‑‑GGATGGTCCCGGCCTCGTTTTTATCACTGTGTTTTTATTACAGCGCGGCGGTAAAAGCccgccg                                        >  2:2910479/1‑91 (MQ=255)
                                            cATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCCGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctt                                        >  1:1617235/1‑92 (MQ=255)
                                            cATTTCTGTAAGAGATAAAAAAGGCCGGGGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGGTAAAAGCctgctg                                        >  1:1551171/1‑93 (MQ=255)
                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                                     >  1:7842964/1‑93 (MQ=255)
                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                                     <  2:1105088/93‑1 (MQ=255)
                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                                     <  1:3235574/93‑1 (MQ=255)
                                                tCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa                                    >  2:2638469/1‑93 (MQ=255)
                                                tCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa                                    >  1:6507236/1‑93 (MQ=255)
                                                  tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                                  >  1:6899229/1‑93 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                                 >  2:6195050/1‑93 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                                 >  2:640071/1‑93 (MQ=255)
                                                    tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCa                                <  1:6300883/93‑1 (MQ=255)
                                                     aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg                               >  1:2252970/1‑93 (MQ=255)
                                                      agagATAAAAAAGTCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt                              <  1:6760627/93‑1 (MQ=255)
                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt                             >  2:2566058/1‑93 (MQ=255)
                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt                             >  2:3787218/1‑93 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           <  2:2381940/90‑1 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           >  2:4200453/4‑93 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           <  1:7361174/90‑1 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           >  1:4642155/4‑93 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           <  2:7369266/90‑1 (MQ=255)
                                                           ataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcg                          >  2:2494154/3‑93 (MQ=255)
                                                           ataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcg                          >  2:3539541/3‑93 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         <  1:5672812/92‑1 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         >  1:7416432/2‑93 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         <  1:896168/92‑1 (MQ=255)
                                                              aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt                       >  1:4474811/1‑93 (MQ=255)
                                                              aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt                       >  2:5487422/1‑93 (MQ=255)
                                                               aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTTATCACTGTGTTTTTATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt                      >  2:7050635/1‑93 (MQ=255)
                                                                aaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc                     <  2:4372326/93‑1 (MQ=255)
                                                                  aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcgc                   <  1:2433621/93‑1 (MQ=255)
                                                                        ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              >  1:4063543/4‑93 (MQ=255)
                                                                         gagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGc             <  1:5076680/91‑1 (MQ=255)
                                                                         gagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGc             >  2:2686147/3‑93 (MQ=255)
                                                                           gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCAt           >  1:4074115/1‑93 (MQ=255)
                                                                           gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCAt           <  2:3881053/93‑1 (MQ=255)
                                                                             atgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAg         >  1:4667034/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        <  2:5781571/93‑1 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  1:5433597/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  1:2242831/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  1:1822623/1‑93 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      <  2:6678557/93‑1 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      >  1:4814933/1‑93 (MQ=255)
                                                                                  ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTc    <  1:3556352/93‑1 (MQ=255)
                                                                                    ggCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCa  <  1:1154658/93‑1 (MQ=255)
                                                                                          |                                                                                       
AGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCA  >  NC_000913/4478437‑4478614

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