Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,503 (A)7→6 intergenic (+12/+47) 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,4970A.99.6% 1280.0 / ‑3.1 276intergenic (+6/+53)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (166/109);  total (166/110)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 309 total aligned reads displayed.

CCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTT  >  NC_000913/4478405‑4478587
                                                                                            |                                                                                          
ccttTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATa                                                                                            <  1:3466716/93‑1 (MQ=255)
ccttTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATa                                                                                            <  2:6179544/93‑1 (MQ=255)
 cttTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaa                                                                                           >  1:6259629/1‑93 (MQ=255)
  tttAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaa                                                                                          >  1:4520558/1‑93 (MQ=255)
  tttAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaa                                                                                          >  1:2951520/1‑93 (MQ=255)
  tttAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaa                                                                                          <  1:6351290/93‑1 (MQ=255)
   ttAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaa                                                                                         <  1:2179200/93‑1 (MQ=255)
     aaTCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaa                                                                                       <  2:1759087/93‑1 (MQ=255)
       tCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaagg                                                                                     <  2:3625785/93‑3 (MQ=255)
        cAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaaggc                                                                                    <  2:5790228/93‑4 (MQ=255)
          aTAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCg                                                                                 >  2:3156896/1‑93 (MQ=255)
          aTAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCg                                                                                 >  1:6012203/1‑93 (MQ=255)
          aTAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCg                                                                                 <  1:2917732/93‑1 (MQ=255)
           tAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCgg                                                                                >  2:6592407/1‑93 (MQ=255)
           tAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCgg                                                                                <  2:4667793/93‑1 (MQ=255)
             aagaagAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGag                                                                              >  1:2859954/1‑92 (MQ=255)
             aagaagAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGag                                                                              >  1:5680281/1‑92 (MQ=255)
              agaagaACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagc                                                                             <  2:519518/93‑3 (MQ=255)
              agaagaACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagc                                                                             >  1:5544811/1‑91 (MQ=255)
              agaagaACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagc                                                                             >  2:147313/1‑91 (MQ=255)
               gaagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagca                                                                            >  1:5790288/1‑90 (MQ=255)
               gaagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagca                                                                            <  1:3872148/93‑4 (MQ=255)
                aagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagcat                                                                           >  2:1525722/1‑89 (MQ=255)
                     caacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTct                                                                     <  2:372284/93‑2 (MQ=255)
                        ccAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGc                                                                  >  1:6140981/1‑93 (MQ=255)
                           gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCtt                                                               <  1:432034/93‑1 (MQ=255)
                            aTAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCttc                                                              >  2:1950994/1‑92 (MQ=255)
                             tAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCttcg                                                             <  2:3426785/93‑3 (MQ=255)
                             tAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCttcg                                                             >  1:3118458/1‑91 (MQ=255)
                              aaaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCttcgt                                                            >  2:6682700/1‑90 (MQ=255)
                              aaaCGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCCTTACCATTTCTGTAAGAGATAAAAAAGGCGGGAGCATGCTCTGGCCttcgt                                                            <  1:908171/93‑4 (MQ=255)
                               aaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCttcgtt                                                           >  2:2769275/1‑89 (MQ=255)
                                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGttt                                                         <  1:1620963/93‑1 (MQ=255)
                                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGttt                                                         <  2:1332139/93‑1 (MQ=255)
                                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGttt                                                         <  1:4135522/93‑1 (MQ=255)
                                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGttt                                                         <  2:2353203/93‑1 (MQ=255)
                                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGttt                                                         <  2:5944612/93‑1 (MQ=255)
                                 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGtttt                                                        >  2:3530504/1‑93 (MQ=255)
                                 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGtttt                                                        >  1:2391870/1‑93 (MQ=255)
                                 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGtttt                                                        >  1:4250759/1‑93 (MQ=255)
                                 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGtttt                                                        >  1:2774825/1‑93 (MQ=255)
                                 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGtttt                                                        >  2:243042/1‑93 (MQ=255)
                                  gCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGtttta                                                       >  1:3349895/1‑92 (MQ=255)
                                    ttcttAGACAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGTCGGAGCATGCTCTGGCCTTCGTTTtcat                                                     >  1:4204764/1‑93 (MQ=255)
                                     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTtcatc                                                    <  1:1036076/93‑1 (MQ=255)
                                     tcttAGCCAGTATTTCACGGGGAGTTATGCATTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTtcatc                                                    <  1:2100361/93‑1 (MQ=255)
                                     tcttAGACAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTtcatc                                                    >  1:1502100/1‑93 (MQ=255)
                                      cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCAATTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTtcatca                                                   >  1:5436117/1‑93 (MQ=255)
                                      cttAGCCAGTATTTCACGGCGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTtcatca                                                   <  1:2289376/93‑1 (MQ=255)
                                       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCAc                                                  >  2:3017562/1‑93 (MQ=255)
                                          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACtgt                                               <  1:3301207/93‑1 (MQ=255)
                                            cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTgtgt                                             <  1:5861922/93‑1 (MQ=255)
                                            cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTgtgt                                             <  2:4333406/93‑1 (MQ=255)
                                              gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGttt                                           >  2:3071518/1‑93 (MQ=255)
                                               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGtttt                                          >  2:5874019/1‑93 (MQ=255)
                                               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGtttt                                          >  1:5909407/1‑93 (MQ=255)
                                                aTTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTg                                         >  2:3845466/1‑93 (MQ=255)
                                                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                      <  2:3495922/93‑1 (MQ=255)
                                                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                      <  1:1175338/93‑1 (MQ=255)
                                                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                      <  1:4406790/93‑1 (MQ=255)
                                                       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                  >  2:3461802/1‑93 (MQ=255)
                                                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                 <  2:5587046/93‑1 (MQ=255)
                                                             ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                            >  2:4124992/1‑93 (MQ=255)
                                                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCTTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                          <  1:3380659/93‑1 (MQ=255)
                                                                tGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGAt                         >  2:4557870/1‑93 (MQ=255)
                                                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                        >  1:6697421/1‑93 (MQ=255)
                                                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                      >  2:6102927/1‑93 (MQ=255)
                                                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                      >  2:1085273/1‑93 (MQ=255)
                                                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                      >  2:2243493/1‑93 (MQ=255)
                                                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                      >  2:2555870/1‑93 (MQ=255)
                                                                      cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGc                   <  1:5000507/93‑1 (MQ=255)
                                                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                  >  2:722129/1‑93 (MQ=255)
                                                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                 <  2:1645746/93‑1 (MQ=255)
                                                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                 <  2:4922502/93‑1 (MQ=255)
                                                                         tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctg                >  2:3579016/1‑93 (MQ=255)
                                                                         tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctg                >  1:5025794/1‑93 (MQ=255)
                                                                          aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc               >  1:203943/1‑93 (MQ=255)
                                                                          aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc               >  2:5969790/1‑93 (MQ=255)
                                                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct              >  1:4758776/1‑93 (MQ=255)
                                                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct              >  2:4513868/1‑93 (MQ=255)
                                                                             aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt            <  2:3922841/93‑1 (MQ=255)
                                                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt           >  1:5416540/1‑93 (MQ=255)
                                                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt           >  1:1796753/1‑93 (MQ=255)
                                                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt           >  2:1717050/1‑93 (MQ=255)
                                                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt           <  1:65745/93‑1 (MQ=255)
                                                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATTCTTGGTCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGCCGATAACAGCCTGCTGTTc          <  1:2667203/93‑1 (MQ=255)
                                                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc          >  2:5428363/1‑93 (MQ=255)
                                                                                 cTGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCaa        >  2:4370876/1‑93 (MQ=255)
                                                                                  tGTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt       <  2:4852822/93‑1 (MQ=255)
                                                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc      >  1:2565911/1‑93 (MQ=255)
                                                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc      >  2:405038/1‑93 (MQ=255)
                                                                                    tAAGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCa     >  1:5251813/1‑93 (MQ=255)
                                                                                     aaGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg    <  1:1243049/93‑1 (MQ=255)
                                                                                     aaGAGATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg    >  1:3742701/1‑93 (MQ=255)
                                                                                      agagATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt   <  2:1017086/93‑1 (MQ=255)
                                                                                      agagATAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGGTCAATCAGt   >  2:5059791/1‑93 (MQ=255)
                                                                                      agagATAAAAAAGGCCGGAACATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt   <  2:1855211/93‑1 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  >  1:1970562/1‑93 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  >  2:4890608/1‑93 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  >  2:4235901/1‑93 (MQ=255)
                                                                                            |                                                                                          
CCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTT  >  NC_000913/4478405‑4478587

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: