Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,505,303 Δ1 bp 5.7% coding (382/1074 nt) PP_3103 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,505,3030C. 5.7% 225.3 / 7.8 70coding (382/1074 nt)PP_3103hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base . (2/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TGGTCGATCGGGTAGTGGTCGCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCG  >  NC_002947/3505154‑3505450
                                                                                                                                                     |                                                                                                                                                    
tggtCGATCGGGTAGTGGTCGCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGc                                                                                                                                                      <  1:1524894/150‑1 (MQ=255)
          ggTAGTGGTCGCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGa                                                                                                                                           <  2:36366/150‑1 (MQ=255)
           gTAGTGGTCGCAATTGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAg                                                                                                                                          >  1:857900/1‑150 (MQ=255)
                 gTCGGAAGGGCTCCGGTTGCGCGTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGGGATCTCGGCCTGCAGCCGAGCGCTcc                                                                                                                                    <  2:96526/150‑1 (MQ=255)
                   cGCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCgg                                                                                                                                  <  2:615065/150‑1 (MQ=255)
                    gCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGc                                                                                                                                                      >  1:561175/1‑130 (MQ=255)
                    gCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGc                                                                                                                                                      <  2:561175/130‑1 (MQ=255)
                      aaTGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCGGCTTGTTGTAGATCTCGGCCTGCAGCCGAGCGCTCCGggcg                                                                                                                               <  1:1506126/150‑1 (MQ=255)
                        tGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGggcggc                                                                                                                             >  2:406131/1‑150 (MQ=255)
                           cTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAAt                                                                                                                          <  2:1625063/150‑1 (MQ=255)
                            tCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCTGCAGGCCGAgatc                                                                                                                                       >  2:892614/1‑133 (MQ=255)
                            tCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCTGCAGGCCGAgatc                                                                                                                                       <  1:892614/136‑4 (MQ=255)
                               ggTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCtt                                                                                                                      <  2:135286/150‑1 (MQ=255)
                                 ttGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCCACAGTTTCTGCTCATTCTTCGCCAACGAGGCCGCCAGCTGCCTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTc                                                                                                                    <  2:1628897/150‑1 (MQ=255)
                                 ttGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTc                                                                                                                    <  1:1446650/150‑1 (MQ=255)
                                            cACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCa                                                                                                         >  1:1240805/1‑150 (MQ=255)
                                                   gaTCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGcc                                                                                                  >  1:576104/1‑150 (MQ=255)
                                                     tCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCgg                                                                                                <  2:857900/150‑1 (MQ=255)
                                                      cGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGc                                                                                               <  2:493797/150‑1 (MQ=255)
                                                        gtgtGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCgg                                                                                                                                  >  1:1081723/1‑113 (MQ=255)
                                                        gtgtGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCgg                                                                                                                                  <  2:1081723/113‑1 (MQ=255)
                                                             gCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGa                                                                                                                                           >  1:1413319/1‑99 (MQ=255)
                                                             gCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGa                                                                                                                                           <  2:1413319/99‑1 (MQ=255)
                                                             gCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATc                                                                                        >  1:632635/1‑150 (MQ=255)
                                                             gCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATc                                                                                        >  2:710361/1‑150 (MQ=255)
                                                             gCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATc                                                                                        >  2:201315/1‑150 (MQ=255)
                                                              ccTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCt                                                                                       <  1:1626566/150‑1 (MQ=255)
                                                                  tttGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGccc                                                                                   <  1:710361/150‑1 (MQ=255)
                                                                      ggCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGa                                                                               <  2:632635/150‑1 (MQ=255)
                                                                         tGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGc                                                                                                                                                      <  1:320676/77‑1 (MQ=255)
                                                                         tGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGc                                                                                                                                                      >  2:320676/1‑77 (MQ=255)
                                                                          gCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTaa                                                                           <  1:201315/150‑1 (MQ=255)
                                                                           cAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAag                                                                          >  1:118512/1‑150 (MQ=255)
                                                                             aCTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAaggt                                                                        <  1:406131/150‑1 (MQ=255)
                                                                              cTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAaggta                                                                       <  2:576104/150‑1 (MQ=255)
                                                                              cTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAaggta                                                                       >  1:1238828/1‑150 (MQ=255)
                                                                                         ccAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATgg                                                            <  1:1615564/150‑1 (MQ=255)
                                                                                          cAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGc                                                           >  1:1020897/1‑150 (MQ=255)
                                                                                            aCAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCGCCGAGTAAGGTAGGTGACGATGGCTg                                                         >  1:5355/1‑150 (MQ=255)
                                                                                                tttCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCaa                                                                                                                           <  1:1204088/80‑1 (MQ=255)
                                                                                                tttCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCaa                                                                                                                           >  2:1204088/1‑80 (MQ=255)
                                                                                                ttgcTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAg                                                     >  1:1641182/4‑150 (MQ=255)
                                                                                                   cTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTaa                                                  <  2:1240805/150‑1 (MQ=255)
                                                                                                    tGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAAc                                                 <  2:1238828/150‑1 (MQ=255)
                                                                                                    tGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAAc                                                 <  2:1298688/150‑1 (MQ=255)
                                                                                                     gCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCATTAAca                                                >  1:504485/1‑150 (MQ=255)
                                                                                                      cTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCt                                                                                       >  2:334227/1‑110 (MQ=255)
                                                                                                      cTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCt                                                                                       <  1:334227/110‑1 (MQ=255)
                                                                                                         aTTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAAcagcag                                            >  1:1360281/1‑150 (MQ=255)
                                                                                                           tcttcGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACg                                                                <  2:866782/128‑1 (MQ=255)
                                                                                                           tcttcGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACg                                                                >  1:866782/1‑128 (MQ=255)
                                                                                                            cttcGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGGGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCAGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCAc                                         <  1:1570496/150‑1 (MQ=255)
                                                                                                             ttcGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTc                                                                                                          <  2:672059/84‑1 (MQ=255)
                                                                                                             ttcGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTc                                                                                                          >  1:672059/1‑84 (MQ=255)
                                                                                                                gCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAAc                                     >  1:1403909/1‑150 (MQ=255)
                                                                                                                 ccAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAACGTAGGTGACGATGGCTGTCAGTAACAgcagcagatcgg                                    >  1:14845/1‑144 (MQ=255)
                                                                                                                   aaCCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGc                                  >  2:1256715/1‑150 (MQ=255)
                                                                                                                         gccgccAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGggc                                                                                                                                <  1:1295194/50‑1 (MQ=255)
                                                                                                                         gccgccAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGggc                                                                                                                                >  2:1295194/1‑50 (MQ=255)
                                                                                                                            gccAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAag                                                                          >  1:508905/1‑101 (MQ=255)
                                                                                                                            gccAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAag                                                                          <  2:508905/101‑1 (MQ=255)
                                                                                                                             ccAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCcgc                        >  1:1046609/1‑150 (MQ=255)
                                                                                                                             ccAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCcgc                        <  2:1360281/150‑1 (MQ=255)
                                                                                                                               aGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGa                      <  1:1256715/150‑1 (MQ=255)
                                                                                                                                       gttgGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCATTAACAGCAGCACCATCTGCTCCTGCACGCCAGcttgctc              <  2:504485/150‑7 (MQ=255)
                                                                                                                                        ttgGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGa                                                                               <  1:487761/84‑1 (MQ=255)
                                                                                                                                        ttgGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGa                                                                               >  2:487761/1‑84 (MQ=255)
                                                                                                                                         tgGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGc            >  2:925212/1‑150 (MQ=255)
                                                                                                                                          tctgaTCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTGGCCATGCTGCTCATGACGCCGGCCTTCATCCGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTgcac           <  1:836077/147‑4 (MQ=255)
                                                                                                                                          ggAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAag                                                                          >  2:413633/1‑87 (MQ=255)
                                                                                                                                          ggAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAag                                                                          <  1:413633/87‑1 (MQ=255)
                                                                                                                                              aTCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCCGCCCCCGAGTAAGGTAGGTGACGATGGATGTAAGTAAGCGCATCACGAACGACGCCGGCCCGCGAGCAAGCTgcaagaga       >  2:836077/1‑143 (MQ=255)
                                                                                                                                               gctcGGCTGCAGGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGt          <  2:620597/145‑1 (MQ=255)
                                                                                                                                               gctcGGCTGCAGGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGt          >  1:620597/2‑146 (MQ=255)
                                                                                                                                                   ggCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCTTCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAgcgcg  >  1:288376/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGTCGATCGGGTAGTGGTCGCAATGGCTCCGGTTGCGCTTGTGCACGAGAGATCGGTGTGGCCTCTTTGGGCTGCAACTCCGGCTGCGCCAACAGTTTCTGCTCATTCTTCGCCAACCAGGCCGCCAGCTGCTTGTTGGAGATCTCGGCCTGCAGCCGAGCGCTCCGGGCGGCAATGCTTTCCATGCTGCTCATGACGCCGGCCTTCATCTGCCCCCGAGTAAGGTAGGTGACGATGGCTGTCAGTAACAGCAGCACCAACTGCTCCTGCCCGCGAGCAAGCTGGCGTGCAGCGCG  >  NC_002947/3505154‑3505450

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: