Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,980,772 T→A 5.9% intergenic (+76/‑115) wecG → / → yifK UDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,7720TA 5.9% 430.8 / 11.2 120intergenic (+76/‑115)wecG/yifKUDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (59/53);  new base A (3/4);  total (63/57)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 129 total aligned reads displayed.

CAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTG  >  NC_000913/3980625‑3980917
                                                                                                                                                   |                                                                                                                                                 
cAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAtt                                                                                                                                                  >  3:1096935/1‑149 (MQ=255)
cAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAt                                                                                                                                                   <  4:522263/148‑1 (MQ=255)
cAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAt                                                                                                                                                   >  3:522263/1‑148 (MQ=255)
                 gCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGc                                                                                                                                 <  2:462428/149‑1 (MQ=255)
                   gTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCa                                                                                                                               >  1:1502892/1‑149 (MQ=255)
                                  tgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGc                                                                                                                                          >  2:1481868/1‑123 (MQ=255)
                                  tgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGc                                                                                                                                          <  1:1481868/123‑1 (MQ=255)
                                         tttACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTa                                                                                                                 >  1:315302/1‑141 (MQ=255)
                                         tttACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTa                                                                                                                 <  2:315302/141‑1 (MQ=255)
                                              gCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCtt                                                                                                                                        >  1:813364/1‑113 (MQ=255)
                                              gCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCtt                                                                                                                                        <  2:813364/113‑1 (MQ=255)
                                               cTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGaaaa                                                                                                           >  4:1588165/1‑141 (MQ=255)
                                                  gCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCaa                                                                                                  <  3:1071678/147‑1 (MQ=255)
                                                      tACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                                                                    >  1:1585626/1‑125 (MQ=255)
                                                      tACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                                                                    <  2:1585626/125‑1 (MQ=255)
                                                           cGGCGACCTATGATTTTCCTCTCTATGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGaaa                                                                                                            <  1:1058900/128‑1 (MQ=255)
                                                                        ttttCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACa                                                                                                         >  3:1789929/1‑118 (MQ=255)
                                                                            ccTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAAt                                                                                                                                            >  1:861000/1‑79 (MQ=255)
                                                                            ccTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAAt                                                                                                                                            <  2:861000/79‑1 (MQ=255)
                                                                                       aaaGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAaca                                                                                                >  4:1169352/1‑112 (MQ=255)
                                                                                       aaaGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAaca                                                                                                <  3:1169352/112‑1 (MQ=255)
                                                                                        aaGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATCTATCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                                                          <  4:1096935/149‑1 (MQ=255)
                                                                                          gCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGc                                                        >  3:510856/1‑149 (MQ=255)
                                                                                            gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                                                          >  3:225027/1‑145 (MQ=255)
                                                                                            gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                                                          <  4:225027/145‑1 (MQ=255)
                                                                                               ggAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                                                          <  3:808569/142‑1 (MQ=255)
                                                                                               ggAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                                                          >  4:808569/1‑142 (MQ=255)
                                                                                                  gTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATaa                                                                                <  1:945307/117‑1 (MQ=255)
                                                                                                  gTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATaa                                                                                >  2:945307/1‑117 (MQ=255)
                                                                                                  gTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaac                                                >  2:1512745/1‑149 (MQ=255)
                                                                                                     gTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACAttt                                                                                                      >  3:931330/1‑92 (MQ=255)
                                                                                                      tCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTAc                                                                                                                <  2:204907/81‑1 (MQ=255)
                                                                                                      tCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTAc                                                                                                                >  1:204907/1‑81 (MQ=255)
                                                                                                      tCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt                                                      >  1:370226/1‑139 (MQ=255)
                                                                                                      tCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt                                                      <  2:370226/139‑1 (MQ=255)
                                                                                                        aCTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACAttt                                                                                                      >  3:189925/1‑89 (MQ=255)
                                                                                                        aCTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACAttt                                                                                                      <  4:189925/89‑1 (MQ=255)
                                                                                                        aCTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt                                                      <  2:997156/137‑1 (MQ=255)
                                                                                                        aCTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt                                                      >  1:997156/1‑137 (MQ=255)
                                                                                                        aCTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa                                                 <  3:413775/142‑1 (MQ=255)
                                                                                                          tCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATaa                                                                                <  2:456235/109‑1 (MQ=255)
                                                                                                          tCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATaa                                                                                >  1:456235/1‑109 (MQ=255)
                                                                                                           ccGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGa                                       <  2:313261/149‑1 (MQ=255)
                                                                                                                gCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGt                                                                     <  4:1289958/114‑1 (MQ=255)
                                                                                                                gCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGa                                       >  3:781945/1‑144 (MQ=255)
                                                                                                                gCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTGATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGt                                                                     >  3:1289958/1‑114 (MQ=255)
                                                                                                                  cTGTCTGCTTTTTTATTACACATAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGaa                                                                                                             <  1:699948/72‑1 (MQ=255)
                                                                                                                  cTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGaa                                                                                                             >  2:699948/1‑72 (MQ=255)
                                                                                                                      cTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCa                                                       <  3:1073163/122‑1 (MQ=255)
                                                                                                                      cTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCa                                                       >  4:1073163/1‑122 (MQ=255)
                                                                                                                       tGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa                                            >  3:865571/1‑132 (MQ=255)
                                                                                                                       tGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa                                            <  4:865571/132‑1 (MQ=255)
                                                                                                                       tGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa                                            <  2:1833165/132‑1 (MQ=255)
                                                                                                                        gCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa                                            >  2:72873/1‑131 (MQ=255)
                                                                                                                        gCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa                                            <  1:72873/131‑1 (MQ=255)
                                                                                                                          ttttttATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt                                                      <  3:1527641/119‑1 (MQ=255)
                                                                                                                          ttttttATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt                                                      >  4:1527641/1‑119 (MQ=255)
                                                                                                                          ttttttATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa                                                 >  2:1559154/1‑124 (MQ=255)
                                                                                                                          ttttttATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGAt                        >  4:316029/1‑149 (MQ=255)
                                                                                                                                attaCACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACaaa                  <  1:1512745/149‑1 (MQ=255)
                                                                                                                                attaCACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACaa                   >  2:226008/1‑148 (MQ=255)
                                                                                                                                attaCACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACaa                   <  1:226008/148‑1 (MQ=255)
                                                                                                                                 ttaCACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa                                                 >  1:394927/1‑117 (MQ=255)
                                                                                                                                 ttaCACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAAc                 >  1:1810076/1‑149 (MQ=255)
                                                                                                                                   acacAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAg                                      >  4:51354/1‑126 (MQ=255)
                                                                                                                                   acacAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAg                                      <  3:51354/126‑1 (MQ=255)
                                                                                                                                   aaataaaGCATTAAAAATTTGAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATa                       >  4:824688/5‑141 (MQ=255)
                                                                                                                                    cacaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGa                                       >  1:1656898/1‑124 (MQ=255)
                                                                                                                                    cacaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGa                                       <  2:1467220/124‑1 (MQ=255)
                                                                                                                                    cacaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCag              >  2:95125/1‑149 (MQ=255)
                                                                                                                                    aataaaGCATTAAAAATTTGAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACAcc                                                               >  4:1269123/4‑100 (MQ=38)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAg                                                         <  3:100894/105‑1 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAg                                                         >  4:100894/1‑105 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa                                                 >  2:1287069/1‑113 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa                                                 <  1:1287069/113‑1 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaaca                                          >  4:1045875/1‑120 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCaga             >  2:475691/1‑149 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCaga             >  2:1719545/1‑149 (MQ=255)
                                                                                                                                     acaAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCaga             <  1:1719545/149‑1 (MQ=255)
                                                                                                                                      caAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTAt                               <  2:1286105/130‑1 (MQ=255)
                                                                                                                                      caAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTAt                               >  1:1286105/1‑130 (MQ=255)
                                                                                                                                       aaaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaca                                                 <  3:1207487/111‑3 (MQ=255)
                                                                                                                                       aaaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa                                                 >  4:1207487/1‑111 (MQ=255)
                                                                                                                                       aaaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGAttt                                 <  1:756347/127‑1 (MQ=255)
                                                                                                                                       aaaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGAttt                                 >  2:756347/1‑127 (MQ=255)
                                                                                                                                       aaaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTAt                               >  3:758701/1‑129 (MQ=255)
                                                                                                                                       aaaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTAt                               <  4:758701/129‑1 (MQ=255)
                                                                                                                                        aaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATaa                      <  1:104076/137‑1 (MQ=255)
                                                                                                                                        aaGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATaa                      >  2:104076/1‑137 (MQ=255)
                                                                                                                                         aGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTa         >  3:1119708/1‑149 (MQ=255)
                                                                                                                                          gCATTAAAAATTTGAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGa                                       <  4:1109888/118‑1 (MQ=255)
                                                                                                                                          gCATTAAAAATTTGAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACaaa                  >  2:903769/1‑139 (MQ=255)
                                                                                                                                          gCATTAAAAATTTGAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACaaa                  <  1:903769/139‑1 (MQ=255)
                                                                                                                                           cATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATg                              <  1:779788/126‑1 (MQ=255)
                                                                                                                                              tCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCagag            >  2:761460/1‑141 (MQ=255)
                                                                                                                                                aaaTTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTg  >  4:284055/1‑149 (MQ=255)
                                                                                                                                                   tttttAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATa                                                                                 <  3:1150224/67‑1 (MQ=255)
                                                                                                                                                   tttttAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCgg                                                             <  1:210441/87‑1 (MQ=255)
                                                                                                                                                   tttttAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCaga             >  3:1488905/1‑135 (MQ=255)
                                                                                                                                                   tttttAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCaga             <  4:1488905/135‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTG  >  NC_000913/3980625‑3980917

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