Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,151,910 G→T 6.0% P21H (CCT→CAT)  PP_5696 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,151,9100GT 6.0% 214.7 / 2.8 66P21H (CCT→CAT) PP_5696hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base T (2/2);  total (30/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 = 6.82e-01

CCATGGTTCGGATTGAGCTGTCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACTCGGGC  >  NC_002947/5151767‑5152058
                                                                                                                                               |                                                                                                                                                    
ccATGGTTCGGATTGAGCTGTCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGcct                                                                                                                                                >  2:858910/1‑149 (MQ=255)
        cGGATTGAGCTGTCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCC‑TCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAgg                                                                                                                                       <  1:974195/150‑1 (MQ=255)
         ggATTGAGCTGTCACGCCGAGCCCTGTTCAAGAAGGTATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAgg                                                                                                                                       <  2:586655/150‑1 (MQ=255)
                 cTGTCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCtg                                                                                                                               <  1:45047/150‑1 (MQ=255)
                  tGTCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCtgt                                                                                                                              >  2:450801/1‑150 (MQ=255)
                    tCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGa                                                                                                                            >  1:1008089/1‑150 (MQ=255)
                    tCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGa                                                                                                                            >  2:1160039/1‑150 (MQ=255)
                     cACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGaa                                                                                                                           <  2:717881/150‑1 (MQ=255)
                         ccGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGgt                                                                                                                       <  2:148903/150‑1 (MQ=255)
                         ccGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGgt                                                                                                                       >  1:655228/1‑150 (MQ=255)
                          cGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGgtg                                                                                                                      >  1:1603238/1‑150 (MQ=255)
                                cTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAga                                                                                                                <  1:450801/150‑1 (MQ=255)
                                 tGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAgat                                                                                                               <  2:97431/150‑1 (MQ=255)
                                     cAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATc                                                                                                           <  1:1403530/150‑1 (MQ=255)
                                      aagaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCa                                                                                                          <  1:1096925/150‑1 (MQ=255)
                                      aagaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCa                                                                                                          <  2:586181/150‑1 (MQ=255)
                                      aagaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCa                                                                                                          >  2:649083/1‑150 (MQ=255)
                                      aagaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCa                                                                                                          <  1:392785/150‑1 (MQ=255)
                                       agaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCtgt                                                                                                                              <  1:1501347/129‑1 (MQ=255)
                                       agaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCag                                                                                                         >  1:193386/1‑150 (MQ=255)
                                       agaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGATGCCAGCGCAAAGGGTGTCCtgt                                                                                                                              >  2:1501347/1‑129 (MQ=255)
                                        gaagGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCaga                                                                                                        <  1:398745/150‑1 (MQ=255)
                                          agGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAgaga                                                                                                      <  1:698797/150‑1 (MQ=255)
                                              atgcatCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCa                                                                                                                                                <  1:126691/104‑1 (MQ=255)
                                              atgcatCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCa                                                                                                                                                >  2:126691/1‑104 (MQ=255)
                                                gcatCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGc                                                                                                                                                  >  2:143948/1‑100 (MQ=255)
                                                gcatCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGc                                                                                                                                                  <  1:143948/100‑1 (MQ=255)
                                                  atCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCTAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTtgtgtg                                                                                              <  1:456385/150‑1 (MQ=255)
                                                  atCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTtgtgtg                                                                                              <  2:1505584/150‑1 (MQ=255)
                                                     cgCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGtc                                                                                           <  1:1303038/149‑1 (MQ=255)
                                                             ggTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCag                                                                                   >  2:572332/1‑150 (MQ=255)
                                                                gACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGc                                                                                <  2:193386/150‑1 (MQ=255)
                                                                   aaTGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACATGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTc                                                                             >  1:578938/1‑150 (MQ=255)
                                                                   aaTGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACtt                                                                                                   <  1:1004621/128‑1 (MQ=255)
                                                                   aaTGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCGGTGAAGGGTGTATAGATGATCAGAGACtt                                                                                                   >  2:1004621/1‑128 (MQ=255)
                                                                        ggaggaTCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTg                                                                        >  1:874863/1‑150 (MQ=255)
                                                                         gaggaTCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGa                                                                       <  1:572332/150‑1 (MQ=255)
                                                                           ggaTCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGaaa                                                                     <  2:1603238/150‑1 (MQ=255)
                                                                              tCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCGAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTc                                                                  >  2:945145/1‑150 (MQ=255)
                                                                              tCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTc                                                                  <  2:655228/150‑1 (MQ=255)
                                                                                 aaTAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGtc                                                               >  2:1282866/1‑150 (MQ=255)
                                                                                         aGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACATGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGtcatcctcat                                                       <  2:578938/150‑1 (MQ=255)
                                                                                              aaTGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCAtcctcatcctca                                                  <  1:1160039/150‑1 (MQ=255)
                                                                                                        tcCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAgat                                                                                                               <  1:1195103/79‑1 (MQ=255)
                                                                                                        tcCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAgat                                                                                                               >  2:1195103/1‑79 (MQ=255)
                                                                                                           aCCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCGTCCTCATCCTCAGGTTGctc                                          <  2:692548/145‑1 (MQ=255)
                                                                                                           aCCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTAGAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCGTCCTCATCCTCAGGTTGctc                                          >  1:692548/1‑145 (MQ=255)
                                                                                                            ccGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCtgt                                                                                                                              <  2:377471/60‑1 (MQ=255)
                                                                                                            ccGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCtgt                                                                                                                              >  1:377471/1‑60 (MQ=255)
                                                                                                                   ccAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCAtg                             >  1:620356/1‑150 (MQ=255)
                                                                                                                         cTGGTAAAAACTCTGGATTACATGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTtgtg                                                                                                >  2:544017/1‑77 (MQ=255)
                                                                                                                         cTGGTAAAAACTCTGGATTACATGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTtgtg                                                                                                <  1:544017/77‑1 (MQ=255)
                                                                                                                            gTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATg                    >  2:561404/1‑150 (MQ=255)
                                                                                                                                    tctGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCaga                                                                                                        >  1:198759/1‑58 (MQ=255)
                                                                                                                                    tctGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCaga                                                                                                        <  2:198759/58‑1 (MQ=255)
                                                                                                                                    tctGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCg            >  2:784572/1‑150 (MQ=255)
                                                                                                                                      tGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTc              <  1:910070/146‑1 (MQ=255)
                                                                                                                                      tGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTc              >  2:910070/1‑146 (MQ=255)
                                                                                                                                      tGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGa          >  1:69367/1‑150 (MQ=255)
                                                                                                                                      tGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCCCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGa          <  2:1008089/150‑1 (MQ=255)
                                                                                                                                         aTTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACt       >  2:1195011/1‑150 (MQ=255)
                                                                                                                                          ttACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACTc      <  2:69367/150‑1 (MQ=255)
                                                                                                                                              aGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAGCTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAAATCGGGc  >  1:195501/1‑150 (MQ=255)
                                                                                                                                              aGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACTCgg    <  2:622777/148‑1 (MQ=255)
                                                                                                                                              aGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACTCgg    >  1:622777/1‑148 (MQ=255)
                                                                                                                                              aGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACTCGGGc  <  2:620356/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCATGGTTCGGATTGAGCTGTCACGCCGAGCCCTGTTCAAGAAGGCATGCATCGGCTTGTAGGTGACAATGGGGAGGATCAAATAGCTGAGCAAAATGTCCTTCTCCACCGGCTCCCAGTGCTGGTAAAAACTCTGGATTACAGGTGCCAGCGCAAAGGGTGTCCTGTGAAGGGTGTATAGATGATCAGAGACTTGTGTGGTCATCACCAGAGCCTCCAGTGAAACTCGTCATCCTCATCCTCAGGTTGCTCTCCTGCGAGCATGTGATAAATGCCGTTCCGGAACTCGGGC  >  NC_002947/5151767‑5152058

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