Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,380,583 T→A 5.8% L71Q (CTG→CAG)  tnpT‑I → cointegrate resolution protein T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,380,5830TA 5.8% 217.5 / 3.4 69L71Q (CTG→CAG) tnpT‑Icointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base A (2/2);  total (34/35)
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.38e-01

GGGCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCG  >  NC_002947/3380434‑3380732
                                                                                                                                                     |                                                                                                                                                     
gggCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCt                                                                                                                                                       <  2:58498/150‑1 (MQ=255)
 ggCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTg                                                                                                                                                      >  2:758741/1‑150 (MQ=255)
    gagaACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTc                                                                                                                                                   <  1:635687/150‑1 (MQ=255)
      gaACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCgg                                                                                                                                                 >  2:1264977/1‑150 (MQ=255)
        aCCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGcc                                                                                                                                               <  1:1093172/150‑1 (MQ=255)
          ccAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCaa                                                                                                                                             >  1:423212/1‑150 (MQ=255)
            aaGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACt                                                                                                                                           <  2:254027/150‑1 (MQ=255)
            aaGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACt                                                                                                                                           >  1:1346372/1‑150 (MQ=255)
                     cgcgGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCa                                                                                                                                  <  1:405735/150‑1 (MQ=255)
                       cgGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGc                                                                                                                                <  2:1032411/150‑1 (MQ=255)
                           aCGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTCCCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTaga                                                                                                                            <  2:849344/150‑1 (MQ=255)
                            cGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTagag                                                                                                                           <  1:1264977/150‑1 (MQ=255)
                            cGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTagag                                                                                                                           <  1:596185/150‑1 (MQ=255)
                              ggTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAgagga                                                                                                                         <  1:758741/150‑1 (MQ=255)
                                        ggCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAg                                                                                                               <  2:317774/150‑1 (MQ=255)
                                               ccGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCg                                                                                                        >  2:588535/1‑150 (MQ=255)
                                                          accaccATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAg                                                                                             >  1:904349/1‑150 (MQ=255)
                                                            caccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgcg                                                                                           >  2:134851/1‑150 (MQ=255)
                                                              ccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGcgcg                                                                                         >  2:18231/1‑150 (MQ=255)
                                                                   cACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaacaa                                                                                    <  2:1079959/150‑1 (MQ=255)
                                                                   cACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaacaa                                                                                    <  1:134851/150‑1 (MQ=255)
                                                                    aCCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaacaac                                                                                   >  2:1351302/1‑150 (MQ=255)
                                                                    aCCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaacaac                                                                                   >  1:1065950/1‑150 (MQ=255)
                                                                        cTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCCGGAGCGCGAACAACTACa                                                                               >  2:210028/1‑150 (MQ=255)
                                                                           ccTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcg                                                                            >  2:412720/1‑150 (MQ=255)
                                                                           ccTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAAcg                                                                            >  2:50776/1‑150 (MQ=255)
                                                                            cTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACt                                                                                                                                           <  2:822055/86‑1 (MQ=255)
                                                                            cTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACt                                                                                                                                           >  1:822055/1‑86 (MQ=255)
                                                                               aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGCCAACTgt                                                                                                                                         <  2:1039832/85‑2 (MQ=38)
                                                                               aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGCCAACTgt                                                                                                                                         >  1:1039832/1‑84 (MQ=38)
                                                                               aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGaa                                                                        <  1:50776/150‑1 (MQ=255)
                                                                               aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGaa                                                                        >  1:76623/1‑150 (MQ=255)
                                                                               aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGaa                                                                        >  1:914297/1‑150 (MQ=255)
                                                                                    aCTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCt                                                                   <  1:1351302/150‑1 (MQ=255)
                                                                                     cTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTc                                                                  <  1:412720/150‑1 (MQ=255)
                                                                                      tCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCg                                                                 >  2:640069/1‑150 (MQ=255)
                                                                                      tCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCg                                                                 >  1:242472/1‑150 (MQ=255)
                                                                                       cGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGa                                                                <  1:417500/150‑1 (MQ=255)
                                                                                        gAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGACACCGTGGCCAAGGAGAGCGCCCCACTACAACGCGAACGGCTCGAt                                                               >  2:1016486/1‑150 (MQ=255)
                                                                                          aGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTa                                                             <  2:1065950/150‑1 (MQ=255)
                                                                                              gCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGaa                                                                        <  1:1277309/135‑1 (MQ=255)
                                                                                              gCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGaa                                                                        >  2:1277309/1‑135 (MQ=255)
                                                                                                gtggcgcGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGaa                                                       >  1:1119517/4‑150 (MQ=255)
                                                                                                           gAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaa                                                                                       <  1:1307453/107‑1 (MQ=255)
                                                                                                           gAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaa                                                                                       >  2:1307453/1‑107 (MQ=255)
                                                                                                            aaaTACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTgg                                                                                                     <  2:935376/92‑1 (MQ=255)
                                                                                                            aaaTACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTgg                                                                                                     >  1:935376/1‑92 (MQ=255)
                                                                                                                aCCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGcc                                       >  2:1038126/1‑150 (MQ=255)
                                                                                                                      gcgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGc                                 <  2:76623/150‑1 (MQ=255)
                                                                                                                       cgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGcc                                <  1:1252166/150‑1 (MQ=255)
                                                                                                                              gCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCca                         >  1:1165659/1‑150 (MQ=255)
                                                                                                                              gCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCca                         <  1:640069/150‑1 (MQ=255)
                                                                                                                                tCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGAt                                                                                                                                    >  2:214500/1‑41 (MQ=255)
                                                                                                                                tCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGAt                                                                                                                                    <  1:214500/41‑1 (MQ=255)
                                                                                                                                 cAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCcagca                      <  2:423212/150‑1 (MQ=255)
                                                                                                                                     tgccAGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTAGCAGaa                                                       >  1:193216/4‑113 (MQ=38)
                                                                                                                                     tgccAGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGaa                                                       <  2:193216/110‑1 (MQ=39)
                                                                                                                                       ccAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGGTGGa                >  1:716568/2‑150 (MQ=255)
                                                                                                                                         aGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCa                                                          <  1:1343778/106‑1 (MQ=255)
                                                                                                                                         aGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCa                                                          >  2:1343778/1‑106 (MQ=255)
                                                                                                                                          gTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCt                                                                   >  2:176046/1‑96 (MQ=255)
                                                                                                                                          gTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCt                                                                   <  1:176046/96‑1 (MQ=255)
                                                                                                                                           ttGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGt            <  1:588535/150‑1 (MQ=255)
                                                                                                                                           ttGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGt            >  1:1041982/1‑150 (MQ=255)
                                                                                                                                            tGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTc           <  2:1041982/150‑1 (MQ=255)
                                                                                                                                            tGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGACGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTc           <  2:869519/150‑1 (MQ=255)
                                                                                                                                             ggCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGGACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCa          <  1:1038126/150‑1 (MQ=255)
                                                                                                                                             ggCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCa          <  2:1119517/150‑1 (MQ=255)
                                                                                                                                                cAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAGGCCGAATCCAGCAGTTGGAAAGTCagag       >  1:1202764/1‑150 (MQ=255)
                                                                                                                                                 aaCCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGagc      >  2:483759/1‑150 (MQ=255)
                                                                                                                                                     tGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCg  <  1:1016486/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCGAGAACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCG  >  NC_002947/3380434‑3380732

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