Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,409,013 C→A 7.2% H322N (CAT→AAT)  ydaN → putative Zn(II) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,409,0130CA 7.2% 186.8 / 2.4 62H322N (CAT→AAT) ydaNputative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base A (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ATGGTCTTTTTACCCAGTACCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTT  >  NC_000913/1408872‑1409160
                                                                                                                                             |                                                                                                                                                   
aTGGTCTTTTTACCCAGTACCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCgta                                                                                                                                               <  4:77113/148‑1 (MQ=255)
aTGGTCTTTTTACCCAGTACCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCg                                                                                                                                                 >  7:256309/1‑146 (MQ=255)
           aCCCAGTACCTTTATGACCTGGTTATTTGGCGTCAACCTTTGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGGGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGtt                                                                                                                                     <  1:312437/147‑1 (MQ=255)
             ccAGTACCTTGCTGCCAGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTa                                                                                                                                  <  1:147485/148‑1 (MQ=255)
                 tACCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACa                                                                                                                               >  1:203130/1‑147 (MQ=255)
                  aCCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACaa                                                                                                                              >  1:266937/1‑147 (MQ=255)
                  aCCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACaa                                                                                                                              >  2:301591/1‑147 (MQ=255)
                      ttCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGc                                                                                                                          >  6:118208/1‑147 (MQ=255)
                                 ttATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACaaa                                                                                                                             >  7:198143/1‑133 (MQ=255)
                                 ttATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACaaa                                                                                                                             <  8:198143/133‑1 (MQ=255)
                                         cGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGta                                                                                                      <  2:203130/148‑1 (MQ=255)
                                            cAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACaaa                                                                                                                             <  8:232365/122‑1 (MQ=255)
                                            cAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACaaa                                                                                                                             >  7:232365/1‑122 (MQ=255)
                                            cAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGtata                                                                                                    >  2:197312/1‑147 (MQ=255)
                                            cAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTACGGCTGTCTGtata                                                                                                    >  6:190472/1‑147 (MQ=255)
                                              aCCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACaa                                                                                                 <  7:137330/148‑1 (MQ=255)
                                               ccTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTTGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGATGTCTTTATACaa                                                                                                 >  4:77871/1‑147 (MQ=255)
                                                cTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACaaaa                                                                                               <  1:301591/148‑1 (MQ=255)
                                                   ggtggGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACaaaaa                                                                                              <  8:256309/146‑1 (MQ=255)
                                                        gATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACaaaa                                                                                               >  2:260974/1‑140 (MQ=255)
                                                        gATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACaaaa                                                                                               <  1:260974/140‑1 (MQ=255)
                                                               ggcggcGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGttt                                                                                 >  1:115763/1‑147 (MQ=255)
                                                                      ggTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCAGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAg                                                                         >  8:235234/1‑148 (MQ=255)
                                                                        tGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGt                                                                        >  8:26779/1‑147 (MQ=255)
                                                                        tGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTc                                                                       <  2:115763/148‑1 (MQ=255)
                                                                        tGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTc                                                                       <  1:197312/148‑1 (MQ=255)
                                                                         ggCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGt                                                                                   >  7:53354/1‑135 (MQ=255)
                                                                         ggCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGt                                                                                   <  8:53354/135‑1 (MQ=255)
                                                                         ggCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTc                                                                       >  6:167924/1‑147 (MQ=255)
                                                                           cAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGa                                                                           >  7:272199/1‑141 (MQ=255)
                                                                           cAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGa                                                                           <  8:272199/141‑1 (MQ=255)
                                                                            aaTTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCatt                                                                    >  3:6427/1‑145 (MQ=255)
                                                                             aTTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATaa                                                                  <  2:266937/148‑1 (MQ=255)
                                                                              ttCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATaaa                                                                 >  8:169079/1‑148 (MQ=255)
                                                                              tgcgttttttCAATTTTTTTTATTCTTTGAGTTTTTCTTATTGTTTGTGTTGCTTTATGGTTTAATCGTAGTAAATGGTTGTAAAAAAATATATTTTTCCGGCTGTCTGTATTCAAAATCGGCTCAAAGTTTGAGCGAAGTCAATaa                                                                  <  3:79540/141‑1 (MQ=255)
                                                                               tCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTAAGCGAAGTCAATAAAc                                                                >  5:56388/1‑148 (MQ=255)
                                                                                 ggATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAActc                                                              <  5:167924/148‑1 (MQ=255)
                                                                                 ggATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTATCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAAct                                                               <  5:118208/147‑1 (MQ=255)
                                                                                  gtttttCAATTTTTTGTATTCTGTTAGTTGTCCTTATTTGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTTTAACAAAAGCAATTTTTCCGGCTGTATGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAActct                                                             <  5:140377/146‑1 (MQ=255)
                                                                                       tCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGATGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAATCTCTCTAcc                                                        >  4:193712/1‑148 (MQ=255)
                                                                                            tttttGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTg                                                                                <  8:158260/119‑1 (MQ=255)
                                                                                            tttttGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTg                                                                                >  7:158260/1‑119 (MQ=255)
                                                                                                     tCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCa                                             >  8:296957/1‑145 (MQ=255)
                                                                                                     tCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCa                                             <  7:296957/145‑1 (MQ=255)
                                                                                                     tCTGTTAGTGGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGGTGTAACAAAAGCAATTTTTCCGTCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAat                                           <  5:190472/147‑1 (MQ=255)
                                                                                                         ttAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATAtctc                                      >  8:19618/1‑148 (MQ=255)
                                                                                                         ttAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATAtc                                        <  7:169079/146‑1 (MQ=255)
                                                                                                                ttCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAATCTCTCTACCCATTCAGGGCAATATCTCTCTTGCa                               <  3:193712/148‑1 (MQ=255)
                                                                                                                    tATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGt                            >  2:169142/1‑147 (MQ=255)
                                                                                                                    tATTGGTGGTGTTGCTTTATGGTTGAATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCtt                                  >  1:65280/1‑141 (MQ=255)
                                                                                                                    tATTGGTGGTGTTGCTTTATGGTTGAATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCtt                                  <  2:65280/141‑1 (MQ=255)
                                                                                                                          tggtgTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACTCCGCAAAGTTTTAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATAACTCTCTTCCAGGTGAATGc                      >  8:310510/1‑147 (MQ=255)
                                                                                                                          tggtgTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCACTAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCa                     <  4:6427/148‑1 (MQ=255)
                                                                                                                          tggtgTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCa                     <  8:184532/148‑1 (MQ=255)
                                                                                                                             tgtTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCa                                                                      <  6:44646/96‑1 (MQ=255)
                                                                                                                             tgtTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCa                                                                      >  5:44646/1‑96 (MQ=255)
                                                                                                                              gtTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGt                 >  2:262435/1‑148 (MQ=255)
                                                                                                                                    ttATGGTTGAATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATAtc                                        >  2:298042/1‑119 (MQ=255)
                                                                                                                                    ttATGGTTGAATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATAtc                                        <  1:298042/119‑1 (MQ=255)
                                                                                                                                     tATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCg           <  7:26779/147‑1 (MQ=255)
                                                                                                                                           tGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCg    >  2:175983/1‑148 (MQ=255)
                                                                                                                                            gCATCGTAGTAAATGGTTTTAACAAAAGCAATTTTTCCGGATGTCTGTATACAAAAACGCCGCAAAGTTTGAGAGAATTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCg    <  3:77871/147‑1 (MQ=255)
                                                                                                                                             cATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGtt  <  7:235234/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ATGGTCTTTTTACCCAGTACCTTTCTGACCGGGTTATTTGGCGTCAACCTTGGTGGGATCCCTGGCGGCGGGTGGCAATTCGGATTTTCAATTTTTTGTATTCTGTTAGTTGTTCTTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTT  >  NC_000913/1408872‑1409160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: