Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,564 A→T 6.0% intergenic (‑179/+49) PP_5038 ← / ← PP_5039 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,5640AT 6.0% 188.8 / 2.6 66intergenic (‑179/+49)PP_5038/PP_5039lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base T (2/2);  total (35/31)
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.52e-01

AGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGCAA  >  NC_002947/5741417‑5741713
                                                                                                                                                   |                                                                                                                                                     
aGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATg                                                                                                                                                     <  1:593530/150‑1 (MQ=255)
 gTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATgg                                                                                                                                                    <  1:914668/150‑1 (MQ=255)
       tGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATACGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCgg                                                                                                                                              <  1:943139/150‑1 (MQ=255)
          aGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCt                                                                                                                                           >  2:670314/1‑150 (MQ=255)
            tGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTg                                                                                                                                         >  1:766787/1‑150 (MQ=255)
             gACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGc                                                                                                                                        <  2:766787/150‑1 (MQ=255)
                       gttaaaaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGc                                                                                                                              >  2:691732/4‑150 (MQ=255)
                       gTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGc                                                                                                                              >  1:264315/1‑150 (MQ=255)
                          aaaaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTg                                                                                                                           <  1:150271/150‑1 (MQ=255)
                          aaaaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAAACTGAAAGCATTACTCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTg                                                                                                                           >  2:304503/1‑150 (MQ=255)
                            aaTCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCCTCCCGCTGCATGcc                                                                                                                         >  2:529775/1‑149 (MQ=255)
                                  tcccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCAcagcag                                                                                                                   >  1:1273881/1‑150 (MQ=255)
                                    ccAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 <  1:670314/150‑1 (MQ=255)
                                        gTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCACCt                                                                                                             >  1:904200/1‑150 (MQ=255)
                                            ccTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGAACAGCAGGCCCCTTggg                                                                                                         >  1:384162/1‑150 (MQ=255)
                                            ccTGCCCCACCCCCTATGTCCCTACACAATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGAATTAATCATGCGCCCAATAAAAAAGGCCCCCCCCAATAGTAGCGGCCTTGCGTTGCGATGCCTCCACAGCAGGCCCCTTagg                                                                                                         >  2:522989/1‑150 (MQ=255)
                                                    aCCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGaa                                                                                                 <  2:264315/150‑1 (MQ=255)
                                                     cccccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAAc                                                                                                >  2:250447/1‑150 (MQ=255)
                                                      ccccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAAcc                                                                                               >  1:587691/1‑150 (MQ=255)
                                                       cccTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCg                                                                                              <  2:587691/150‑1 (MQ=255)
                                                             gTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGAccc                                                                                        >  2:273252/1‑150 (MQ=255)
                                                                       cATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGc                                                                              >  2:20780/1‑150 (MQ=255)
                                                                         ttGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTc                                                                                  <  1:156915/144‑1 (MQ=255)
                                                                         ttGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTc                                                                                  >  2:156915/1‑144 (MQ=255)
                                                                          tGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                           <  1:529775/150‑1 (MQ=255)
                                                                            gTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTg                                                                                                                           <  1:1182050/100‑1 (MQ=255)
                                                                            gTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTg                                                                                                                           >  2:1182050/1‑100 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCa                                                                                >  2:433296/1‑142 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCa                                                                                <  1:433296/142‑1 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                           >  1:235890/1‑147 (MQ=255)
                                                                             tCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaa                                                                           <  2:235890/147‑1 (MQ=255)
                                                                             tCGGATGGCCACTACTGCGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAg                                                                        <  2:1083188/150‑1 (MQ=255)
                                                                                gATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaaaa                                                                         >  1:18363/1‑146 (MQ=255)
                                                                                gATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCaaaaa                                                                         <  2:18363/146‑1 (MQ=255)
                                                                                     ccACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAtctc                                                                >  2:882959/1‑150 (MQ=255)
                                                                                      cACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCTCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATctct                                                               >  1:245622/1‑150 (MQ=255)
                                                                                         tactGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAAGCTCCAGCAAAAAGTCCATCTCTTCg                                                            >  1:903306/1‑150 (MQ=255)
                                                                                              tgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTg                                                         >  1:1182633/1‑148 (MQ=255)
                                                                                              tgCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTg                                                         <  2:1182633/148‑1 (MQ=255)
                                                                                                    tAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCgg                                                 >  2:887/1‑150 (MQ=255)
                                                                                                     aaTCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTcc                                                                     >  2:841073/1‑129 (MQ=255)
                                                                                                     aaTCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTcc                                                                     <  1:841073/129‑1 (MQ=255)
                                                                                                       tCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCg                                              >  1:376394/1‑150 (MQ=255)
                                                                                                          caAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTc                                           >  2:367841/1‑150 (MQ=255)
                                                                                                             aGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCGCGGCCGCTCCAg                                        <  2:384162/150‑1 (MQ=255)
                                                                                                                    atcgagcgcCGAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAgg                                 <  1:1231196/145‑1 (MQ=255)
                                                                                                                        tGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 <  2:182178/66‑1 (MQ=255)
                                                                                                                        tGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGc                                                                                                                 >  1:182178/1‑66 (MQ=255)
                                                                                                                        tGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgc                             <  1:882959/150‑1 (MQ=255)
                                                                                                                         gcgcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcg                            <  1:489507/150‑1 (MQ=255)
                                                                                                                          cgcCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcgg                           >  2:805309/1‑150 (MQ=255)
                                                                                                                             ccAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGa                                                                                                                                 <  2:1092487/45‑1 (MQ=255)
                                                                                                                             ccAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGa                                                                                                                                 >  1:1092487/1‑45 (MQ=255)
                                                                                                                              cAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGACCCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggtgcg                       <  1:522989/150‑1 (MQ=255)
                                                                                                                               aaTAAAAAAGGCCGCCCCCAATGGGAGCTGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggtgcgg                      <  1:887/150‑1 (MQ=255)
                                                                                                                                   aaaaaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAg                  <  1:273252/150‑1 (MQ=255)
                                                                                                                                      aaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggtgcgg                      >  1:447007/1‑143 (MQ=255)
                                                                                                                                      aaGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTtgcggtgcgg                      <  2:447007/143‑1 (MQ=255)
                                                                                                                                       aGGCCGCTCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAt                                                                   <  2:564347/97‑1 (MQ=255)
                                                                                                                                       aGGCCGCTCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCAt                                                                   >  1:564347/1‑97 (MQ=39)
                                                                                                                                          ccgctCCCATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGa           >  1:234529/1‑150 (MQ=255)
                                                                                                                                              tcccATTGGGGGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCg       <  2:234529/149‑1 (MQ=255)
                                                                                                                                              ccccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCg       >  1:536710/1‑150 (MQ=255)
                                                                                                                                               cccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTGGAGCGGCCGAGCACCAgcgaag                             <  1:1007702/127‑6 (MQ=38)
                                                                                                                                                ccAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCggg     <  1:20780/150‑1 (MQ=255)
                                                                                                                                                 cAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGc    >  2:1273276/1‑150 (MQ=255)
                                                                                                                                                   aTGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGCaa  <  2:245622/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTGGCCTGCAGTGACCTTCGACGTGAAAATCCCTCCCAGGTTCCCTGCCCCACCCCCTATGTCCCTACACCATTGGTCGGATGGCCACTACTGTGCGATTAATCACAAAGCATTAATCATGCGCCCAATAAAAAAGGCCGCCCCCAATGGGAGCGGCCTTGCGTCGCGATGCCTGCACAGCAGGCCCCTTGGGCCTCAGAACCGGGACCCAGGCTCCAGCAAAAAGTCCATCTCTTCGCTGGTGCTCGGCCGCTCCAGCACCAGGTTGCGGTGCGGGAAGCGCCCGAAGCGGGCAA  >  NC_002947/5741417‑5741713

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