Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,750,775 C→G 15.0% A31G (GCC→GGC)  PP_1551 → replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,750,7750CG15.0% 163.9 / 21.2 60A31G (GCC→GGC) PP_1551replication protein O
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base G (5/4);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

AAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGT  >  NC_002947/1750630‑1750920
                                                                                                                                                   |                                                                                                                                                 
aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTcccggcccgg                                                                                                                                               >  1:1007077/1‑150 (MQ=255)
     aaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGAt                                                                                                                                          <  2:915053/150‑1 (MQ=255)
        cGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGACGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGaa                                                                                                                                       >  2:1747303/1‑150 (MQ=255)
            tGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                   <  2:1060338/150‑1 (MQ=255)
               cAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgtg                                                                                                                                <  1:72656/150‑1 (MQ=255)
                aGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGc                                                                                                                               <  1:214363/150‑1 (MQ=255)
                 gggCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCt                                                                                                                              <  1:350769/150‑1 (MQ=255)
                 gggCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCt                                                                                                                              <  1:1211524/150‑1 (MQ=255)
                        tttCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                       <  2:1352751/150‑1 (MQ=255)
                        tttCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                       <  2:1323523/150‑1 (MQ=255)
                            tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGGCCGGGAGAtc                                                                                                                                         >  2:1784189/1‑127 (MQ=255)
                            tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGGCCGGGAGAtc                                                                                                                                         <  1:1784189/128‑2 (MQ=255)
                            tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                   >  2:1794415/1‑134 (MQ=255)
                            tctACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                   <  1:1794415/134‑1 (MQ=255)
                                 aGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcc                                                                                                              >  2:1624306/1‑150 (MQ=255)
                                  gCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                             <  2:1808338/150‑1 (MQ=255)
                                     ttAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCTGGAGATGAAGATTGTGCTGTACGTGGCCAAGGccacca                                                                                                          <  1:1728803/150‑1 (MQ=255)
                                     ttAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGccacca                                                                                                          <  2:300747/150‑1 (MQ=255)
                                        gAACGGGTTAAAGCATGACCAACCTCGCCAGAACTAACAAGGCCAGTGGGGTCAGCCGGCTGGACAACCAGCTACACCAATCACCATGC‑‑CTTACGGAACTCCCGGCCCGGGAGATGAAGATTGTGCTGTCCGTGGCCAACGCCACCATCa                                                                                                       >  1:1877100/1‑150 (MQ=255)
                                        gAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacgtaca                                                                                                       <  1:1439717/150‑6 (MQ=255)
                                        gAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacgtaca                                                                                                       <  1:4683/150‑6 (MQ=255)
                                         aaCGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCaa                                                                                                      <  1:1624090/150‑1 (MQ=255)
                                              gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGGTGAAGATTGTGCTGTACGTGGCCAAGGCCACCAACAACTTCg                                                                                                 >  2:1723286/1‑150 (MQ=255)
                                              gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCg                                                                                                 >  2:588925/1‑150 (MQ=255)
                                              gTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCg                                                                                                 >  2:1857340/1‑150 (MQ=255)
                                               ttAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGGCCGGGAGAtcga                                                                                                                                       >  1:1673999/1‑108 (MQ=255)
                                               ttAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGGCCGGGAGAtcga                                                                                                                                       <  2:1673999/111‑4 (MQ=255)
                                                    gCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCgg                                                                                           <  2:1877100/150‑1 (MQ=255)
                                                            aaCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCaa                                                                                   <  2:1712498/150‑1 (MQ=255)
                                                                 cGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCAt                                                                              <  1:957158/150‑1 (MQ=255)
                                                                  gTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgtg                                                                                                                                <  2:1359947/99‑1 (MQ=255)
                                                                  gTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgtg                                                                                                                                >  1:1359947/1‑99 (MQ=255)
                                                                  gTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATc                                                                             >  1:1733755/1‑150 (MQ=255)
                                                                   tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                       >  1:1734574/1‑107 (MQ=255)
                                                                   tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTg                                                                                                                       <  2:1734574/107‑1 (MQ=255)
                                                                   tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGc                                                                                             >  1:1879655/1‑133 (MQ=255)
                                                                   tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGc                                                                                             <  2:1879655/133‑1 (MQ=255)
                                                                     ccACTTGACAAGTCCAGGGGGGTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCAGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCg                                                                          <  2:1865237/150‑1 (MQ=255)
                                                                      cACTTGACAGGTCCACGGGGTTGTCCCTGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCgg                                                                         <  1:1554009/150‑1 (MQ=255)
                                                                              aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGAc                                                                 >  2:1112447/1‑150 (MQ=255)
                                                                                             aCCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                  <  1:842734/118‑1 (MQ=255)
                                                                                             aCCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                  >  2:842734/1‑118 (MQ=255)
                                                                                                           ccAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGa                                    >  2:676326/1‑150 (MQ=255)
                                                                                                              gCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcgc                                                            <  2:920714/123‑1 (MQ=255)
                                                                                                              gCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcgc                                                            >  1:920714/1‑123 (MQ=255)
                                                                                                              gCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGacac                                 >  2:712140/1‑150 (MQ=255)
                                                                                                                t‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCg                               >  1:259560/1‑150 (MQ=255)
                                                                                                                  cATGGATGGCC‑TGCGGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCCGACACCgt                              <  2:1733755/150‑1 (MQ=255)
                                                                                                                      gatggCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTcc                          >  2:14918/1‑150 (MQ=255)
                                                                                                                        tggCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCaa                        <  1:1112447/150‑1 (MQ=255)
                                                                                                                         ggCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAg                       >  1:1200061/1‑150 (MQ=255)
                                                                                                                            c‑tgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCg                    >  2:775811/1‑150 (MQ=255)
                                                                                                                                    cTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCa             <  1:1031302/150‑1 (MQ=255)
                                                                                                                                      cttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                             >  2:1683143/5‑132 (MQ=255)
                                                                                                                                      cttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtg                             <  1:1683143/128‑1 (MQ=255)
                                                                                                                                       ttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGTGGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGcc          >  1:1076574/4‑150 (MQ=255)
                                                                                                                                          aTCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGTCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgc       >  1:1104255/1‑150 (MQ=255)
                                                                                                                                             tcCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGgc                                                                                         <  2:1043612/65‑1 (MQ=255)
                                                                                                                                             tcCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGgc                                                                                         >  1:1043612/1‑65 (MQ=255)
                                                                                                                                             tcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCtgctgc    >  1:1785206/1‑150 (MQ=255)
                                                                                                                                             tcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCtgctgc    >  2:1621514/1‑150 (MQ=255)
                                                                                                                                               ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATGCCGGCTACCGACATCGCGAAAGCCATCTACACTAATCTTTTCACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGt  >  1:1867726/1‑150 (MQ=255)
                                                                                                                                               ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCgg  >  1:374081/1‑149 (MQ=255)
                                                                                                                                               ccggcccggGAGATGAAGATTGTGCTGTACGCGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAGTGCATCCCGGATCACGTCACTGCTAAAGCTATCCACGCTCACCCTGACACCGTGTCCAAGGCGGGTCCCAGCCTGCTGCGt  >  1:611480/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCT‑CATGGATGGCC‑TGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGT  >  NC_002947/1750630‑1750920

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