Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,019,431 T→A 5.1% W91R (TGG→AGG)  ssnA → putative chlorohydrolase/aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,019,4310TA 5.1% 254.4 / 3.5 78W91R (TGG→AGG) ssnAputative chlorohydrolase/aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTG  >  NC_000913/3019284‑3019577
                                                                                                                                                   |                                                                                                                                                  
cccGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCtg                                                                                                                                                   <  4:7712/149‑1 (MQ=255)
       ccAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTggcgg                                                                                                                                            >  6:355073/1‑149 (MQ=255)
        cAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTggcg                                                                                                                                             >  8:356798/1‑147 (MQ=255)
        cAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTggcg                                                                                                                                             <  7:356798/147‑1 (MQ=255)
          gCTTCAAAGTGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCAGTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATATCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTc                                                                                                                                         <  1:58272/149‑1 (MQ=255)
          gCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTc                                                                                                                                         <  8:352145/149‑1 (MQ=255)
              cAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGAt                                                                                                                                      <  4:311681/148‑1 (MQ=255)
              cAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATc                                                                                                                                     >  6:221213/1‑149 (MQ=255)
               aaaGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATcg                                                                                                                                    >  3:242379/1‑149 (MQ=255)
                aaGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATcgc                                                                                                                                   >  5:255800/1‑149 (MQ=255)
                     aTGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCAGGTGGCGGCTCGATCGCGCCCt                                                                                                                              >  8:20913/1‑149 (MQ=255)
                      tGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCtt                                                                                                                             >  4:199536/1‑149 (MQ=255)
                         aTGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTtgat                                                                                                                          <  6:150187/149‑1 (MQ=255)
                            gCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGaa                                                                                                                       <  6:255800/149‑1 (MQ=255)
                              cGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAaga                                                                                                                     <  5:355073/149‑1 (MQ=255)
                                gATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGt                                                                                                                   <  6:174423/149‑1 (MQ=255)
                                  ttGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCg                                                                                                                 >  3:295058/1‑149 (MQ=255)
                                           cAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTa                                                                                                         >  3:185963/1‑148 (MQ=255)
                                           cAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTAc                                                                                                        >  2:157635/1‑149 (MQ=255)
                                            aGGGATTGTCTGCTCGCACCACCATTTTTATTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACa                                                                                                       <  8:192584/149‑1 (MQ=255)
                                             gggATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAg                                                                                                      >  2:147055/1‑149 (MQ=255)
                                              ggATTGTCTGCTCGCACAACCAGTTTTACTAGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCACCGCTGAAAAATCTCGGGTGGCGGCTCGATCGCGCCCTTGAGGAAGAGTCGCTCTATTACa                                                                                                       >  8:10295/1‑147 (MQ=255)
                                                       gCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTa                                                                                                         >  8:377709/1‑136 (MQ=255)
                                                       gCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTa                                                                                                         <  7:377709/136‑1 (MQ=255)
                                                        cTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGAttt                                                                                           <  5:221213/149‑1 (MQ=255)
                                                         tCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGATCGATCGCTCCCTTGATGAAGAGTCGCTCTATTAAAGCGGACTGATTTg                                                                                          >  3:205656/1‑149 (MQ=255)
                                                           gCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCAGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGtt                                                                                        >  2:2996/1‑149 (MQ=255)
                                                             acaACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTcc                                                                                      >  4:160681/1‑149 (MQ=255)
                                                              caACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTccc                                                                                     <  1:157635/149‑1 (MQ=255)
                                                              caACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTccc                                                                                     <  2:98512/149‑1 (MQ=255)
                                                              caACCATTTTACCTCGGTCCTATCCCGCGCAATTATGGCCAACATCGCCCCTTGCCCGGCTTTCATCTCAACGCTGAAAACTCTCTGGTGGCTGCTCGATCGCGCCCTTGATGCAGAGTCCCTCTAATACCGCGGACTGATTTGTTccc                                                                                     <  2:199465/149‑1 (MQ=255)
                                                               aaCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCt                                                                                    >  8:173735/1‑149 (MQ=255)
                                                               aaCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCt                                                                                    >  5:319401/1‑149 (MQ=255)
                                                                aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTg                                                                                   <  3:199536/149‑1 (MQ=255)
                                                                   aTTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAACCATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTgg                                                                                  >  6:281277/1‑147 (MQ=255)
                                                                   aTTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTgg                                                                                  >  4:5632/1‑147 (MQ=255)
                                                                   aTTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGac                                                                                >  4:171721/1‑148 (MQ=255)
                                                                    tttttACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGc                                                                                                     >  4:180274/1‑127 (MQ=255)
                                                                    tttttACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGc                                                                                                     <  3:180274/127‑1 (MQ=255)
                                                                           tCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGAtt                                                                          >  7:231002/1‑147 (MQ=255)
                                                                           tCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCACTTGCCCGGATTTCATCTCAACGCAGAAAAATCTCTGGTGGCGCCTCGATCGCGCCCTTGATCAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGCAGCGATTaa                                                                        >  3:63722/1‑149 (MQ=255)
                                                                             ggggCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGTTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTaa                                                                        >  3:291149/1‑147 (MQ=255)
                                                                             ggggCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGTTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTaa                                                                        >  3:296075/1‑147 (MQ=255)
                                                                              gggCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAagag                                                                     >  6:7597/1‑149 (MQ=255)
                                                                               ggCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGc                                                                    <  5:281277/149‑1 (MQ=255)
                                                                                gCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCAGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCg                                                                   <  7:20913/149‑1 (MQ=255)
                                                                                 cTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCgg                                                                  <  4:295058/149‑1 (MQ=255)
                                                                                    tCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGAt                                                                <  3:160681/148‑1 (MQ=255)
                                                                                     cccGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATg                                                               <  6:319401/148‑1 (MQ=255)
                                                                                      ccGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTa                                                             >  7:325038/1‑149 (MQ=255)
                                                                                       cgcgGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATGAAGAGCGGATGTa                                                             >  4:256597/1‑148 (MQ=255)
                                                                                        gcgGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCt                                                                                    >  1:283474/1‑124 (MQ=255)
                                                                                        gcgGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCt                                                                                    <  2:283474/124‑1 (MQ=255)
                                                                                           gAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCg                                                        >  3:189578/1‑149 (MQ=255)
                                                                                           gAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCg                                                        <  6:230435/149‑1 (MQ=255)
                                                                                              ttATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGtt                                                     <  7:173735/149‑1 (MQ=255)
                                                                                                 tGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTAt                                                   <  4:242379/148‑1 (MQ=255)
                                                                                                 tGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATc                                                  >  4:86233/1‑149 (MQ=255)
                                                                                                 tGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCAGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTAt                                                   <  1:2996/148‑1 (MQ=255)
                                                                                                     aaaCATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATc                                              <  3:5632/149‑1 (MQ=255)
                                                                                                       aCATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCAc                                            >  3:37194/1‑149 (MQ=255)
                                                                                                         aTCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCa                                          >  6:219974/1‑149 (MQ=255)
                                                                                                              cccTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCct                                     >  1:368188/1‑149 (MQ=255)
                                                                                                                     ccGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGa                                                                 <  6:233574/114‑1 (MQ=255)
                                                                                                                     ccGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGa                                                                 >  5:233574/1‑114 (MQ=255)
                                                                                                                     ccGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGc                              <  3:227871/149‑1 (MQ=255)
                                                                                                                      cGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCa                                          >  5:217073/1‑136 (MQ=255)
                                                                                                                      cGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGACGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGCGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCg                             >  6:154456/1‑149 (MQ=255)
                                                                                                                      cGGATGTCATCTAAAAGATGAAAAATCTATGTTGGCGGCTCGATCGCGCCCTTGATGAATAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCa                                          <  6:217073/136‑1 (MQ=255)
                                                                                                                            gcATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTTATTTGTTACCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTatat                        <  4:205656/147‑1 (MQ=255)
                                                                                                                             cATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATcg                      >  5:176791/1‑149 (MQ=255)
                                                                                                                                  cAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGTGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCggg                 >  4:205182/1‑149 (MQ=255)
                                                                                                                                   aaCGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                <  5:7597/149‑1 (MQ=255)
                                                                                                                                   aaCGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                <  4:37194/149‑1 (MQ=255)
                                                                                                                                   aaCGCTGAAAAATCTCTGGTGGCGGCTCGATAGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATcggcgg                  <  4:185963/147‑1 (MQ=255)
                                                                                                                                       cTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGct            <  8:325038/149‑1 (MQ=255)
                                                                                                                                       cTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGct            <  3:86233/149‑1 (MQ=255)
                                                                                                                                          aaaaaTCTCTGGTTGCGGCTCGATCGCGCGCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTaa                                                                        <  3:67846/86‑1 (MQ=255)
                                                                                                                                           aaaaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGGtctc         >  5:311470/1‑148 (MQ=255)
                                                                                                                                                 tctGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTg  >  7:111447/1‑149 (MQ=255)
                                                                                                                                                   tggtggCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACAtt   <  1:147055/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCGACGCCAGCTTCAAAGAGATGCATGGCCGGATTGTGATGCCAGGGATTGTCTGCTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTG  >  NC_000913/3019284‑3019577

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