Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,341,059 G→C 10.0% P88P (CCC→CCG) ‡ PP_2937 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,341,0590GC10.0% 248.6 / 16.2 80P88P (CCC→CCG) ‡PP_2937integrase
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base C (3/5);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TCGTCAGTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCTTTGCTAGG  >  NC_002947/3340913‑3341208
                                                                                                                                                  |                                                                                                                                                     
tcgtcAGTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCgggg                                                                                                                                                    <  2:1868399/150‑1 (MQ=255)
      gTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACCCCGCGTTCATc                                                                                                                                              <  2:637960/150‑1 (MQ=255)
       tCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATcc                                                                                                                                             <  2:1724498/150‑1 (MQ=255)
         ttgggccaattttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCca                                                                                                                                           >  1:1606016/1‑150 (MQ=255)
          tgggccaattttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCcag                                                                                                                                          <  1:957219/150‑1 (MQ=255)
          tgggccaattttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCcag                                                                                                                                          <  1:1228732/150‑1 (MQ=255)
          tgggccaattttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCcag                                                                                                                                          <  1:198792/150‑1 (MQ=255)
               caattttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCGAGTTCTAGGAGGAGCTGATCCATTTCGCGTACGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGGGATGAACGCGGGGTTCATCCcacgcatc                                                                                                                                     >  2:366792/1‑150 (MQ=255)
                 attttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCcaggcatcag                                                                                                                                   <  1:2394560/150‑1 (MQ=255)
                 attttgggccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCcaggcatcag                                                                                                                                   <  1:2394229/150‑1 (MQ=255)
                        gccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAgg                                                                                                                            <  2:1096164/150‑1 (MQ=255)
                        gccaattCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAgg                                                                                                                            >  1:1960644/1‑150 (MQ=255)
                             ttCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGc                                                                                                                       <  1:965870/150‑1 (MQ=255)
                              tCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCg                                                                                                                      <  2:1358064/150‑1 (MQ=255)
                                gACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTa                                                                                                                    <  2:1775690/150‑1 (MQ=255)
                                gACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTa                                                                                                                    <  1:1982268/150‑1 (MQ=255)
                                   tCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGgt                                                                                                                 >  1:233018/1‑150 (MQ=255)
                                   tCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGgt                                                                                                                 >  2:920704/1‑150 (MQ=255)
                                               cTTCTCCAGCTCCGTAGTGTCTGTGGAGTGGCGTATCCATTGCTCTTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTc                                                                                                     <  2:263930/150‑1 (MQ=255)
                                                ttCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGTGTAGGTGTGGCGGGTGTCg                                                                                                    >  1:842484/1‑150 (MQ=255)
                                                ttCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCg                                                                                                    >  2:1460696/1‑150 (MQ=255)
                                                     cAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAgg                                                                                                                            <  1:1917437/121‑1 (MQ=255)
                                                     cAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAgg                                                                                                                            >  2:1917437/1‑121 (MQ=255)
                                                        ctccctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACg                                                                                            >  2:2165036/1‑150 (MQ=255)
                                                          ccctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATcc                                                                                                                                             >  2:379288/1‑99 (MQ=255)
                                                          ccctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATcc                                                                                                                                             <  1:379288/99‑1 (MQ=255)
                                                           cctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGc                                                                                          <  1:1559770/149‑1 (MQ=255)
                                                           cctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGc                                                                                          >  2:1559770/1‑149 (MQ=255)
                                                            ctccAGTCCGAGGAGGAGCTGCTCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGt                                                                                        >  2:437076/1‑150 (MQ=255)
                                                                        gaggagCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttc                                                                            >  1:217625/1‑150 (MQ=255)
                                                                             gCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTgg                                                                                                                        <  2:1885360/101‑1 (MQ=255)
                                                                             gCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTgg                                                                                                                        >  1:1885360/1‑101 (MQ=255)
                                                                              cTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttcagcgcc                                                                      >  1:875296/1‑150 (MQ=255)
                                                                                 aTCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttcagcgccttc                                                                   >  2:858059/1‑150 (MQ=255)
                                                                                    cATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttcagcgccttcagc                                                                <  2:1750918/150‑1 (MQ=255)
                                                                                    cATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttcagcgccttcagc                                                                >  1:1191577/1‑150 (MQ=255)
                                                                                     aTTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCTAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttcagcgccttcagcg                                                               >  1:171090/1‑150 (MQ=255)
                                                                                     aTTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          >  2:1364083/1‑107 (MQ=255)
                                                                                     aTTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          <  1:1364083/107‑1 (MQ=255)
                                                                                          gcgTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGAt                                                          >  2:1500449/1‑150 (MQ=255)
                                                                                            gtaggtagAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGAttt                                                        >  2:973704/1‑150 (MQ=255)
                                                                                                gtagAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          >  2:1011663/1‑96 (MQ=255)
                                                                                                gtagAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          <  1:1011663/96‑1 (MQ=255)
                                                                                                gtagAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAg                                                    >  1:978762/1‑150 (MQ=255)
                                                                                                     gagCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTa                                                                                                  >  2:2154467/1‑99 (MQ=255)
                                                                                                     gagCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTa                                                                                                  <  1:2154467/99‑1 (MQ=255)
                                                                                                       gCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGa                                             <  1:437076/150‑1 (MQ=255)
                                                                                                        cAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGAt                                            >  2:1540650/1‑150 (MQ=255)
                                                                                                         aaCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTc                                                                                                     >  2:599991/1‑92 (MQ=255)
                                                                                                         aaCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTc                                                                                                     <  1:599991/92‑1 (MQ=255)
                                                                                                         aaCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGAtt                                           <  1:1460696/150‑1 (MQ=255)
                                                                                                          aCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTg                                          >  2:2304470/1‑150 (MQ=255)
                                                                                                               tcGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGCAaca                                     >  1:848999/1‑150 (MQ=255)
                                                                                                                 gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTa                                                                                                                    >  2:1122003/1‑69 (MQ=255)
                                                                                                                 gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTa                                                                                                                    <  1:1122003/69‑1 (MQ=255)
                                                                                                                 gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCCCGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACt                                   >  2:491008/1‑150 (MQ=255)
                                                                                                                 gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACt                                   <  1:2053702/150‑1 (MQ=255)
                                                                                                                  aCGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACtt                                  <  2:1960644/150‑1 (MQ=255)
                                                                                                                  aCGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACtt                                  <  2:1646293/150‑1 (MQ=255)
                                                                                                                    gCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTgg                                <  1:973704/150‑1 (MQ=255)
                                                                                                                    gCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTgg                                <  2:875296/150‑1 (MQ=255)
                                                                                                                    gCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTgg                                >  2:856607/1‑150 (MQ=255)
                                                                                                                     cTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGc                               <  2:1606016/150‑1 (MQ=255)
                                                                                                                       cctatgcctgatgcctggGATGAACCCCGCGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCAc                                                                                    <  2:111810‑M1/77‑1 (MQ=255)
                                                                                                                       cctatgcctgatgcctggGATGAACCCCGCGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCAc                                                                                    >  1:111810‑M1/19‑95 (MQ=255)
                                                                                                                           ccAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGAGCCTGCAGCCCCTAGTGTAGGGTGCGGATTGTAACACTTGGCTCGttt                         >  2:1600369/1‑149 (MQ=255)
                                                                                                                            cAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          >  1:611212/1‑68 (MQ=255)
                                                                                                                            cAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          <  2:611212/68‑1 (MQ=255)
                                                                                                                                     tCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGcc               >  1:1842293/1‑150 (MQ=255)
                                                                                                                                     tCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGcc               >  2:1527900/1‑150 (MQ=255)
                                                                                                                                     tCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGAAGACGTGCACGGATGTTCAGCGCCTTCAGAGCCGATTTGAATGGGCGGATTGTAACACGTGGCTCGTTGATCCACAGcc               >  1:1364857/1‑150 (MQ=255)
                                                                                                                                      tggGATGAACCCCGCGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          >  2:898369‑M1/4‑58 (MQ=255)
                                                                                                                                      tggGATGAACCCCGCGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCgg                                                                                                          <  1:898369‑M1/55‑1 (MQ=255)
                                                                                                                                         GATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCATTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCAcc           <  1:451057‑M1/150‑1 (MQ=255)
                                                                                                                                         GATGAACCCCGCGTTCATCCCAGGCATCAGGCATAGGGTGGCGTGGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCAcc           <  2:2013173‑M1/150‑1 (MQ=255)
                                                                                                                                         GATGAACCCCGCGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGAc                                                                                             >  1:1231905‑M1/1‑68 (MQ=255)
                                                                                                                                         GATGAACCCCGCGTTCAGCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGAc                                                                                             <  2:1231905‑M1/68‑1 (MQ=255)
                                                                                                                                            GAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGTCTATGTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCttt        >  2:1191577‑M1/1‑150 (MQ=255)
                                                                                                                                            GAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCttt        <  1:1107429‑M1/150‑1 (MQ=255)
                                                                                                                                            GAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCttt        >  1:858059‑M1/1‑150 (MQ=255)
                                                                                                                                              ACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCTTTGc      >  1:1600369‑M1/1‑150 (MQ=255)
                                                                                                                                                gcgGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGtt                          >  2:1958996/1‑128 (MQ=255)
                                                                                                                                                gcgGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGtt                          <  1:1958996/128‑1 (MQ=255)
                                                                                                                                                  ggggTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGCTCCACAGCCCACCTTTGCTAgg  <  2:217625/150‑1 (MQ=255)
                                                                                                                                                  ggggTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCTTTGCTAgg  <  1:1857959/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCGTCAGTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCTTTGCTAGG  >  NC_002947/3340913‑3341208

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