Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,559 A→G 5.5% intergenic (‑382/+15) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,5590AG 5.5% 214.1 / 3.2 73intergenic (‑382/+15)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base G (2/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGC  >  NC_000913/699410‑699698
                                                                                                                                                     |                                                                                                                                           
gTGCGAGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCa                                                                                                                                             <  2:1232524/150‑1 (MQ=255)
   cGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGAc                                                                                                                                          <  2:302225/150‑1 (MQ=255)
     aTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGGAGACGt                                                                                                                                        >  1:754812/1‑150 (MQ=255)
        tgtgCATGCCGGATGCGGCGTGTACGCCTTATCCGGCGTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGt                                                                                                                                     <  2:1091363/150‑1 (MQ=21)
         gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGCCGTGGtt                                                                                                                                    >  2:392594/1‑150 (MQ=255)
         gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGtt                                                                                                                                    <  1:977384/150‑1 (MQ=255)
         gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGtt                                                                                                                                    >  1:1105610/1‑150 (MQ=255)
         gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGtt                                                                                                                                    <  1:877246/150‑1 (MQ=255)
         gTTAATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGtt                                                                                                                                    >  2:608703/1‑150 (MQ=255)
          ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGTGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTg                                                                                                                                   >  1:1080040/1‑150 (MQ=255)
           tcctGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTgg                                                                                                                                  <  2:849519/147‑2 (MQ=255)
           tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCAGGAAAAATCAATAAATTGCAGGAGCTGCGTAGGCCGGATAAGCGTAGCGCATCAGGCAGTTTGGCGGTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGc                                                                                                                                  <  2:56582/150‑1 (MQ=21)
                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCTTCAGAGCCAACCACGTCCGCAGACGGGGTTGCTATTc                                                                                                                             >  1:1203017/1‑150 (MQ=255)
                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                             >  2:602378/1‑150 (MQ=255)
                  ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAg                                                                                                                           <  2:745745/150‑1 (MQ=255)
                       cggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAAc                                                                                                                      >  1:750668/1‑150 (MQ=255)
                                cGCCTTATCCGGCCTACAACAGCCTGAAACTTCAATATAGTGCAGGAGCGGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATGTc                                                                                                             <  1:507046/150‑1 (MQ=25)
                                  ccTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                           >  1:98714/1‑150 (MQ=255)
                                  ccTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                           >  1:1070056/1‑150 (MQ=255)
                                   cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGc                                                                                                          >  2:995569/1‑150 (MQ=255)
                                   cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGc                                                                                                          <  1:1062302/150‑1 (MQ=255)
                                           gCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACg                                                                                                  <  1:645212/150‑1 (MQ=255)
                                           gCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACg                                                                                                  >  2:771811/1‑150 (MQ=255)
                                               aCAATAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGGGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                                                              <  2:1011978/150‑1 (MQ=255)
                                               aCAAAAGCATGTAAATTCAATATATTGCTGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGCCGTGGTTGCTATTCAGATAACGTCGATTTCCGCGACTGCCCGGtt                                                                                              >  1:999882/1‑149 (MQ=37)
                                               aCAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                                                              <  1:1102001/150‑1 (MQ=255)
                                               aCAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                                                              <  1:1111561/150‑1 (MQ=255)
                                                   aaGCATGAAAAATCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATc                                                                                          <  1:481885/150‑1 (MQ=255)
                                                    aGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATcc                                                                                         <  1:513650/150‑1 (MQ=255)
                                                      cATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGGCGGGTAAATCCAg                                                                                       >  2:834167/1‑150 (MQ=255)
                                                       aTGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGc                                                                                      <  2:1105610/150‑1 (MQ=255)
                                                           aaaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTggg                                                                                  <  2:1314342/150‑1 (MQ=255)
                                                            aaatCAATATATTGCAGGAGCGTTGTAGGCCTGATAAGCGTAGCTCATCAGGCAGGTTGTCGTTTGTAATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTgggg                                                                                 <  2:255644/147‑1 (MQ=37)
                                                             aTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGCCGTGGTTGCTATTCAGATAACGTCGATTTCAGCTACTGACGGGTAAATCCAGCTGGGGc                                                                                >  1:971815/1‑150 (MQ=255)
                                                             aTTCAAAAAATTGCAGCACCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTGGCCGTTTGGCATCAGAGCCAACCACGCCCGCAGACGTGGTTGCTATCCAGATAACGCCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGc                                                                                <  2:1093414/150‑1 (MQ=38)
                                                              ttCAATATATTGCAGGAGCTGCGTAGGCATGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCg                                                                               <  1:758272/150‑1 (MQ=255)
                                                               tCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCTGCGACTGCCGGGTAAATCCAGCTGGGGCgg                                                                              >  2:401690/1‑150 (MQ=255)
                                                                  atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGCCTGCCGGGTAAATCCAGCTGGGGCGGaaa                                                                           >  2:984873/1‑150 (MQ=255)
                                                                    atatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                         >  2:1192810/1‑150 (MQ=255)
                                                                    atatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                         <  1:587643/150‑1 (MQ=255)
                                                                     tatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGc                                                                        >  1:769216/1‑150 (MQ=255)
                                                                        tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATa                                                                     >  2:222160/1‑150 (MQ=255)
                                                                         gCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATAc                                                                    <  1:185563/150‑1 (MQ=255)
                                                                           aGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGGTGGGGCGGAAAGGCATACTg                                                                  >  1:795633/1‑150 (MQ=255)
                                                                            ggAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGt                                                                 >  1:867530/1‑150 (MQ=255)
                                                                                cTGCGTATGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCATCCACTTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGat                                                             <  1:188393/150‑1 (MQ=255)
                                                                                     tAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgt                                                        >  1:474733/1‑150 (MQ=255)
                                                                                      aGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtc                                                       <  1:608703/150‑1 (MQ=255)
                                                                                           tGATAAGGGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCATCCACGTCCGCAGACGTGGTTGCTATTCAGATACCGTCGTTTTGAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTAGAGcg                                                  >  2:676513/1‑150 (MQ=39)
                                                                                           tGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGGTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCCAGcg                                                  >  1:856191/1‑150 (MQ=255)
                                                                                                aGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGCCGaa                                             >  2:191837/1‑150 (MQ=255)
                                                                                                 gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGaaa                                            >  2:1030561/1‑150 (MQ=255)
                                                                                                  cGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                           <  1:635600/150‑1 (MQ=255)
                                                                                                  cGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                           <  2:361586/150‑1 (MQ=255)
                                                                                                         cATCAGGCAGTTTGGCGGTTGTCATCAGAGCCAACCACGTACGCAGACGTTGTTGCTATTCAGATAACGTCGCTTTCAGCGCCGGACGGGTAAATACAGCTGGCGCGTATAGGCATACTCTCGATATCGTAGAGCGACGAAACACGAGaa                                    >  2:251331/1‑150 (MQ=255)
                                                                                                                  gTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGa                           >  1:112340/1‑150 (MQ=255)
                                                                                                                        cGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGAGGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTAGATATCGTCGCGCGGCGAAACCCCAGAAAGCACCAGAATCGtg                     >  1:1161855/1‑149 (MQ=255)
                                                                                                                        cGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGCAAGCACCAGACTCGtc                     >  1:704150/1‑150 (MQ=255)
                                                                                                                           ttGTCATCAGCGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCTTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTcc                  <  2:384654/150‑1 (MQ=255)
                                                                                                                            tGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTccc                 >  1:672428/1‑149 (MQ=255)
                                                                                                                            tGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCAAACTGTCGATATCGTCGAGCGCCGAAACAGCAGAAAGCACCAGAATCGTCTccc                 >  2:1235489/1‑149 (MQ=255)
                                                                                                                             gTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTTTTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGATAGGCATACTGTCGATATCGTCGAGCGCCGAAACACCAGCAAGCACCAGAATCGTCGCCAg                >  2:800040/1‑150 (MQ=255)
                                                                                                                              tcatcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTAAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCAGACTGTCGATATCGTCGAGCGCCGAAACACCAGAAAGCCCCAGAATCGTCGCCAgc               >  1:1180869/1‑149 (MQ=255)
                                                                                                                               catcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAc              <  2:1209142/150‑1 (MQ=255)
                                                                                                                                atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcc             <  1:222160/150‑1 (MQ=255)
                                                                                                                                    gagCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGCCGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcctgcc         >  1:30852/1‑150 (MQ=255)
                                                                                                                                    gagCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGCAACACCAGAAAGCACCAGAATCGTCTCCAGAcctgcc         >  2:870003/1‑150 (MQ=255)
                                                                                                                                    gagCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcctgcc         <  1:357674/150‑1 (MQ=255)
                                                                                                                                      gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcctgcctg       <  2:750668/150‑1 (MQ=255)
                                                                                                                                       ccAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTgg      >  1:139821/1‑150 (MQ=255)
                                                                                                                                        cAACCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGa     >  2:1008468/1‑150 (MQ=255)
                                                                                                                                         aaCCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGac    >  1:1222294/1‑149 (MQ=255)
                                                                                                                                          aCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAg   >  1:1053251/1‑150 (MQ=255)
                                                                                                                                           ccACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGGCGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGc  <  1:1268096/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGC  >  NC_000913/699410‑699698

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