Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,033,546 G→C 5.8% T401S (ACC→AGC) ‡ PP_3556 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,033,5460GC 5.8% 266.2 / 5.5 85T401S (ACC→AGC) ‡PP_3556membrane protein
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base C (2/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGGA  >  NC_002947/4033401‑4033691
                                                                                                                                                 |                                                                                                                                                 
ccGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCgg                                                                                                                                               <  2:1726051/150‑1 (MQ=255)
ccGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCgg                                                                                                                                               <  1:1831823/150‑1 (MQ=255)
  ggCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGc                                                                                                                                             <  1:591019/150‑1 (MQ=255)
   gCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGcc                                                                                                                                            >  2:1558370/1‑150 (MQ=255)
       gggCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAgg                                                                                                                                        <  1:1558370/150‑1 (MQ=255)
       gggCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGACCTCGGGcctgaa                                                                                                                                        <  2:1003329/150‑5 (MQ=255)
            gACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGa                                                                                                                                   <  2:1915272/150‑1 (MQ=255)
            gACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGa                                                                                                                                   <  2:542848/150‑1 (MQ=255)
            gACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGACTTCTACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGa                                                                                                                                   <  1:1550984/150‑1 (MQ=255)
               gCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTc                                                                                                                                <  2:1286987/150‑1 (MQ=255)
                cccTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTcc                                                                                                                               <  2:1266177/150‑1 (MQ=255)
                  cgccccGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCAGCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGa                                                                                                                             <  2:591307/148‑1 (MQ=255)
                  cTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGa                                                                                                                             >  1:430140/1‑150 (MQ=255)
                   tCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAg                                                                                                                            <  2:331659/150‑1 (MQ=255)
                   tCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAg                                                                                                                            <  1:309942/150‑1 (MQ=255)
                   tCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAg                                                                                                                            >  1:1580096/1‑150 (MQ=255)
                         ccTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGc                                                                                                                      >  2:931541/1‑150 (MQ=255)
                                    cATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATc                                                                                                           <  2:1075751/150‑1 (MQ=255)
                                      tCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGc                                                                                                         <  1:1149610/150‑1 (MQ=255)
                                           gCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAAc                                                                                                    <  1:1428463/150‑1 (MQ=255)
                                              aaCCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCa                                                                                                 >  1:921986/1‑150 (MQ=255)
                                              aaCCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCACGTTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCa                                                                                                 >  1:312273/1‑150 (MQ=255)
                                              aaCCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGATGAGAATTCGCTTTCAGGCCCGAGGTCGGGCACCGGGTCGGCTTACACACCAACGAACCGGCATCGCCCAACGCa                                                                                                 >  1:513038/1‑150 (MQ=255)
                                                 cAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCTCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGccc                                                                                                       >  1:94230/1‑141 (MQ=255)
                                                 cAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCTCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGccc                                                                                                       <  2:94230/141‑1 (MQ=255)
                                                 cAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAgg                                                                                                                                        >  1:1715658/1‑108 (MQ=255)
                                                 cAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAgg                                                                                                                                        <  2:1715658/108‑1 (MQ=255)
                                                 cAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCg                                                                                              >  1:1285940/1‑150 (MQ=255)
                                                  aGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCgg                                                                                             >  2:1122653/1‑150 (MQ=255)
                                                   gCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGc                                                                                            >  1:646359/1‑150 (MQ=255)
                                                           ccATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCga                                                                                    <  2:316598/150‑1 (MQ=255)
                                                            cATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGGTGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCgat                                                                                   >  2:15205/1‑150 (MQ=255)
                                                            cATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCa                                                                                                                                          >  1:1025903/1‑95 (MQ=255)
                                                            cATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCa                                                                                                                                          <  2:1025903/95‑1 (MQ=255)
                                                            cATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCgat                                                                                   >  2:1176966/1‑150 (MQ=255)
                                                            cATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCgat                                                                                   >  1:642134/1‑150 (MQ=255)
                                                                   ggCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGcac                                                                            >  1:336523/1‑150 (MQ=255)
                                                                        gTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGt                                                                       >  1:306816/1‑150 (MQ=255)
                                                                        gTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGt                                                                       >  1:306648/1‑150 (MQ=255)
                                                                         ttGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTa                                                                      <  1:414066/150‑1 (MQ=255)
                                                                         ttGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTa                                                                      <  1:954215/150‑1 (MQ=255)
                                                                         ttGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTa                                                                      <  2:642134/150‑1 (MQ=255)
                                                                         ttGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCACGAGGTCGGGCCCAGGGCCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGCACCCATGCTGCACCCCGTa                                                                      >  2:287476/1‑150 (MQ=255)
                                                                               tGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAAc                                                                <  1:1742584/150‑1 (MQ=255)
                                                                                     aaGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCg                                                          >  2:1107055/1‑150 (MQ=255)
                                                                                         tGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCg                                                                                                          <  2:261616/98‑1 (MQ=255)
                                                                                         tGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCg                                                                                                          >  1:261616/1‑98 (MQ=255)
                                                                                           ggCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAg                                                       <  2:1749548/147‑1 (MQ=255)
                                                                                           ggCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAg                                                       >  1:1749548/1‑147 (MQ=255)
                                                                                               ccacGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAATGGGTCGAAGCCCAGtt                                                >  2:975065/1‑150 (MQ=255)
                                                                                               ccacGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGtt                                                <  2:1580096/150‑1 (MQ=255)
                                                                                                cacGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGtt                                                <  1:1567778/149‑1 (MQ=255)
                                                                                                cacGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGtt                                                >  2:1567778/1‑149 (MQ=255)
                                                                                                   ggTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTgcac                                            >  1:44352/1‑150 (MQ=255)
                                                                                                       tCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTgcacgcac                                        >  1:1397268/1‑150 (MQ=255)
                                                                                                         ccAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGc                                                                                                              <  2:75474/78‑1 (MQ=255)
                                                                                                         ccAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGc                                                                                                              >  1:75474/1‑78 (MQ=255)
                                                                                                         ccAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTc                                                           >  1:860099/1‑129 (MQ=255)
                                                                                                         ccAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTc                                                           <  2:860099/129‑1 (MQ=255)
                                                                                                         ccAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAg                                      >  1:942644/1‑150 (MQ=255)
                                                                                                           aggtgtgGTGCTGCTGACCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACAGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTg                                    <  2:974138/147‑1 (MQ=255)
                                                                                                              tgtgGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCa                                 <  2:921986/150‑1 (MQ=255)
                                                                                                                  gtggtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGaa                             <  2:430140/150‑1 (MQ=255)
                                                                                                                  gtggtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGaa                             >  1:652080/1‑150 (MQ=255)
                                                                                                                     gtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGAt                          >  1:824316/1‑150 (MQ=255)
                                                                                                                      tgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATc                         >  2:798058/1‑150 (MQ=255)
                                                                                                                      tgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATc                         >  2:1133013/1‑150 (MQ=255)
                                                                                                                              aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACggg                                                             <  2:1371235/106‑1 (MQ=255)
                                                                                                                              aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACggg                                                             >  1:1371235/1‑106 (MQ=255)
                                                                                                                              aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACgg                                                              >  1:1203569/1‑105 (MQ=255)
                                                                                                                              aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACgg                                                              >  1:1203568/1‑105 (MQ=255)
                                                                                                                              aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACgg                                                              <  2:1203569/105‑1 (MQ=255)
                                                                                                                              aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACgg                                                              <  2:1203568/105‑1 (MQ=255)
                                                                                                                               tCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGAAAAAGAATATCTAAAGTCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAgg                >  2:1910746/1‑150 (MQ=255)
                                                                                                                                ccACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGa               >  1:333062/1‑150 (MQ=255)
                                                                                                                                 cACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGAt              >  2:453717/1‑150 (MQ=255)
                                                                                                                                  aCGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTgcac                                            <  1:1230514/119‑1 (MQ=255)
                                                                                                                                  aCGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTgcac                                            >  2:1230514/1‑119 (MQ=255)
                                                                                                                                   cGTGCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCGAACGCGCCGGCAGTACCGACGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCGCGCACAGTGGCACGAAGATCGCGGCCAGGATcc            <  2:400807/150‑2 (MQ=255)
                                                                                                                                    ccctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTg                                    >  2:722030/6‑125 (MQ=255)
                                                                                                                                    ccctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTg                                    <  1:722030/120‑1 (MQ=255)
                                                                                                                                     cctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGc                    <  1:1583187/136‑1 (MQ=255)
                                                                                                                                     cctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGc                    >  2:1583187/5‑140 (MQ=255)
                                                                                                                                       cAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCAc        <  1:1122653/150‑1 (MQ=255)
                                                                                                                                             cGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGga  <  1:1050513/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGGA  >  NC_002947/4033401‑4033691

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