Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,975,678 T→A 6.4% intergenic (+622/‑522) PP_5727 → / → yifB hypothetical protein/ATP‑dependent chelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,975,6780TA 6.4% 181.9 / 3.0 63intergenic (+622/‑522)PP_5727/yifBhypothetical protein/ATP‑dependent chelatase
Reads supporting (aligned to +/- strand):  ref base T (26/33);  new base A (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

TTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAT  >  NC_002947/5975534‑5975825
                                                                                                                                                |                                                                                                                                                    
ttGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAg                                                                                                                                                   <  2:50650/148‑1 (MQ=255)
  gTAAAAATAAGAGCGGAGGATTCTCGCGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTAAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGata                                                                                                                                                >  1:10495/1‑149 (MQ=255)
       aaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg                                                                                                                                           <  8:25197/149‑1 (MQ=255)
        aTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATATATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg                                                                                                                                           >  3:173720/1‑148 (MQ=255)
         tAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCCTGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCa                                                                                                                                         <  1:293234/149‑1 (MQ=255)
          aacagCGGAGGATTCTCCAGACAGTAAAATAAATTCAGCATGCACCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGCAGACCTGCAAAATATCTTATAGATGTTTTGCAg                                                                                                                                        <  8:321882/146‑1 (MQ=255)
               cGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAg                                                                                                                                   >  4:225359/1‑149 (MQ=255)
                 gaggaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTg                                                                                                                                 >  7:182357/1‑149 (MQ=255)
                      ttCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCg                                                                                                                            >  3:54775/1‑149 (MQ=255)
                      ttCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCg                                                                                                                            <  8:182357/149‑1 (MQ=255)
                             gagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGc                                                                                                                     <  1:287562/149‑1 (MQ=255)
                                agGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGctc                                                                                                                   >  3:273720/1‑148 (MQ=255)
                                 ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGACGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGctct                                                                                                                  <  6:118370/148‑1 (MQ=255)
                                 ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGctct                                                                                                                  <  6:57563/148‑1 (MQ=255)
                                 ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa                                                                                                                 <  3:225359/149‑1 (MQ=255)
                                  gAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa                                                                                                                 >  4:249970/1‑148 (MQ=255)
                                   aaaaTAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTaca                                                                                                               <  4:255608/149‑1 (MQ=255)
                                   aaaaTAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTAAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTaca                                                                                                               <  2:10495/149‑1 (MQ=255)
                                    aaataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTacac                                                                                                              <  6:22463/149‑1 (MQ=255)
                                    aaataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTaca                                                                                                               >  1:287483/1‑148 (MQ=255)
                                     aataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATATATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca                                                                                                             <  4:173720/149‑1 (MQ=255)
                                      ataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTaca                                                                                                               >  5:366497/1‑146 (MQ=255)
                                          atTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTCCACAGaa                                                                                                          >  3:225976/1‑147 (MQ=255)
                                             cAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAgg                                                                                                      <  4:82283/148‑1 (MQ=255)
                                              aGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGc                                                                                                    <  1:282893/149‑1 (MQ=255)
                                               gCATTCTACATCAATTCAAAATTGAAATTATACATATTGAACCTGGCAAATATACGATGCTTATAAGCGTTCCATACGAAGACCTGCAAAACATATATAAGATATTTTGCAGTTAAGTGCCTCTGAGTTGCTATACACAGAAAAGGGcc                                                                                                   <  5:59882/149‑1 (MQ=255)
                                                    atCCATCAATTCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtt                                                                                              <  4:125150/148‑1 (MQ=255)
                                                       cATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTg                                                                                           >  4:366756/1‑149 (MQ=255)
                                                         tCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGat                                                                                         <  7:280051/149‑1 (MQ=255)
                                                         tCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGat                                                                                         <  4:54775/149‑1 (MQ=255)
                                                         tCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGa                                                                                          >  5:256868/1‑148 (MQ=255)
                                                             tGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc                                                                                     <  5:222833/149‑1 (MQ=255)
                                                              gCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAggg                                                                                                     >  4:191491/1‑132 (MQ=255)
                                                              gCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAggg                                                                                                     <  3:191491/132‑1 (MQ=255)
                                                               cAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGAtcatca                                                                                    <  1:362679/148‑1 (MQ=255)
                                                               cAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGAtcatca                                                                                    <  6:366497/148‑1 (MQ=255)
                                                                   aGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACtt                                                                               >  5:136589/1‑149 (MQ=255)
                                                                      gAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAg                                                                           >  2:331814/1‑149 (MQ=255)
                                                                                atatTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAggg                                                                                                     <  7:63403/114‑1 (MQ=255)
                                                                                atatTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAggg                                                                                                     >  8:63403/1‑114 (MQ=255)
                                                                                   ttGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATa                                                              <  4:12082/149‑1 (MQ=255)
                                                                                   ttGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATa                                                              <  3:249970/149‑1 (MQ=255)
                                                                                        ccTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACtt                                                                               >  1:266717/1‑128 (MQ=255)
                                                                                        ccTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACtt                                                                               <  2:266717/128‑1 (MQ=255)
                                                                                         cTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGaa                                                        <  2:287483/149‑1 (MQ=255)
                                                                                             cAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta                                                    >  7:94620/1‑149 (MQ=255)
                                                                                               aaGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAatctagc                                                  <  1:331814/149‑1 (MQ=255)
                                                                                                       aTGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACt                                          >  7:245942/1‑149 (MQ=255)
                                                                                                        tGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACt                                          >  6:303744/1‑148 (MQ=255)
                                                                                                              taAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATgggg                                   <  3:366756/149‑1 (MQ=255)
                                                                                                                 gCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtg                                >  4:106471/1‑149 (MQ=255)
                                                                                                                      ccATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTc                           >  5:213738/1‑149 (MQ=255)
                                                                                                                        aTACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCg                         <  6:136589/149‑1 (MQ=255)
                                                                                                                                  cTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAAc               >  4:133117/1‑149 (MQ=255)
                                                                                                                                  cTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAAc               <  6:256868/149‑1 (MQ=255)
                                                                                                                                   tGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACg              >  7:92597/1‑149 (MQ=255)
                                                                                                                                    gCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGAAATAGAAGAAGCTAGCTTGAGACTGATTGGGCTGTGCTCAGGTACCGAAACg              >  3:74780/1‑148 (MQ=255)
                                                                                                                                     cAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCa            <  8:92597/149‑1 (MQ=255)
                                                                                                                                      aaaaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCa            <  3:106471/148‑1 (MQ=255)
                                                                                                                                        aaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGACGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCt                                 <  1:259138/126‑1 (MQ=255)
                                                                                                                                        aaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCt                                 >  2:259138/1‑125 (MQ=255)
                                                                                                                                        aaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCAtt          >  2:242981/1‑148 (MQ=255)
                                                                                                                                           aTCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAAc       >  7:3780/1‑148 (MQ=255)
                                                                                                                                              tataAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAAc       >  6:337491/1‑145 (MQ=255)
                                                                                                                                               ataAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATg                                      >  2:173444/1‑113 (MQ=255)
                                                                                                                                               ataAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATg                                      <  1:173444/113‑1 (MQ=255)
                                                                                                                                                taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAg    >  4:133966/1‑146 (MQ=255)
                                                                                                                                                taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAg    <  3:133966/146‑1 (MQ=255)
                                                                                                                                                taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt  <  6:213738/148‑1 (MQ=255)
                                                                                                                                                taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt  >  3:367083/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAT  >  NC_002947/5975534‑5975825

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: