Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,601,224 A→G 6.1% K220E (AAG→GAG)  PP_2279 → head‑to‑tail joining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,601,2240AG 6.1% 196.3 / 2.3 65K220E (AAG→GAG) PP_2279head‑to‑tail joining protein
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base G (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

GCAACGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAA  >  NC_002947/2601075‑2601370
                                                                                                                                                     |                                                                                                                                                  
gCAACGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCa                                                                                                                                                    >  1:1117523/1‑150 (MQ=255)
 cAACGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCaa                                                                                                                                                   <  1:252082/150‑1 (MQ=255)
   aCGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCTAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATACCCAGATCGTCTCGATCTACACCCACATCACCCTCAAgg                                                                                                                                                 <  1:743430/150‑1 (MQ=255)
   aCGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATCAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAgg                                                                                                                                                 >  1:1022763/1‑150 (MQ=255)
        cTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAAGGTCTCGATCTACACCAACATCACCTTCAAGGGTgag                                                                                                                                            >  2:515331/1‑150 (MQ=255)
          ggACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGAc                                                                                                                                          >  1:741065/1‑150 (MQ=255)
            aCATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCt                                                                                                                                        <  1:800913/150‑1 (MQ=255)
                 gTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAAc                                                                                                                                   >  2:575681/1‑150 (MQ=255)
                         ggaggaggTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCa                                                                                                                           >  1:1128644/1‑150 (MQ=255)
                          gaggaggTCTCCCTGGCGGTACTGCCAGATGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCag                                                                                                                          <  1:588299/150‑1 (MQ=255)
                          gaggaggTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAAc                                                                                                                                   <  1:670237/141‑1 (MQ=255)
                          gaggaggTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAAc                                                                                                                                   >  2:670237/1‑141 (MQ=255)
                             gaggTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCaggag                                                                                                                       <  2:1184962/150‑1 (MQ=255)
                                    cccTGGCGGTACTGCCAGAGGAGGCCCGGTCATCAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAgg                                                                                                                <  2:1022763/150‑1 (MQ=255)
                                      cTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGc                                                                                                              >  1:340349/1‑150 (MQ=255)
                                                tGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTTGAGGGTgatgtgg                                                                                                                                        <  1:290281/114‑6 (MQ=255)
                                                tGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTTGAGGGTgatgtgg                                                                                                                                        >  2:290281/1‑109 (MQ=255)
                                                  ccAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAgg                                                                                                                                                 >  2:28134/1‑103 (MQ=255)
                                                  ccAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAgg                                                                                                                                                 <  1:28134/103‑1 (MQ=255)
                                                  ccAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTg                                                                                                       <  1:344671/145‑1 (MQ=255)
                                                  ccAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTg                                                                                                       >  2:344671/1‑145 (MQ=255)
                                                   cAGAGGAGGCCCGGTCATTAGTCTAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCg                                                                                                 <  1:515331/150‑1 (MQ=255)
                                                     gaggagGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGc                                                                                               <  2:609060/150‑1 (MQ=255)
                                                         aggCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGCGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTccc                                                                                           <  2:741065/150‑1 (MQ=255)
                                                          ggCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAgg                                                                                                                                                 >  2:1034948/1‑95 (MQ=255)
                                                          ggCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAgg                                                                                                                                                 <  1:1034948/95‑1 (MQ=255)
                                                            cccGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAg                                                                                        <  2:984725/150‑1 (MQ=255)
                                                              cGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGc                                                                                      <  1:575681/150‑1 (MQ=255)
                                                                gTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTgg                                                                                                                                      <  2:1203222/100‑1 (MQ=255)
                                                                gTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTgg                                                                                                                                      >  1:1203222/1‑100 (MQ=255)
                                                                     tAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTAcc                                                                               <  2:1117523/150‑1 (MQ=255)
                                                                       gTCGAGCCTGTCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTgg                                                                                                                                      >  2:857374/1‑93 (MQ=255)
                                                                       gTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTgg                                                                                                                                      <  1:857374/93‑1 (MQ=255)
                                                                        tCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGtt                                                                                                     <  2:1225732/125‑1 (MQ=255)
                                                                        tCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGAGTATCAGGAGGTCAAGGGCCAGATTGtt                                                                                                     >  1:1225732/1‑125 (MQ=255)
                                                                          gAGCCTGGCGATGACTCCGGCGACACCCTTCGAGACGATAAGGAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGaaa                                                                          <  2:340349/150‑1 (MQ=255)
                                                                                ggCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCgg                                                                                                >  2:62574/1‑122 (MQ=255)
                                                                                ggCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCgg                                                                                                <  1:62574/122‑1 (MQ=255)
                                                                                 gCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGt                                                                   >  1:511457/1‑150 (MQ=255)
                                                                                     tGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTgcgc                                                               <  2:533172/150‑1 (MQ=255)
                                                                                      gACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGcc                                                              <  2:1128644/150‑1 (MQ=255)
                                                                                         tCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCAc                                                                                    >  2:356033/1‑125 (MQ=255)
                                                                                         tCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCAc                                                                                    <  1:356033/125‑1 (MQ=255)
                                                                                         tCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTgcgc                                                               >  1:1069355/1‑146 (MQ=255)
                                                                                         tCCGGCGACACCCCCCGACACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACTTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTgcgc                                                               <  2:1069355/146‑1 (MQ=255)
                                                                                                   ccccccGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgattc                                                 <  2:479234/150‑1 (MQ=255)
                                                                                                    cccccGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgattcg                                                >  2:816130/1‑150 (MQ=255)
                                                                                                     ccccGAGCCCATAACAAGATCGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgattcga                                               <  1:218109/150‑1 (MQ=255)
                                                                                                       ccGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGcc                                                              <  2:1131811/133‑1 (MQ=255)
                                                                                                       ccGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGcc                                                              >  1:1131811/1‑133 (MQ=255)
                                                                                                               aTAACAAGAACGTCTCGATCACCACCCACCTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCTTGTCGCCCTTGCTAGGCACTATTGCGCCTGGCTCCCGATTCGATTCGTCAAGa                                     <  1:542944/150‑1 (MQ=255)
                                                                                                                    aaGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGc                                                                                      <  2:645795/96‑1 (MQ=255)
                                                                                                                    aaGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGc                                                                                      >  1:645795/1‑96 (MQ=255)
                                                                                                                    aaGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGa                                                                        >  2:630818/1‑110 (MQ=255)
                                                                                                                    aaGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGa                                                                        <  1:630818/110‑1 (MQ=255)
                                                                                                                               gATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTa                     <  1:816130/150‑1 (MQ=255)
                                                                                                                                        aggtctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGcc                >  1:695501/6‑146 (MQ=255)
                                                                                                                                        aggtctCACCCTTGAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGcc                <  2:695501/141‑1 (MQ=255)
                                                                                                                                            aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATgg                  <  2:1082418/140‑1 (MQ=255)
                                                                                                                                            aTCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATgg                  >  1:1082418/1‑140 (MQ=255)
                                                                                                                                                 ccTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGcc                                                              <  2:577734/91‑1 (MQ=255)
                                                                                                                                                 ccTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGcc                                                              >  1:577734/1‑91 (MQ=255)
                                                                                                                                                  cTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGaa  >  2:1176941/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgattcga                                               >  2:125212/1‑104 (MQ=255)
                                                                                                                                                   tCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCcgattcga                                               <  1:125212/104‑1 (MQ=255)
                                                                                                                                                    cAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTAt                    >  2:687335/1‑130 (MQ=255)
                                                                                                                                                    cAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTAt                    <  1:687335/130‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCAACGTCCTGGACATCGTGGTCAAGGAGGAGGTCTCCCTGGCGGTACTGCCAGAGGAGGCCCGGTCATTAGTCGAGCCTGGCGATGACTCCGGCGACACCCCCCGAGACCATAACAAGAACGTCTCGATCTACACCCACATCACCCTCAAGGGTGAGACCTGGAACGTGTATCAGGAGGTCAAGGGCCAGATTGTTCCCGGCTCCCGAGGCACGTACCCGAAAGACAAGTGCGCCTGGCTCCCGATTCGATTCGTCAAGATCGACGGTGAGAACTATGGCCGCTCCTACGTGGAA  >  NC_002947/2601075‑2601370

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