Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,880 A→C 5.5% intergenic (+97/+55) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8800AC 5.5% 890.2 / 5.1 307intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (115/158);  new base C (19/4);  total (137/170)
Fisher's exact test for biased strand distribution p-value = 2.73e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504731‑3505029
                                                                                                                                                     |                                                                                                                                                      
 cttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGACCACACCCAAAACCAAACCCCaa                                                                                                                                                       >  1:3097307/1‑150 (MQ=38)
 cttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaa                                                                                                                                                       <  1:4227767/150‑1 (MQ=255)
       ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACACCCCAAAGGGGc                                                                                                                                                 >  2:466035/1‑150 (MQ=255)
         gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                               <  1:5534449/150‑1 (MQ=255)
           gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAgg                                                                                                                                             <  1:3627101/150‑1 (MQ=255)
           gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                             <  1:4636423/150‑1 (MQ=255)
            ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAggg                                                                                                                                            <  1:7223597/150‑1 (MQ=39)
             cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACCACTACCCCAAGGGGGGAgggg                                                                                                                                           >  1:4974904/1‑150 (MQ=38)
              aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                                                          <  2:1286599/150‑1 (MQ=255)
                tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑cc                                                                                                                                       <  2:4570511/150‑1 (MQ=255)
                tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCCCACCCGAAACCACAACCCAAAAGGGCCAGGGGc‑cc                                                                                                                                       >  2:4961078/1‑150 (MQ=38)
                    aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTc                                                                                                                                   <  1:5439652/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTa                                                                                                                               <  2:4628358/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTa                                                                                                                               <  1:5659219/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTa                                                                                                                               >  2:5886799/1‑150 (MQ=255)
                          ttAACTATTGATAATACAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGCCCCGGACCCTTACCCCCTCCCCAAAGGGGCGGGGGGC‑CCGTTCCGGACg                                                                                                                             >  2:5059547/1‑150 (MQ=38)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGTCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTTCCCGCACACAAAAACCCACCCCCAACGGGCGAGGGGA‑CAGCTCATTAcc                                                                                                                             >  2:2481416/1‑149 (MQ=37)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCACACCCCAAAGGGGCGAGGGGG‑CCGGTAAGTACg                                                                                                                             >  2:6687863/1‑150 (MQ=38)
                            aaCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTAAATCCCGTACGTTACCCACACCCAAACCCACACCCCAAACGGGGGCGGGGG‑CCGGTCAGTAAGGa                                                                                                                           >  1:3414453/1‑150 (MQ=37)
                             aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGaa                                                                                                                          <  1:6936707/150‑1 (MQ=255)
                               tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCGAACCCACAACCCAAAGGGGCGAGGGGG‑CCGGTAAGTAAgggaat                                                                                                                        >  1:5770614/1‑146 (MQ=37)
                               tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAAtt                                                                                                                        <  1:4560154/150‑1 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CTGTTCAGTACGGAATTTa                                                                                                                      <  1:5039704/150‑1 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTa                                                                                                                      <  2:5537982/150‑1 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTa                                                                                                                      <  1:2747703/150‑1 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTa                                                                                                                      <  2:6126097/150‑1 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTa                                                                                                                      <  2:1221020/150‑1 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACTCCCCCAAGGGGCGAGGGGC‑CCGTTCAGTTCGGGATTTa                                                                                                                      >  1:7175792/1‑150 (MQ=38)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCAAACCCAAACCCACAACCCCAACGGGCGCGGGGA‑CCGGTCAGTAAGGCATTTa                                                                                                                      >  1:3483017/1‑150 (MQ=37)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCCATGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CTGTTCAGTACGGAATTTa                                                                                                                      <  1:5987566/150‑1 (MQ=255)
                                    ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAACCCACTCCCCCAACGGGCGAGGGGG‑CCCGTCAGTACGGCATTTATCt                                                                                                                   >  1:3806155/1‑150 (MQ=37)
                                         aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACg                                                                                                              <  1:6750616/150‑1 (MQ=255)
                                          aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACgg                                                                                                             <  1:1735675/150‑1 (MQ=255)
                                           aaGACCTTAACCGTAGGCCGGCTAAGATGCGCCAGCATCGCATCCGGCGATGCTCGCGCGTTGAATTGTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                                            <  2:3962773/150‑1 (MQ=39)
                                            aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  1:6194102/150‑1 (MQ=255)
                                            aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGGTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           >  2:6020155/1‑150 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCACCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACAACCCGTACGTTGCCCTCCCCCTAACCCACTCCCCCCAGGGGCGAGGGGC‑CCCGTCAGGACCGCAGTTAGCGTGCCGGGGCtctc                                                                                                      >  2:4280666/1‑150 (MQ=25)
                                                  tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                                     <  1:4944121/150‑1 (MQ=255)
                                                  tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCACACCCAAACCCACAACCCCAACGGGCGGGGGGACCCGTACGTAACGGAA‑TTAGCATTGCGGTGCTAtac                                                                                                     >  1:4786795/1‑148 (MQ=19)
                                                    aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                                   <  2:3436312/150‑1 (MQ=255)
                                                         aGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTa                                                                                              <  2:6236711/150‑1 (MQ=255)
                                                           gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCCAACCCACAACCCCAAGGGGCGGGGGGA‑CCGGTCAGTGAGGAATTTTGCGTTAACGTGCTATACAggcataaa                                                                                            >  2:5557593/1‑143 (MQ=255)
                                                           gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATCCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                            <  2:1898422/150‑1 (MQ=255)
                                                              ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAc                                                                                         <  2:5036601/150‑1 (MQ=255)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                        <  1:6356852/150‑1 (MQ=255)
                                                                aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  1:1640660/150‑1 (MQ=255)
                                                                aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  1:61683/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGACTTTTACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:1367709/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:6100861/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:485901/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:3510358/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:1384859/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:1912830/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGGTCCCCTCACCCTAACCCGCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:2320487/150‑1 (MQ=255)
                                                                  ttccatcgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGc                                                                                     <  2:6861332/144‑1 (MQ=255)
                                                                    gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCAAACCCTAAAACACAACCCAAAAGGGCGAGGGGG‑CCGTTAAGTACCGAATTTAGCTTGGAGGTTCTATAAAGCATTAAAAACAGGCaa                                                                                   >  2:6210415/1‑150 (MQ=255)
                                                                      tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                                                 <  1:701671/150‑1 (MQ=255)
                                                                         gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                                              <  2:220490/150‑1 (MQ=255)
                                                                               atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCTCTCCCCACAGGGGCGAGGGGA‑CCGGTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCACTACTCTCACGc                                                                        >  2:884745/1‑150 (MQ=255)
                                                                                      ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                                                 <  1:4671044/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               >  2:4823404/1‑150 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTGAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  1:15930/150‑1 (MQ=255)
                                                                                             tGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAc                                                          <  2:200585/150‑1 (MQ=255)
                                                                                              gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                                                         <  1:229623/150‑1 (MQ=255)
                                                                                                  gcgcGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtct                                                     <  2:5702520/150‑1 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  2:2993867/150‑1 (MQ=255)
                                                                                                          aaTTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCg                                             <  2:814764/150‑1 (MQ=255)
                                                                                                            ttttACATCCCGTACGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACCGAAATTAGCTTGACCGTGCTAACCAGCATTAAAAACTGGCAATACACTCACGCAAAATCCGCACCACCCTCACATCAccccg                                           >  2:5287672/1‑147 (MQ=37)
                                                                                                             tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  2:3501812/150‑1 (MQ=255)
                                                                                                                     ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  <  2:6515076/150‑1 (MQ=255)
                                                                                                                       gTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                                <  1:930421/150‑1 (MQ=255)
                                                                                                                            ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                           <  2:6502071/150‑1 (MQ=255)
                                                                                                                             tCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                          <  2:1616055/150‑1 (MQ=255)
                                                                                                                              ccccTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCCCGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCg                         <  2:5619397/150‑1 (MQ=255)
                                                                                                                               cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa                        <  2:2915609/150‑1 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCCACGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       <  1:6020155/150‑1 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       <  1:6605962/150‑1 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:3570648/1‑150 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:3738202/1‑150 (MQ=255)
                                                                                                                                 cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGGAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                      >  1:2399131/1‑150 (MQ=255)
                                                                                                                                    ccccTAACCCTCTCCCCAACGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                   <  1:6741913/149‑1 (MQ=255)
                                                                                                                                    accctcaccctcTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCGCGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                   <  1:840600/150‑1 (MQ=255)
                                                                                                                                     cccTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                  <  2:642019/150‑1 (MQ=255)
                                                                                                                                     cccTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGca                  <  1:868553/150‑1 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:4424079/1‑150 (MQ=255)
                                                                                                                                       cTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  1:108884/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa               >  1:3865273/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCACCCAAAAGGGGCGAGCCGA‑CAGTTCATTACGGAATTTAGCTTGACGGAGCTATCCAGCATAAAAAACAGGCAATACTCTCACGCAAAAACCGCAGCACCTTCACTTCACCGCTGCGATACGGGTAAAACCGAATCGCACaa               >  2:5431972/1‑150 (MQ=255)
                                                                                                                                         tccccTCGCCCCTTTGGGGAGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAAt              >  2:7090274/3‑150 (MQ=255)
                                                                                                                                         aaCCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAAt              <  2:3260191/150‑1 (MQ=255)
                                                                                                                                           cccTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  1:4096977/1‑150 (MQ=255)
                                                                                                                                           cccTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:3353881/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           <  1:4418404/150‑1 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  1:6643924/150‑1 (MQ=255)
                                                                                                                                              tctcCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGGTAATGCGAATCGCACAATGTTCt         >  2:2343113/1‑150 (MQ=255)
                                                                                                                                                tcCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGATGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGa       >  1:2926786/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:5969733/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:4730037/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:5240724/1‑150 (MQ=255)
                                                                                                                                                     aaaGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  >  2:2080368/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGA‑CCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504731‑3505029

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑