Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,310,453 G→C 6.2% H362Q (CAC→CAG)  PP_1144 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,310,4530GC 6.2% 229.5 / 4.8 75H362Q (CAC→CAG) PP_1144membrane protein
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base C (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAC  >  NC_002947/1310307‑1310602
                                                                                                                                                  |                                                                                                                                                      
gtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCCGGCTGGTGGCCACGAGCCACGGTGCGCCACCAAAGCGCGCAACCAGGTCGT‑gg                                                                                                                                                    >  1:1373889/1‑150 (MQ=255)
gtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCAGGCAGAACTCGGCGCCACCAAAGCGCGCAACCAGGTCGT‑gg                                                                                                                                                    >  1:1349178/1‑150 (MQ=255)
gtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTATAATCCGTTGCATGATGCGCTGGGCCAGGTGGAGGTTTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑gg                                                                                                                                                    <  1:890548/150‑1 (MQ=255)
    tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCT‑GGTTGc                                                                                                                                                <  2:53710/150‑1 (MQ=255)
      aTGACCCTGCGCAGCCCGCCCAGGTCGGTGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCCCGCTTCGTCGCGGCTGTTCAGGCTCGTTACGAAACAGAACTCGTCGCCACCACAGCGCGCAACTCCGTCGT‑GGTTGCgg                                                                                                                                              >  2:530613/1‑150 (MQ=255)
      aTGACCATGCGCCGCCCGCCCAGGTCGCTGGGGTCGTTCATCCGTTGCATGGTGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCACCCAGGTCGG‑GGTTGCgg                                                                                                                                              >  2:549523/1‑150 (MQ=255)
      aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCTTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCGCCAAAGCGCGCAACCAGGTCGT‑GGTTGCgg                                                                                                                                              <  2:1118048/150‑1 (MQ=255)
      aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGGCGT‑GGGTGCgg                                                                                                                                              >  2:1384569/1‑150 (MQ=255)
                   gcccgcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCT‑GGTTGCgc                                                                                                                                              >  2:638825/1‑136 (MQ=255)
                   gcccgcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCT‑GGTTGCgc                                                                                                                                              <  1:638825/137‑2 (MQ=255)
                    cccgcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGAt                                                                                                                                <  1:823939/150‑1 (MQ=255)
                    cccgcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTGCAGGCTGGTGACCACACAGAACTCGTCGCTACCAAGGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGAt                                                                                                                                <  1:549523/150‑1 (MQ=255)
                     ccgcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGAtg                                                                                                                               >  2:1480598/1‑150 (MQ=255)
                       gcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGAtgtg                                                                                                                             <  1:267410/150‑1 (MQ=255)
                       gcccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGAtgtg                                                                                                                             <  1:319258/150‑1 (MQ=255)
                         ccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGt                                                                                                                           >  1:22877/1‑150 (MQ=255)
                          cAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGtt                                                                                                                          <  1:482690/150‑1 (MQ=255)
                             gTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTCTTTGGTAGCGGCCTTGATGTGGTTGGc                                                                                                                       <  1:1485076/150‑1 (MQ=255)
                              tCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGAGGGTAGCGGCATTGATGTGGTTGGCg                                                                                                                      >  1:1898906/1‑150 (MQ=255)
                               cGATGGGGGCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGa                                                                                                                     <  1:565648/150‑1 (MQ=255)
                                gATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGAt                                                                                                                    <  2:128227/150‑1 (MQ=255)
                                      gTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                              <  1:570597/150‑1 (MQ=255)
                                      gTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                              >  1:1479721/1‑150 (MQ=255)
                                         tttCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGag                                                                                                           >  1:1358351/1‑150 (MQ=255)
                                         tttCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAG‑GCTGAAGAGCATCCTAGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGag                                                                                                           >  2:1483604/1‑150 (MQ=255)
                                           tCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagca                                                                                                         >  2:1465322/1‑150 (MQ=255)
                                             aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGAGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                       >  1:1760136/1‑150 (MQ=255)
                                              tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCt                                                                                                      >  1:250945/1‑150 (MQ=255)
                                              tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCt                                                                                                      >  1:1441345/1‑150 (MQ=255)
                                                cGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCg                                                                                                    >  2:153722/1‑150 (MQ=255)
                                                    gCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCAc                                                                                                <  1:1414187/150‑1 (MQ=255)
                                                     cATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCAcc                                                                                               <  1:1121273/150‑1 (MQ=255)
                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                              >  2:1037471/1‑150 (MQ=255)
                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGAGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                              <  2:1760136/150‑1 (MQ=255)
                                                             gcTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTc                                                                                                     <  1:1489535/136‑1 (MQ=255)
                                                             gcTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTc                                                                                                     >  2:1489535/1‑136 (MQ=255)
                                                             gcTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGAGGTTGGCGATCACCGTGAGCATCGCGTCACCGGcgtcgt                                                                                       >  1:746182/1‑150 (MQ=255)
                                                                     aGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAgg                                                                               <  1:1384569/150‑1 (MQ=255)
                                                                     aGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAgg                                                                               <  2:1373889/150‑1 (MQ=255)
                                                                                ttCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCgg                                                                                                                                              >  2:1324081/1‑76 (MQ=255)
                                                                                ttCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCgg                                                                                                                                              <  1:1324081/76‑1 (MQ=255)
                                                                                      gcgGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCtt                                                              >  1:1660115/1‑150 (MQ=255)
                                                                                        ggCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCa                                                                                                 <  2:1777841/113‑1 (MQ=255)
                                                                                        ggCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCa                                                                                                 >  1:1777841/1‑113 (MQ=255)
                                                                                        ggCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGa                                                            >  1:847868/1‑150 (MQ=255)
                                                                                          cTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaaa                                                          <  2:847868/150‑1 (MQ=255)
                                                                                            gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATg                                                        <  2:250945/150‑1 (MQ=255)
                                                                                             ttCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATgg                                                       <  1:153722/150‑1 (MQ=255)
                                                                                                  gCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGAt                                                  >  1:1855283/1‑150 (MQ=255)
                                                                                                   cTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATg                                                 >  1:1511967/1‑150 (MQ=255)
                                                                                                    tGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAGCCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGt                                                >  2:1600682/1‑150 (MQ=255)
                                                                                                     ggTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTc                                               <  1:164876/150‑1 (MQ=255)
                                                                                                         aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTg                                                             >  2:1846124/1‑132 (MQ=255)
                                                                                                         aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTg                                                             <  1:1846124/132‑1 (MQ=255)
                                                                                                         aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAgg                                           <  1:1252670/150‑1 (MQ=255)
                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                              >  1:730736/1‑81 (MQ=255)
                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                              <  2:730736/81‑1 (MQ=255)
                                                                                                            acacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                      <  1:1479891/136‑1 (MQ=255)
                                                                                                            acacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                      >  2:1479891/1‑136 (MQ=255)
                                                                                                             cacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACa                                       >  2:109972/1‑150 (MQ=255)
                                                                                                                aGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGa                                    <  2:1368811/150‑1 (MQ=255)
                                                                                                                aGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGa                                    <  2:822001/150‑1 (MQ=255)
                                                                                                                    cTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGc                                <  2:1479721/150‑1 (MQ=255)
                                                                                                                           ccaccaAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                        >  2:1111520/1‑135 (MQ=255)
                                                                                                                           ccaccaAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                        <  1:1111520/135‑1 (MQ=255)
                                                                                                                               caAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtttt                     <  1:1282975/150‑1 (MQ=255)
                                                                                                                                   gcgcgcAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCt                 <  2:22877/150‑1 (MQ=255)
                                                                                                                                   gcgcgcAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGCTGGCGATCACCTTGAGCAGATCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCt                 <  2:452799/150‑1 (MQ=255)
                                                                                                                                    cgcgcAACCAGGTCGT‑GGGTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTg                <  2:1443905/150‑1 (MQ=255)
                                                                                                                                      cgcAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGc              <  2:156549/150‑1 (MQ=255)
                                                                                                                                      cgcAACCACGACCT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGc              >  1:1158777/1‑150 (MQ=255)
                                                                                                                                             aGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGt       <  2:1430016/150‑1 (MQ=255)
                                                                                                                                              ggTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                              >  1:1830143/1‑46 (MQ=255)
                                                                                                                                              ggTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                              <  2:1830143/46‑1 (MQ=255)
                                                                                                                                               gTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTc     >  1:1652858/1‑150 (MQ=255)
                                                                                                                                                  gT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAc  <  2:1855283/150‑1 (MQ=255)
                                                                                                                                                  gT‑GGGTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGCACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAc  <  2:746182/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGT‑GGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAC  >  NC_002947/1310307‑1310602

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