Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,756,563 A→G 5.8% intergenic (‑52/+113) aldB ← / ← yiaY aldehyde dehydrogenase B/L‑threonine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,756,5630AG 5.8% 660.9 / 20.3 191intergenic (‑52/+113)aldB/yiaYaldehyde dehydrogenase B/L‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (89/91);  new base G (6/5);  total (95/96)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

Only 100 of 221 total aligned reads displayed.

TCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTT  >  NC_000913/3756417‑3756710
                                                                                                                                                  |                                                                                                                                                   
   ccgccAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 <  2:257670/148‑1 (MQ=255)
     gccAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCg                                                                                                                                               <  3:452744/148‑1 (MQ=255)
      ccAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGc                                                                                                                                              >  4:717481/1‑148 (MQ=255)
      ccAATAAAGTTATCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCg                                                                                                                                               >  1:861874/1‑147 (MQ=255)
               ttGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 >  2:955652/1‑136 (MQ=255)
                tGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaaa                                                                                                                                                     <  3:686966/131‑1 (MQ=255)
                tGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaaa                                                                                                                                                     >  4:686966/1‑131 (MQ=255)
                     tAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 <  2:290516/130‑1 (MQ=255)
                      aGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAAc                                                                                                                                   >  1:942491/1‑143 (MQ=255)
                              tttAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCttt                                                                                                                             >  2:511870/1‑141 (MQ=255)
                                     tGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCt                                                                                                                                                  <  3:748006/113‑1 (MQ=255)
                                         ggggAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCgg                                                                                                                                       >  1:1079785/1‑120 (MQ=255)
                                         ggggAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCgg                                                                                                                                       <  2:1079785/120‑1 (MQ=255)
                                          gggAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 >  4:513937/1‑109 (MQ=255)
                                           ggAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAAttt                                                                                                         >  1:356762/1‑148 (MQ=255)
                                             aaaCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTGACGACAGCTATCACGAATTTAc                                                                                                       >  1:813984/1‑148 (MQ=255)
                                             aaaCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCAGGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTAGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTAc                                                                                                       <  2:813984/148‑1 (MQ=255)
                                               aCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 >  1:617474/1‑104 (MQ=255)
                                                   tACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                             >  3:530762/1‑136 (MQ=255)
                                                   tACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTa                                                                                                        <  1:132407/141‑1 (MQ=255)
                                                       cGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCg                                                                                                                                               <  3:428615/98‑1 (MQ=255)
                                                       cGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCg                                                                                                                                               >  4:428615/1‑98 (MQ=255)
                                                       cGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGATCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGt                                                                                             >  4:1080916/1‑148 (MQ=255)
                                                              cTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACg                                                                                                      >  1:304345/1‑132 (MQ=255)
                                                               ttAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCa                                                                                                  <  1:587195/135‑1 (MQ=255)
                                                                   tCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGt                                                                                 <  2:356762/148‑1 (MQ=255)
                                                                    cTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAAc                                                                                                                                   <  1:438581/97‑1 (MQ=255)
                                                                    cTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAAc                                                                                                                                   >  2:438581/1‑97 (MQ=255)
                                                                     tgtgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCg                                                                                                                                               >  3:329482/1‑84 (MQ=255)
                                                                         cTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCa                                                                                                                <  3:529859/111‑1 (MQ=255)
                                                                           gAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTa                                                                                                                                     <  4:472180/88‑1 (MQ=255)
                                                                                 ggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                             >  3:190312/1‑106 (MQ=255)
                                                                                 ggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                             <  4:190312/106‑1 (MQ=255)
                                                                                   attattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCACTCGGTAAc                                                                                                                                   <  3:1030174/82‑1 (MQ=255)
                                                                                   attattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCACTCGGTAAc                                                                                                                                   >  4:1030174/1‑82 (MQ=255)
                                                                                    ttattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 >  3:1141467/1‑67 (MQ=255)
                                                                                    ttattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 <  4:1141467/67‑1 (MQ=255)
                                                                                     tattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACag                                                                                        >  2:515345/1‑123 (MQ=255)
                                                                                        tGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTg                                                                               >  4:810193/1‑129 (MQ=255)
                                                                                           tCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGATTTCGCCAg                                                                                                                                           <  3:815142/66‑1 (MQ=38)
                                                                                           tCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                             >  4:200745/1‑96 (MQ=255)
                                                                                           tCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                             <  3:200745/96‑1 (MQ=255)
                                                                                           tCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACg                                                                                                      <  4:383067/103‑1 (MQ=255)
                                                                                           tCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACg                                                                                                      >  3:383067/1‑103 (MQ=255)
                                                                                           tCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGcc                                                         <  3:1080916/148‑1 (MQ=255)
                                                                                           gcATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAAc                                                                                                                                   <  3:104305/73‑1 (MQ=255)
                                                                                           acATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAAc                                                                                                                                   >  4:104305/2‑74 (MQ=255)
                                                                                            cATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGt                                                                                                                                          >  2:747245/1‑66 (MQ=255)
                                                                                            cATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCa                                                        >  1:1119757/1‑148 (MQ=255)
                                                                                                aTCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGaa                                                                                                            >  4:529920/1‑92 (MQ=255)
                                                                                                aTCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaa                                                                             >  4:531213/1‑123 (MQ=255)
                                                                                                aTCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTg                                                           <  3:883275/141‑1 (MQ=255)
                                                                                                aTCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCa                                                        >  2:481422/1‑144 (MQ=255)
                                                                                                  cGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                                 >  2:305727/1‑53 (MQ=255)
                                                                                                  cGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCa                                                                                                                <  3:918748/86‑1 (MQ=255)
                                                                                                    tctcCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCgg                                                                                                                                       >  1:512191/1‑61 (MQ=255)
                                                                                                    tctcCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAg                                                                                                >  4:342821/1‑100 (MQ=255)
                                                                                                     ctcCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTAc                                                                                                       >  1:479494/1‑92 (MQ=255)
                                                                                                      tcCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                         <  4:908050/121‑1 (MQ=255)
                                                                                                       ccTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTa                                                                                                                            >  4:86800/1‑69 (MQ=255)
                                                                                                        cTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTtg                                                 <  3:51338/143‑1 (MQ=255)
                                                                                                         ttGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGaaa                                                                                                                                                     >  2:615003/1‑42 (MQ=255)
                                                                                                         ttGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTg                                                           <  2:870070/132‑1 (MQ=255)
                                                                                                                gtgAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                         <  3:865312/111‑1 (MQ=255)
                                                                                                                gtgAGGTATAACAAGGGTAGACGCGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAggg                                                                                              >  1:898154/1‑90 (MQ=255)
                                                                                                                  gAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTa                                                                                           >  1:508978/1‑91 (MQ=255)
                                                                                                                  gAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGt                                                                                 <  4:926538/101‑1 (MQ=255)
                                                                                                                       ataACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGc                                                                                                                     >  2:899088/1‑60 (MQ=255)
                                                                                                                       ataACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGt                                                                                 <  1:733240/96‑1 (MQ=255)
                                                                                                                       ataACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGt                                                                                 >  2:733240/1‑96 (MQ=255)
                                                                                                                         aaCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                         <  3:326059/102‑1 (MQ=255)
                                                                                                                         aaCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                         >  4:326059/1‑102 (MQ=255)
                                                                                                                          aCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGc                                                      >  1:659118/1‑120 (MQ=255)
                                                                                                                          aCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGc                                                      <  2:659118/120‑1 (MQ=255)
                                                                                                                           cAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGcgca                         >  3:476027/1‑145 (MQ=255)
                                                                                                                             aGGGTAGACGTGACTAGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCg                                                                       >  4:851568/1‑100 (MQ=255)
                                                                                                                               ggTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGt                                                                                 <  4:1001064/88‑1 (MQ=255)
                                                                                                                               ggTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGt                                                                                 >  3:1001064/1‑88 (MQ=255)
                                                                                                                                 tAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGAt                                                                                    >  3:275271/1‑83 (MQ=255)
                                                                                                                                 tAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGAt                                                                                    <  4:275271/83‑1 (MQ=255)
                                                                                                                                         aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCg                                                                       <  2:769284/88‑1 (MQ=255)
                                                                                                                                         aCTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaa                                                                             >  3:585828/1‑82 (MQ=255)
                                                                                                                                          cTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGaa                                            <  3:520553/114‑1 (MQ=255)
                                                                                                                                           tGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTa                                                                   >  4:234305/1‑90 (MQ=255)
                                                                                                                                           tGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTa                                                                   <  3:234305/90‑1 (MQ=255)
                                                                                                                                           tGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCg                                                                       >  1:215588/1‑86 (MQ=255)
                                                                                                                                           tGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGcc                                                     >  1:759331/1‑104 (MQ=255)
                                                                                                                                            ggCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAg                                                                                                <  2:177787/60‑1 (MQ=255)
                                                                                                                                            ggCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAatgat                                       >  4:511665/1‑117 (MQ=255)
                                                                                                                                            ggCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACaa        >  4:603856/1‑148 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTa                                                                                                        >  4:656674/1‑48 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTa                                                                                                        <  3:656674/48‑1 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATa                                                                                   <  2:580989/69‑1 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGtt                                                  >  2:1128244/1‑102 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGt                                                   >  4:525665/1‑101 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcaca         <  3:365155/143‑1 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcaca         >  4:365155/1‑143 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATttc    >  3:951124/1‑148 (MQ=255)
                                                                                                                                                aaaTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATttc    >  1:557115/1‑148 (MQ=255)
                                                                                                                                                  aTCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATttctt  >  2:1114640/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTT  >  NC_000913/3756417‑3756710

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