Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,351,342 A→C 11.7% intergenic (+21/+55) PP_4706 → / ← PP_4707 hypothetical protein/OsmY‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,351,3420AC11.7% 171.5 / 14.2 68intergenic (+21/+55)PP_4706/PP_4707hypothetical protein/OsmY‑like protein
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base C (4/4);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TGCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAA  >  NC_002947/5351198‑5351488
                                                                                                                                                 |                                                                                                                                                  
tgCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACggg                                                                                                                                                >  7:333140/1‑149 (MQ=255)
     cgTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCACAAGAAAAAGCCCTGGCACGGGCCAgg                                                                                                                                           <  1:198607/149‑1 (MQ=255)
       tATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCATCCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCCCGTGCCAGGGc                                                                                                                                         <  5:213564/149‑1 (MQ=255)
         tGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCtt                                                                                                                                       <  3:105287/149‑1 (MQ=255)
             gcgcATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGAt                                                                                                                                   >  4:86268/1‑149 (MQ=255)
                     ccAGACTGCGCATTTTCGCGAGAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGt                                                                                                                             <  8:270103/147‑1 (MQ=255)
                     ccAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCg                                                                                                                                     <  8:65641/139‑1 (MQ=255)
                     ccAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCg                                                                                                                                     >  7:65641/1‑139 (MQ=255)
                     ccAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGt                                                                                                                             >  7:270103/1‑147 (MQ=255)
                     ccAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTgg                                                                                                                           >  3:312076/1‑149 (MQ=255)
                       aGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGAACCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGc                                                                                                                          >  1:255245/1‑148 (MQ=255)
                          cTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCaaa                                                                                                                      <  7:69267/149‑1 (MQ=255)
                             cgcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCACACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCaaa                                                                                                                      <  1:107178/147‑1 (MQ=255)
                             cgcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATc                                                                                                                    >  5:364741/1‑148 (MQ=255)
                                     cgcgACAGCCGGCCATGGGGGGCGCCTCTGAGCCGTGAGCAGTCATGCTTCGAAAGCGCTGTGAGCGCGCGCAACCCGGGGTATTGAGCGCAAGACACCGCCCGGGCCCGGGCCAGGGCGTCGATGAAAGTGGCCAAATCCATTTGGc                                                                                                            <  7:247616/148‑1 (MQ=255)
                                       cgACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGaaa                                                                                                                               >  2:194350/1‑127 (MQ=255)
                                       cgACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTCGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGa                                                                                                         >  5:144346/1‑149 (MQ=255)
                                       cgACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGAGGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGaaa                                                                                                                               <  1:194350/127‑1 (MQ=255)
                                        gACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGAc                                                                                                        >  8:217415/1‑149 (MQ=255)
                                           aGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACgg                                                                                                      >  2:59958/1‑148 (MQ=255)
                                              cGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTg                                                                                                              >  4:122827/1‑137 (MQ=255)
                                              cGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTg                                                                                                              <  3:122827/137‑1 (MQ=255)
                                                tCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAttt                                                                                                               <  6:58413/134‑1 (MQ=255)
                                                tCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAttt                                                                                                               >  5:58413/1‑134 (MQ=255)
                                                  cATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTc                                                                                               <  6:364741/148‑1 (MQ=255)
                                                   aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTcc                                                                                              <  1:226490/148‑1 (MQ=255)
                                                   aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTcc                                                                                              >  2:226490/1‑148 (MQ=255)
                                                   aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCt                                                                                             >  7:274274/1‑149 (MQ=255)
                                                   aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCt                                                                                             >  5:142697/1‑149 (MQ=255)
                                                        tggAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTc                                                                                        <  3:86268/149‑1 (MQ=255)
                                                        tggAGCCTCTGAACCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTc                                                                                        <  2:255245/149‑1 (MQ=255)
                                                             ccTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAGATCCATTTGGCCGACGGGACATTCCTTATTCGGCt                                                                                    <  4:312076/148‑1 (MQ=255)
                                                             ccTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAACCCCTGGCACGGGCCAAGGCTTCGATGAAAATTGCCAAATCCATTTTGCCGAAGGGGCATTCCATATTAGGcat                                                                                   >  3:153271/1‑146 (MQ=255)
                                                               tctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCttt                                                                                  >  2:113592/1‑148 (MQ=255)
                                                                ctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCttt                                                                                  >  1:304633/1‑147 (MQ=255)
                                                                   aTCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAgg                                                                              <  8:333140/148‑1 (MQ=255)
                                                                   aTCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGc                                                                             >  2:234567/1‑149 (MQ=255)
                                                                      cTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTCGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGt                                                                          >  1:183587/1‑149 (MQ=255)
                                                                         gAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCa                                                                        <  2:304633/148‑1 (MQ=255)
                                                                          aGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAg                                                                       <  7:217415/148‑1 (MQ=255)
                                                                           gCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGAc                                                                     <  8:274274/149‑1 (MQ=255)
                                                                                    cTTCGAAAGCTATGTCATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGcc                                                                            <  4:182618/133‑1 (MQ=255)
                                                                                    cTTCGAAAGCGATGTGATCGCGCGAAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGcc                                                                            >  3:182618/1‑133 (MQ=255)
                                                                                           aGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTa                                                                 >  2:361844/1‑137 (MQ=255)
                                                                                           aGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTa                                                                 <  1:361844/137‑1 (MQ=255)
                                                                                           aGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGc                                                     <  1:234567/149‑1 (MQ=255)
                                                                                                     tCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGtt                                           >  6:260960/1‑149 (MQ=255)
                                                                                                                cccGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATg                                 >  3:160823/1‑149 (MQ=255)
                                                                                                                         ttGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCt                        >  6:107484/1‑149 (MQ=255)
                                                                                                                         ttGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCt                        >  1:75108/1‑149 (MQ=255)
                                                                                                                           gAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCtt                       <  5:107484/148‑1 (MQ=255)
                                                                                                                              cgcAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCtnn                   <  2:183587/149‑3 (MQ=255)
                                                                                                                                 aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAAt                                  >  2:323958/1‑131 (MQ=255)
                                                                                                                                 aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAAt                                  <  1:323958/131‑1 (MQ=255)
                                                                                                                                 aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCt                     <  6:133427/144‑1 (MQ=255)
                                                                                                                                 aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCt                     >  5:133427/1‑144 (MQ=255)
                                                                                                                                    tcgaaGCCCTGGCCCGTGCCAGGGCTTTGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCttt                >  4:239028/4‑146 (MQ=255)
                                                                                                                                    tcgaaGCCCTGGCCCGTGCCAGGGCTTTGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCttt                <  3:239028/143‑1 (MQ=255)
                                                                                                                                    aaaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTg                                       <  6:125102/123‑1 (MQ=255)
                                                                                                                                    aaaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTg                                       >  5:125102/1‑123 (MQ=255)
                                                                                                                                     aaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCa                                                                        >  1:260358/1‑89 (MQ=255)
                                                                                                                                     aaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCa                                                                        <  2:260358/89‑1 (MQ=255)
                                                                                                                                       aaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc               <  2:82660/144‑1 (MQ=255)
                                                                                                                                       aaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc               >  1:82660/1‑144 (MQ=255)
                                                                                                                                       aaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGa          >  8:164053/1‑149 (MQ=255)
                                                                                                                                       aaGCCCTCGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGa          <  6:144346/149‑1 (MQ=255)
                                                                                                                                          cccTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTc       >  4:238432/1‑149 (MQ=255)
                                                                                                                                               gCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTaa  >  2:53236/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAA  >  NC_002947/5351198‑5351488

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: