Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,341,057 G→C 9.9% A89G (GCG→GGG)  PP_2937 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,341,0570GC 9.9% 244.9 / 16.7 81A89G (GCG→GGG) PP_2937integrase
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base C (3/5);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

CCCTGTCGTCAGTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCTTTGC  >  NC_002947/3340908‑3341204
                                                                                                                                                     |                                                                                                                                                   
cccTGTCGTCAGTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAAcg                                                                                                                                                     >  1:1358064/1‑150 (MQ=255)
     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                                                                                                             >  2:920704/1‑150 (MQ=255)
                                        tCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGgt                                                                                                             >  1:233018/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                                                                                                                        >  2:1917437/1‑121 (MQ=255)
                                                          cAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAgg                                                                                                                        <  1:1917437/121‑1 (MQ=255)
                                                             ctccctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACg                                                                                        >  2:2165036/1‑150 (MQ=255)
                                                               ccctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATcc                                                                                                                                         >  2:379288/1‑99 (MQ=255)
                                                               ccctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATcc                                                                                                                                         <  1:379288/99‑1 (MQ=255)
                                                                cctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGc                                                                                      >  2:1559770/1‑149 (MQ=255)
                                                                cctccAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGc                                                                                      <  1:1559770/149‑1 (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                                                            >  1:1191577/1‑150 (MQ=255)
                                                                                         cATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGttcagcgccttcagc                                                            <  2:1750918/150‑1 (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                                                                                                 <  1:599991/92‑1 (MQ=255)
                                                                                                              aaCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTc                                                                                                 >  2:599991/1‑92 (MQ=255)
                                                                                                              aaCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGAtt                                       <  1:1460696/150‑1 (MQ=255)
                                                                                                               aCATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTg                                      >  2:2304470/1‑150 (MQ=255)
                                                                                                                    tcGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGCAaca                                 >  1:848999/1‑150 (MQ=255)
                                                                                                                      gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTa                                                                                                                <  1:1122003/69‑1 (MQ=255)
                                                                                                                      gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTa                                                                                                                >  2:1122003/1‑69 (MQ=255)
                                                                                                                      gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCCCGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACt                               >  2:491008/1‑150 (MQ=255)
                                                                                                                      gACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACt                               <  1:2053702/150‑1 (MQ=255)
                                                                                                                       aCGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACtt                              <  2:1646293/150‑1 (MQ=255)
                                                                                                                       aCGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACtt                              <  2:1960644/150‑1 (MQ=255)
                                                                                                                         gCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTgg                            >  2:856607/1‑150 (MQ=255)
                                                                                                                         gCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTgg                            <  1:973704/150‑1 (MQ=255)
                                                                                                                         gCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTgg                            <  2:875296/150‑1 (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                      <  1:1958996/128‑1 (MQ=255)
                                                                                                                                                     gcgGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGtt                      >  2:1958996/1‑128 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCTGTCGTCAGTCTTGGGCCAATTTTGGGCCAATTCGACTCGGGGCGGCAGCTTCTCCAGCTCCCTCCAGTCCGAGGAGGAGCTGATCCATTTCGCGTAGGTAGAGAGCAACATCTCGACGCTGTGGCCAAGCTGGCTCGCGATGAACGCGGGGTTCATCCCAGGCATCAGGCATAGGGTGGCGTAGGTGTGGCGGGTGTCGTACTGACGGCGTTCACGGATGTTCAGCGCCTTCAGCGCCGATTTGAAGTGGCGGATTGTAACACTTGGCTCGTTGATCCACAGCCCACCTTTGC  >  NC_002947/3340908‑3341204

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: