Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,315 G→T 6.3% P169H (CCC→CAC)  PP_3505 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,3150GT 6.3% 251.5 / 5.2 80P169H (CCC→CAC) PP_3505hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base T (3/2);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

AACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGTGCCGCGCCT  >  NC_002947/3975168‑3975462
                                                                                                                                                   |                                                                                                                                                   
aaCTTCTCTGTTGCTGTTGTGGAAACTTTCCCAATTGCGACAAAGTGTCGCAACCGGCCGGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGgc                                                                                                                                                   >  2:1303111/1‑150 (MQ=255)
     ctctGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCCCGTGCGcctg                                                                                                                                              <  2:344363/150‑3 (MQ=255)
          ttgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCa                                                                                                                                         <  2:1136604/150‑1 (MQ=255)
                ttgtGGCACCTTTCCCAATTGCGAAAAAGTGTCGCCACCTGCCCGTTCAATCCATAACCCTGAATTGTTGCAAGTGCTGGTAATCAGACTGCAATTGTTCTGCCAGAACCTGCACGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCgc                                                                                                                                   >  2:1002999/1‑149 (MQ=255)
                 tgtgGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAc                                                                                                                                  <  1:196065/150‑1 (MQ=255)
                   tgGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAc                                                                                                                                  <  1:1277522/148‑1 (MQ=255)
                   tgGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAc                                                                                                                                  >  2:1277522/1‑148 (MQ=255)
                     gCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaac                                                                                                                              >  2:963617/1‑150 (MQ=255)
                     gCAACTTTCCCAATTGCGACAAAGTGTAGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaac                                                                                                                              >  1:332730/1‑150 (MQ=255)
                       aaCTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGAGGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTGCCATTTCGACcaaaaa                                                                                                                            >  1:139749/1‑147 (MQ=255)
                            tCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGGCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTgg                                                                                                                       <  2:218525/150‑1 (MQ=255)
                                aaTTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCa                                                                                                                   >  1:57962/1‑150 (MQ=255)
                                 aTTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCCACAACGTGgcgccc                                                                                                                  <  2:766281/150‑3 (MQ=255)
                                  ttGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACg                                                                                                                 <  2:245434/150‑1 (MQ=255)
                                   tGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACgg                                                                                                                <  2:851769/150‑1 (MQ=255)
                                      gACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCATGAATTGTTGCAAGGGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCAt                                                                                                             >  1:281859/1‑150 (MQ=255)
                                              gtCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCa                                                                                                     <  2:482807/150‑1 (MQ=255)
                                                  taaCCTGCCAGTCCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGt                                                                                                 <  1:1002999/149‑1 (MQ=255)
                                                  cAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCg                                                                                                           <  2:1157293/140‑1 (MQ=255)
                                                  cAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCg                                                                                                           >  1:1157293/1‑140 (MQ=255)
                                                   aaCCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaacaa                                                                                                                            >  2:46391/1‑122 (MQ=255)
                                                   aaCCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaacaa                                                                                                                            <  1:46391/122‑1 (MQ=255)
                                                     ccTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTAtt                                                                                              <  2:364387/150‑1 (MQ=255)
                                                      cTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGTCATCGGCTGCACGGTATtg                                                                                             <  1:963617/150‑1 (MQ=255)
                                                      cTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATtg                                                                                             >  1:746813/1‑150 (MQ=255)
                                                       tGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATtgc                                                                                            >  2:941493/1‑150 (MQ=255)
                                                        gCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATtgct                                                                                           >  2:1056755/1‑150 (MQ=255)
                                                             ttCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATtgctgtgct                                                                                      >  2:808242/1‑150 (MQ=255)
                                                                aaTCCATCCCCCTGTATTGTCGCAAGTGCTGGTAATCAGCCTGCAATCGTTCTGCCAGACGCTGCCCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACCACGTGGCGCACGGCATCGGCTGCACGCTATTGCTGTGCTGCa                                                                                   <  2:11283/150‑1 (MQ=255)
                                                                aaTCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCa                                                                                   >  1:1145456/1‑150 (MQ=255)
                                                                aaTCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCa                                                                                   >  1:177181/1‑150 (MQ=255)
                                                                  tCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaacaa                                                                                                                            >  2:1326044/1‑107 (MQ=255)
                                                                  tCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaacaa                                                                                                                            <  1:1326044/107‑1 (MQ=255)
                                                                   ccATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGc                                                                                >  2:490935/1‑150 (MQ=255)
                                                                    cATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGgcgc                                                                                                                                                 >  2:490038/1‑84 (MQ=255)
                                                                    cATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGgcgc                                                                                                                                                 <  1:490038/84‑1 (MQ=255)
                                                                    cATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCg                                                                               >  2:279455/1‑150 (MQ=255)
                                                                      tCACCCTGAATTGTTGCAAGGGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGAc                                                                             >  2:703441/1‑150 (MQ=255)
                                                                            tGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCgg                                                                       <  1:315504/150‑1 (MQ=255)
                                                                             gAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGt                                                                      <  1:1056755/150‑1 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCCCGTGCGcctgg                                                                                                                                             <  2:42807/73‑4 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCCCGTGCGcctgg                                                                                                                                             >  1:42807/1‑70 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTg                                                                                                                        <  2:753377/94‑1 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTg                                                                                                                        >  1:753377/1‑94 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCg                                                                           <  1:304052/139‑1 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCg                                                                           >  2:304052/1‑139 (MQ=255)
                                                                                   ttgCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAacac                                                                >  2:920173/1‑150 (MQ=255)
                                                                                             tGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGc                                                      >  2:1313247/1‑150 (MQ=255)
                                                                                               gTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCaa                                                    >  1:610660/1‑150 (MQ=255)
                                                                                                     aGCCTGCAAGTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGTCAGTTCCATATCGACCAACAACGTCGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCctc                                              <  2:338659/150‑1 (MQ=255)
                                                                                                            aaTTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGa                                       <  1:279455/150‑1 (MQ=255)
                                                                                                               tGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATgg                                    <  2:281859/150‑1 (MQ=255)
                                                                                                                gTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATggg                                   <  2:1270578/150‑1 (MQ=255)
                                                                                                                  tCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCATTTTCATATCGACCAACAACGTGGCGCCCGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGCAAGGGCAAGGCCTCTTGCGGATGGGCt                                 <  1:1046497/150‑1 (MQ=255)
                                                                                                                    tGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCg                               >  1:803295/1‑150 (MQ=255)
                                                                                                                     gCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATtgct                                                                                           <  2:784409/89‑1 (MQ=255)
                                                                                                                     gCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATtgct                                                                                           >  1:784409/1‑89 (MQ=255)
                                                                                                                      ccAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTTCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCt                             >  1:471552/1‑150 (MQ=255)
                                                                                                                      ccAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAg                                                                                  >  2:671700/1‑97 (MQ=255)
                                                                                                                      ccAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAg                                                                                  <  1:671700/97‑1 (MQ=255)
                                                                                                                      ccAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCt                             >  2:844779/1‑150 (MQ=255)
                                                                                                                       cAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTg                            <  1:437669/150‑1 (MQ=255)
                                                                                                                        aGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGt                                                                                                                         <  2:842656/56‑1 (MQ=255)
                                                                                                                        aGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGt                                                                                                                         >  1:842656/1‑56 (MQ=255)
                                                                                                                         gACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCaa                                                    <  2:10297/124‑1 (MQ=255)
                                                                                                                         gACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCaa                                                    >  1:10297/1‑124 (MQ=255)
                                                                                                                            gCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCcgcg                       >  2:1146677/1‑150 (MQ=255)
                                                                                                                             cTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCgcgc                      <  2:332730/150‑1 (MQ=255)
                                                                                                                              tGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATgg                                    >  1:334346/1‑135 (MQ=255)
                                                                                                                              tGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATgg                                    <  2:334346/135‑1 (MQ=255)
                                                                                                                                cgcgcgTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTg                   >  2:1418738/1‑150 (MQ=255)
                                                                                                                                    cgTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGc               <  1:920173/150‑1 (MQ=255)
                                                                                                                                      tCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGGTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCa             >  1:297509/1‑150 (MQ=255)
                                                                                                                                       aactggCGCCCGTGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCAGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAg            >  2:479729/5‑150 (MQ=255)
                                                                                                                                        ccAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGt           >  2:675779/1‑150 (MQ=255)
                                                                                                                                        ccAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGt           >  1:140347/1‑150 (MQ=255)
                                                                                                                                         ctggCGCCCGTGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGTg          >  1:974231/3‑150 (MQ=255)
                                                                                                                                              gcACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGTGCcgcg     >  1:1276797/1‑150 (MQ=255)
                                                                                                                                               cACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGTGCCgcgc    <  1:675779/150‑1 (MQ=255)
                                                                                                                                                 cGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGTGCCGCGCCt  >  2:695642/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATCGGCTGCACGGTATTGCTGTGCTGCAGGCGACCGTCGGTGAACACCAGGCAACGCAAGGCCTCTTGCGGATGGGCTCGCTGCCGCGCCTGTAGCCAGTGCCGCGCCT  >  NC_002947/3975168‑3975462

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: