Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,305,914 G→T 8.3% intergenic (‑97/+57) mtr ← / ← deaD tryptophan:H(+) symporter Mtr/ATP‑dependent RNA helicase DeaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,305,9140GT 8.3% 317.0 / 12.4 97intergenic (‑97/+57)mtr/deaDtryptophan:H(+) symporter Mtr/ATP‑dependent RNA helicase DeaD
Reads supporting (aligned to +/- strand):  ref base G (47/41);  new base T (2/6);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

CCACGCCGCCAAGCAGCGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCAT  >  NC_000913/3305766‑3306062
                                                                                                                                                    |                                                                                                                                                    
ccacGCCGCCAAGCAGCGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGt                                                                                                                                                     >  2:1102401/1‑150 (MQ=255)
ccacGCCGCCAAGCAGCGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGt                                                                                                                                                     <  2:923102/150‑1 (MQ=255)
       gccAAGCAGCGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATc                                                                                                                                              <  1:221099/150‑1 (MQ=255)
            gcagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTAGACGAGTACAAAAGACAGAAAAAAAGCCCCGATTTTTACCATCGGGGc                                                                                                                                         <  2:1270491/150‑1 (MQ=255)
            gcagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATTTTTACCATCGGGGc                                                                                                                                         <  2:205902/150‑1 (MQ=255)
            gcagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATTTTTACCATCGGGGc                                                                                                                                         <  2:8379/150‑1 (MQ=255)
            gcagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGc                                                                                                                                         <  1:696001/150‑1 (MQ=255)
             cagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATTTTTACCATCGGGGCt                                                                                                                                        <  2:1216663/150‑1 (MQ=255)
             cagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCgg                                                                                                                                            <  2:1102336/146‑1 (MQ=255)
             cagcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCgg                                                                                                                                            >  1:1102336/1‑146 (MQ=255)
              agcGACGGTGTCGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCTGTTAAAAATATTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATTTTTAGCATCGGGGCtt                                                                                                                                       <  2:1023560/150‑2 (MQ=255)
               gcGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGt                                                                                                                                      <  2:349534/150‑1 (MQ=255)
                cGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGta                                                                                                                                     >  2:1332432/1‑150 (MQ=255)
                 gACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGtat                                                                                                                                    >  2:324111/1‑150 (MQ=255)
                 gACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGtat                                                                                                                                    <  2:1353587/150‑1 (MQ=255)
                   cGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGTCAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGtatat                                                                                                                                  <  2:1309656/150‑1 (MQ=255)
                    ggTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATTTTTACCATCGGGGCtttttttc                                                                                                                                 <  1:1481350/150‑8 (MQ=255)
                              gggTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACa                                                                                                                       <  1:1225446/150‑1 (MQ=255)
                              gggTGGGGGTTAGTGTTGCCTTGACTGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTTTGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGCTAAAAATCGGGGCTGTATATATTATTTTACa                                                                                                                       <  2:1255510/150‑1 (MQ=255)
                               ggtggtggtTAGTGGTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAACCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAg                                                                                                                      <  1:735864/150‑1 (MQ=255)
                                        tAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTc                                                                                                             <  2:649968/150‑1 (MQ=255)
                                        tAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTc                                                                                                             <  1:1102401/150‑1 (MQ=255)
                                          gtgtTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGc                                                                                                           <  2:969237/150‑1 (MQ=255)
                                          gtgtTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGc                                                                                                           <  1:1332432/150‑1 (MQ=255)
                                          gtgtTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGc                                                                                                           >  2:591531/1‑150 (MQ=255)
                                          gtgtTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGc                                                                                                           >  1:1104774/1‑150 (MQ=255)
                                             ttGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGt                                                                                                        <  1:694777/150‑1 (MQ=255)
                                              tGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGtt                                                                                                       <  2:1104774/150‑1 (MQ=255)
                                               gCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGGAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTc                                                                                                      >  2:641935/1‑150 (MQ=255)
                                               gCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATtgtg                                                                                                                <  2:796090/140‑1 (MQ=255)
                                               gCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATtgtg                                                                                                                >  1:796090/1‑140 (MQ=255)
                                                 cATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAg                                                                                                    <  1:293895/150‑1 (MQ=255)
                                                   tGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCg                                                                                                  >  1:1436649/1‑150 (MQ=255)
                                                    gAGGGCTTGTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCAGTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCga                                                                                                 <  2:1515859/150‑1 (MQ=255)
                                                      gggCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCgatg                                                                                               <  2:1083296/150‑1 (MQ=255)
                                                       ggCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTAtttt                                                                                                                          >  1:1485873/1‑122 (MQ=255)
                                                       ggCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTAtttt                                                                                                                          <  2:1485873/122‑1 (MQ=255)
                                                       ggCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCgatga                                                                                              <  1:544439/150‑1 (MQ=255)
                                                         cTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGAtt                                                                                            >  2:341865/1‑150 (MQ=255)
                                                            ctctcCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACg                                                                                         <  1:558934/150‑1 (MQ=255)
                                                                   gTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATcacc                                                                                  >  1:908665/1‑150 (MQ=255)
                                                                   gTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATcacc                                                                                  >  2:377882/1‑150 (MQ=255)
                                                                    tGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGAGATGATTACGCATcacca                                                                                 >  2:893512/1‑150 (MQ=255)
                                                                     gAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATcaccac                                                                                >  1:948805/1‑150 (MQ=255)
                                                                              tgtgACTGCGTTGTTTTGCATTGCTCTGTACCTTTACACGAGGACATAATACAGAAATAAAGCCCCGATTGTAAAAATCGGGGCGGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCTATGATTACGCATCACCACCGAAacgac                                                                       <  2:1523018/147‑1 (MQ=255)
                                                                               gCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATCGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATAACCACCGAAacgacg                                                                      >  2:278859/1‑150 (MQ=255)
                                                                                cGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCAATGATTACGc                                                                                        >  1:48744/1‑131 (MQ=255)
                                                                                cGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCAATGATTACGc                                                                                        <  2:48744/131‑1 (MQ=255)
                                                                                  aCTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAacgacgacg                                                                   <  1:1326952/150‑1 (MQ=255)
                                                                                     gCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATcacca                                                                                 >  1:794988/1‑133 (MQ=255)
                                                                                     gCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATcacca                                                                                 <  2:794988/133‑1 (MQ=255)
                                                                                           ttATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGa                                                          >  2:1053310/1‑150 (MQ=255)
                                                                                            tATGCATTGCACTGTACCAGTACACGAGTACAAAAGTCAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGaa                                                         >  2:1394726/1‑150 (MQ=255)
                                                                                            tATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGaa                                                         >  2:850781/1‑150 (MQ=255)
                                                                                               gCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCACGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAatca                                                      >  2:688351/1‑150 (MQ=255)
                                                                                                   tGCACTGTACCAGTACACTAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCac                                                  <  1:351537/150‑1 (MQ=255)
                                                                                                    gCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCacg                                                 >  1:168565/1‑150 (MQ=255)
                                                                                                      aCTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCacgac                                               >  1:607514/1‑150 (MQ=255)
                                                                                                        tGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGAcgc                                             <  1:1394726/150‑1 (MQ=255)
                                                                                                        tGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGAcgc                                             <  1:1078/150‑1 (MQ=255)
                                                                                                         gTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCgg                                                              >  1:730463/1‑132 (MQ=255)
                                                                                                         gTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCgg                                                              <  2:730463/132‑1 (MQ=255)
                                                                                                         gTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGAcgcg                                            >  2:231381/1‑150 (MQ=255)
                                                                                                         gTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGAcgcg                                            >  2:612280/1‑150 (MQ=255)
                                                                                                             cAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAgc                                        >  2:1234731/1‑150 (MQ=255)
                                                                                                                 acacGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCgg                                    >  1:29938/1‑150 (MQ=255)
                                                                                                                 acacGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTAAGCATAACAACAGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCgg                                    >  2:347310/1‑150 (MQ=255)
                                                                                                                  cacGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGc                                   >  2:363396/1‑150 (MQ=255)
                                                                                                                  cacGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGc                                   >  2:1307514/1‑150 (MQ=255)
                                                                                                                   acGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGcc                                  <  2:534726/150‑1 (MQ=255)
                                                                                                                   acGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGcc                                  >  2:665291/1‑150 (MQ=255)
                                                                                                                    cGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCt                                 >  2:651395/1‑150 (MQ=255)
                                                                                                                    cGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCt                                 <  2:1311019/150‑1 (MQ=255)
                                                                                                                    cGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCt                                 >  2:962665/1‑150 (MQ=255)
                                                                                                                           aaaaGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGAGATGATTACGCATAAACACCGAAACGACGACGACCGGTAGAATCATCACGACCCCGAGCGCGGCCTTCacgac                          >  2:925287/1‑150 (MQ=255)
                                                                                                                             aaGACAGGAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGt                        <  1:641935/150‑1 (MQ=255)
                                                                                                                               gACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTc                      >  2:558592/1‑150 (MQ=255)
                                                                                                                                 cAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgc                    >  1:898823/1‑150 (MQ=255)
                                                                                                                                 cAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgc                    <  1:893512/150‑1 (MQ=255)
                                                                                                                                 cAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgc                    >  1:1404454/1‑150 (MQ=255)
                                                                                                                                 cAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgc                    >  1:787265/1‑150 (MQ=255)
                                                                                                                                 cAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgc                    >  1:606482/1‑150 (MQ=255)
                                                                                                                                   gAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAacgacg                                                                      <  2:1504007/98‑1 (MQ=255)
                                                                                                                                   gAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAacgacg                                                                      >  1:1504007/1‑98 (MQ=255)
                                                                                                                                   gAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgccgc                 >  1:551765/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGCCCCGATGGTAAAAATCGTGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgccgc                 <  2:633835/150‑1 (MQ=255)
                                                                                                                                    aaaaaaaGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgccgc                 >  1:964727/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTcgccgc                 >  1:356921/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCt                <  2:908665/150‑1 (MQ=255)
                                                                                                                                        acaGCCCCGATTTTTACCATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACAGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTaaa             >  1:462067/3‑150 (MQ=255)
                                                                                                                                         caGCCCCGATTTTTACCATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTaaaa            >  2:1053904/2‑150 (MQ=255)
                                                                                                                                         aaGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTaaaa            <  1:1096246/150‑1 (MQ=255)
                                                                                                                                          aGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAac           >  2:1315380/1‑150 (MQ=255)
                                                                                                                                           gCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAacg          >  1:132953/1‑150 (MQ=255)
                                                                                                                                           gCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAacg          <  2:1190870/150‑1 (MQ=255)
                                                                                                                                             cccGATGGTAAAAATCGGGGCTTTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAacgac        <  1:1284857/150‑1 (MQ=255)
                                                                                                                                               cGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAAcgacga      <  1:1315380/150‑1 (MQ=255)
                                                                                                                                                   ggTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCAt  <  2:132953/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCACGCCGCCAAGCAGCGACGGTGACGTTTGGGTGGTGGTTAGTGTTGCCATGAGGGCTTCTCTCCAGTGAAAAATAGTGCGACTGCGTTGTTATGCATTGCACTGTACCAGTACACGAGTACAAAAGACAGAAAAAAAGCCCCGATGGTAAAAATCGGGGCTGTATATATTATTTTACAGATTGTGTTCGCTGTTCAGCGATGATTACGCATCACCACCGAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCAT  >  NC_000913/3305766‑3306062

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