Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,975,684 A→G 5.4% intergenic (+628/‑516) 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,6840AG 5.4% 216.4 / 2.1 75intergenic (+628/‑516)PP_5727/yifBhypothetical protein/ATP‑dependent chelatase
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base G (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

AAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAGC  >  NC_002947/5975540‑5975829
                                                                                                                                                 |                                                                                                                                                 
aaaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATAtttt                                                                                                                                               >  4:330817/1‑149 (MQ=255)
 aaTAAGAGCTGAGGATTCTCGATAGAGGAAAATAAATTCATCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAAACTGTCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAATATCTTATAGATGTTTTg                                                                                                                                              <  7:248433/149‑1 (MQ=255)
        gCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTaa                                                                                                                                       <  8:300011/149‑1 (MQ=255)
            aggaTTCTCGAGAGAGAAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGc                                                                                                                                   >  2:52509/1‑149 (MQ=255)
              gaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCt                                                                                                                                 >  8:273345/1‑149 (MQ=255)
                ttCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCg                                                                                                                               <  6:142192/149‑1 (MQ=255)
                ttCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCg                                                                                                                               <  8:115395/149‑1 (MQ=255)
                     gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTg                                                                                                                          >  8:348782/1‑149 (MQ=255)
                     gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTg                                                                                                                          >  1:55080/1‑149 (MQ=255)
                     gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTg                                                                                                                          >  6:371106/1‑149 (MQ=255)
                      agagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATTTTTTGCAGTTAAGTGCCTCGGAGTgg                                                                                                                         <  1:208620/149‑1 (MQ=255)
                      agagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTgg                                                                                                                         <  6:64869/149‑1 (MQ=255)
                       gagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGc                                                                                                                        >  2:73313/1‑149 (MQ=255)
                              aaataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTaca                                                                                                                  >  2:187237/1‑148 (MQ=255)
                                ataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAg                                                                                                               <  7:326387/149‑1 (MQ=255)
                                ataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAg                                                                                                               <  5:371106/149‑1 (MQ=255)
                                 taaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAg                                                                                                               >  3:192389/1‑148 (MQ=255)
                                    atTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGaaaa                                                                                                           <  2:142920/149‑1 (MQ=255)
                                     ttCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAg                                                                                                          <  1:73313/149‑1 (MQ=255)
                                          cATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGcc                                                                                                      >  6:335460/1‑148 (MQ=255)
                                              cTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtt                                                                                                 <  1:187237/149‑1 (MQ=255)
                                              cTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtt                                                                                                 >  7:68159/1‑149 (MQ=255)
                                                ccATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtttt                                                                                               >  4:282591/1‑149 (MQ=255)
                                                 cATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtttt                                                                                               >  5:101343/1‑148 (MQ=255)
                                                 cATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtttt                                                                                               <  3:330817/148‑1 (MQ=255)
                                                 cATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAtttt                                                                                               >  1:76889/1‑148 (MQ=255)
                                                  aTCAATGCAAAAGTGAAATTATACATTTTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAAttttta                                                                                             >  8:70929/1‑147 (MQ=255)
                                                  aTCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGa                                                                                             >  8:269408/1‑149 (MQ=255)
                                                  aTCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGa                                                                                             >  7:308130/1‑149 (MQ=255)
                                                   tCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGat                                                                                            >  6:204189/1‑149 (MQ=255)
                                                         cAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcat                                                                                         >  7:22852/1‑146 (MQ=255)
                                                         cAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGAtcatca                                                                                       <  5:204189/148‑1 (MQ=255)
                                                         cAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAg                                                                                      <  3:282591/149‑1 (MQ=255)
                                                          aaaaGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGa                                                                                     <  7:348782/149‑1 (MQ=255)
                                                                   aTTATACATATTGAACCTCGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCAc                                                                            <  4:279812/149‑1 (MQ=255)
                                                                    ttATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCAGACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCAc                                                                            >  1:255566/1‑148 (MQ=255)
                                                                      ataCATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGa                                                                          <  2:55080/148‑1 (MQ=255)
                                                                          atatTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCAAACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGcc                                                                     <  2:255566/149‑1 (MQ=255)
                                                                           tatTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACAATGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGc                                                                      >  1:96176/1‑148 (MQ=255)
                                                                                 aCCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTACAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACtt                                                                                  >  4:266940/1‑129 (MQ=255)
                                                                                 aCCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTACAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACtt                                                                                  <  3:266940/129‑1 (MQ=255)
                                                                                   cTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAACGGGCCGAATTTTGATCTTCAGACTTGTGCGCGAGTGCCATAGctg                                                              >  4:39345/1‑145 (MQ=255)
                                                                                   cTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGaa                                                            <  7:273345/149‑1 (MQ=255)
                                                                                       cAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta                                                        <  3:49511/149‑1 (MQ=255)
                                                                                          aGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagctagct                                                     >  4:300026/1‑149 (MQ=255)
                                                                                                    cGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGa                                            <  6:101343/148‑1 (MQ=255)
                                                                                                           gCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGt                                                                        >  8:309287/1‑113 (MQ=255)
                                                                                                           gCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGt                                                                        <  7:309287/113‑1 (MQ=255)
                                                                                                            cGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgt                                   <  1:52509/149‑1 (MQ=255)
                                                                                                                ccATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTc                               <  5:335460/149‑1 (MQ=255)
                                                                                                                ccATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTc                               <  7:269408/149‑1 (MQ=255)
                                                                                                                 cATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTcc                              >  7:78133/1‑149 (MQ=255)
                                                                                                                 cATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTcc                              >  4:115457/1‑149 (MQ=255)
                                                                                                                 cATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTcc                              >  4:25224/1‑149 (MQ=255)
                                                                                                                      gAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg                                  <  4:120247/140‑1 (MQ=255)
                                                                                                                      gAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg                                  >  3:120247/1‑140 (MQ=255)
                                                                                                                      gAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTAc                         >  7:345801/1‑149 (MQ=255)
                                                                                                                         gAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCg                       >  1:241166/1‑148 (MQ=255)
                                                                                                                            acTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGACTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGaaa                    >  7:29561/2‑149 (MQ=255)
                                                                                                                             cTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAAc                   <  8:308130/149‑1 (MQ=255)
                                                                                                                             cTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAAc                   <  3:236384/149‑1 (MQ=255)
                                                                                                                             cTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGATTGCGACTGATGGGGCTGTGCTCCGGTACCGAAAc                   >  5:177530/1‑149 (MQ=255)
                                                                                                                              tGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACg                  >  1:179966/1‑149 (MQ=255)
                                                                                                                              tGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCTaa                     >  1:138388/1‑146 (MQ=255)
                                                                                                                                cAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGc                 <  2:96176/148‑1 (MQ=255)
                                                                                                                                 aaaaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCAt               <  2:76889/149‑1 (MQ=255)
                                                                                                                                           taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCg                       <  3:251637/131‑1 (MQ=255)
                                                                                                                                           taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCg                       >  4:251637/1‑131 (MQ=255)
                                                                                                                                           taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCa         <  8:22852/145‑1 (MQ=255)
                                                                                                                                           taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt      >  7:229332/1‑148 (MQ=255)
                                                                                                                                           taAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt      <  8:68159/148‑1 (MQ=255)
                                                                                                                                           ataGATGTTTTGCAGTTAAGTGCCTCGGAGTGACTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt      <  8:29561/146‑1 (MQ=255)
                                                                                                                                             aGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGa    <  6:177530/148‑1 (MQ=255)
                                                                                                                                             aGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAg   >  1:187216/1‑149 (MQ=255)
                                                                                                                                              gATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAgc  >  7:94939/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGAC‑CTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAGC  >  NC_002947/5975540‑5975829

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: ‑