Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,880 A→C 7.7% 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 7.7% 739.4 / 13.2 266intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (85/145);  new base C (21/3);  total (110/156)
Fisher's exact test for biased strand distribution p-value = 2.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504731‑3505029
                                                                                                                                                      |                                                                                                                                                     
  ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCaaa                                                                                                                                                     <  2:5043758/150‑1 (MQ=255)
  ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGCCCCACACCCCAAACCGC‑TCCCCaaa                                                                                                                                                     >  2:4536064/1‑150 (MQ=39)
   tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAg                                                                                                                                                    <  2:3347277/150‑1 (MQ=255)
      gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGCTGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCGC‑ACCCCCAAgggg                                                                                                                                                 >  2:2129067/1‑150 (MQ=255)
       ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAACCCAC‑ACCCCCAAggggg                                                                                                                                                >  1:1526787/1‑149 (MQ=38)
       ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGc                                                                                                                                                <  2:680432/150‑1 (MQ=255)
         gcgcCAATCGCATCACCTTAACTATTGATAATAAACGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGa                                                                                                                                              <  2:3705295/150‑1 (MQ=255)
          cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCCTCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAg                                                                                                                                             <  2:4043012/150‑1 (MQ=255)
          cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAg                                                                                                                                             <  2:4895258/150‑1 (MQ=255)
            ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑GCCCCTTTGGGGAGAggg                                                                                                                                           <  1:1325643/150‑1 (MQ=39)
             cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCCCCCGAACCCAA‑ACCCCCAAGGGGGGCgggg                                                                                                                                          >  1:5618390/1‑150 (MQ=38)
             cAATCGCATCACCTTAACTATAGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAACCCCA‑AACCCCAACGGGCGCgggg                                                                                                                                          >  1:1023411/1‑150 (MQ=37)
              aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGa                                                                                                                                         <  1:3289539/150‑1 (MQ=255)
               ctCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGGCCCCTCACCCTCACCCTC‑TCCCCAACGGGGCGAGGGGAc                                                                                                                                        <  1:971414/149‑1 (MQ=39)
                tCGCATCAACTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑GCCCCTTTGGGGAGAgggttag                                                                                                                                       <  1:114381/150‑5 (MQ=38)
                   cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGtt                                                                                                                                    <  1:3564088/150‑1 (MQ=255)
                     tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCa                                                                                                                                  <  1:1561201/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                               <  1:2797985/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                               <  1:2694468/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                              <  1:3940024/150‑1 (MQ=255)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCGAAAACAC‑AACCCAAAAGGGCGAGGGGGCCCTTCAGTaag                                                                                                                             >  1:558965/1‑148 (MQ=37)
                               tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                                                        <  1:2999522/150‑1 (MQ=255)
                                   gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                                    <  1:4696576/150‑1 (MQ=255)
                                    ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCCAACCCCC‑ACCCCAAACGGGCGAGGGGACCCGTCAGTTAGGAATTTAGCt                                                                                                                   >  2:1677015/1‑150 (MQ=38)
                                     taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCtt                                                                                                                  <  1:4378355/150‑1 (MQ=255)
                                       ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGa                                                                                                                <  2:4190035/150‑1 (MQ=255)
                                        taaAAGACATTAACCGTAGGCCGGATTAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                               <  2:61824/150‑1 (MQ=255)
                                        taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                               <  1:3999789/150‑1 (MQ=255)
                                          aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAAACCACAACCCCAAAGGGGC‑AGGGGGCCCTTAAGGACCGAAATTTGcatttacg                                                                                                             >  1:4090397/1‑143 (MQ=21)
                                          aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCGCACCCAAACCCCA‑ACCCCCAACGGGCGCGGGGCCCCGTCAGTACCGCAATTAGCTTGACgg                                                                                                             >  2:5507382/1‑150 (MQ=25)
                                           aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                                            <  2:5800287/150‑1 (MQ=255)
                                           aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                                            <  1:1898992/150‑1 (MQ=255)
                                              aCATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCt                                                                                                         <  1:2346341/150‑1 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      <  1:5053715/150‑1 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACACAAAACCAC‑TAACCAAACGGGCGAGGGGGCCCGTCAGTAACGAAATTAGCATGAAGGTGCAAta                                                                                                      >  2:2056546/1‑149 (MQ=25)
                                                   caCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcaa                                                                                                    <  1:3336597/149‑3 (MQ=255)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  2:1866576/150‑1 (MQ=255)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  2:1061740/150‑1 (MQ=255)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCAC‑ACCCCCAAGGGGCGCGGGGACCGTTCAGTACGGAAATTAGCATGACGGTGCTATCCa                                                                                                    >  1:5216705/1‑150 (MQ=38)
                                                      cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                 <  2:1151408/150‑1 (MQ=255)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCAC‑ACCCCAAAGGGGCGAGGGGGCCGTTCAGTACGGGAATTACCTTTACGGTGCGATCCAGCATTAAATACt                                                                                        >  2:3584163/1‑150 (MQ=255)
                                                                aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCAC‑AACCCCAAGGGGCGAGGGGCCCCGTCAGTTCGGAAATTAGCATGACGGGGCTATACAGCCTTAAATAcat                                                                                       >  2:3485508/1‑148 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:3374999/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:4290578/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:763385/150‑1 (MQ=255)
                                                                  aaGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGc                                                                                     <  1:2108326/150‑1 (MQ=255)
                                                                   aGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCGCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCa                                                                                    <  2:4731502/150‑1 (MQ=255)
                                                                    aacgcgcCAGCATCGCCGGATGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTGCCCGCACCCAAACACCC‑ACCCCCAAAGGGCGCGGGGGCCCGTGAGTTCGGACATTATCGTGACGGGGCGATACAGCCTTAAAAAAAGGCaa                                                                                   >  2:5569428/4‑150 (MQ=255)
                                                                      tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                                                 <  1:255848/150‑1 (MQ=255)
                                                                            aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                           <  1:4399071/150‑1 (MQ=255)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:4260403/150‑1 (MQ=255)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:3136362/150‑1 (MQ=255)
                                                                                  gcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaaa                                                                     <  2:5530525/150‑1 (MQ=255)
                                                                                   catcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGTCGATGGGACCGTTCAGTACTGAATTTAGCTTGACGTTGCTATCCAGCATTAAATACTGGCAATACTCTCACGAaaaa                                                                    <  2:2845061/150‑1 (MQ=39)
                                                                                     tcCGGCTATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTTACCCTC‑TCCCCTAAGGGGCGAGGGGACAGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATTCTGGCAATACTCTCACGCAAAATc                                                                  <  2:5259195/150‑1 (MQ=255)
                                                                                       cGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                                                <  1:4536064/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  2:2613722/150‑1 (MQ=255)
                                                                                         gCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgca                                                              <  1:5688094/150‑1 (MQ=255)
                                                                                         gCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgca                                                              <  2:3887223/150‑1 (MQ=255)
                                                                                         gCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACAGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgca                                                              <  2:1541410/150‑1 (MQ=255)
                                                                                            aTGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑ACCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATACgcagca                                                           >  2:4826450/1‑150 (MQ=255)
                                                                                              gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCTATACTCTCACGCAAAATCCGCAGCAcc                                                         <  2:12664/150‑1 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  2:2710229/150‑1 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  1:649259/150‑1 (MQ=255)
                                                                                                            ttttACATCCCGTACGTTCCCCTCACCCTAAACCAA‑AACCCAAATGGGCGAGGGGGCCGTTCAGTACGGAAATTAGCATTACGGTGCTATGCAGCAATAAATACAGGCAATAATGTCACGCAAAAACCGCAGCACCGTCTCATCAACGCt                                           >  2:3024258/1‑150 (MQ=255)
                                                                                                            ttttACATCCCGTACGTTACCCTCACCCTAACCCAC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTATCATGAAGGTGCAATACAGCAATAAATACAGGCAATACTCTCACGCAAAAAACGCACCACCCTCACTTCACCCCt                                           >  2:3143878/1‑150 (MQ=25)
                                                                                                                   tCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc                                    <  1:2109923/150‑1 (MQ=255)
                                                                                                                    cccGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                                   <  2:699399/150‑1 (MQ=255)
                                                                                                                    cccGTACGGTCCCCTCACCCTAACCCTC‑TCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                                   <  1:374839/150‑1 (MQ=255)
                                                                                                                     ccGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGTGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  <  1:5748512/150‑1 (MQ=255)
                                                                                                                     ccGTACGGTCCCCTCACCCTCACCCTC‑TCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  <  1:3623966/150‑1 (MQ=255)
                                                                                                                     ccGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  <  1:5507382/150‑1 (MQ=255)
                                                                                                                      cGTACGTTCCCCGCACCCTAACCCTC‑TACCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACAGGCAATACTCTCACGCAAAATCCGCACCACCGTCTCTTCACCGCTGCGATTCggg                                 >  2:171123/1‑150 (MQ=255)
                                                                                                                           gTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAt                            <  2:1450804/150‑1 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  2:1320761/1‑150 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:2200582/1‑150 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAT‑CAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGGTAATGCGAAt                      >  2:3983499/1‑150 (MQ=255)
                                                                                                                                   cACCCTAACCCTC‑TCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCg                    <  2:1989500/150‑1 (MQ=255)
                                                                                                                                      ccTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  1:2225418/1‑150 (MQ=255)
                                                                                                                                       ggtccccTC‑GCCCCTTTGGGGAGAGGGTACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  2:2181600/5‑150 (MQ=255)
                                                                                                                                        tAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa               >  2:2886723/1‑150 (MQ=255)
                                                                                                                                           cccTC‑GCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:512840/1‑150 (MQ=255)
                                                                                                                                            ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGACTCGCACACTGtt           >  2:3302240/1‑150 (MQ=255)
                                                                                                                                            ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  1:4216891/1‑150 (MQ=255)
                                                                                                                                            ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           <  1:347139/150‑1 (MQ=255)
                                                                                                                                            ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  1:5273296/1‑150 (MQ=255)
                                                                                                                                             ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  2:4301374/150‑1 (MQ=255)
                                                                                                                                             ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:4857871/1‑150 (MQ=255)
                                                                                                                                             ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  1:4867191/150‑1 (MQ=255)
                                                                                                                                             ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:1987334/1‑150 (MQ=255)
                                                                                                                                             ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:215701/1‑150 (MQ=255)
                                                                                                                                                 tcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGa       >  1:5619137/1‑150 (MQ=255)
                                                                                                                                                 tcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGa       >  1:4187559/1‑150 (MQ=255)
                                                                                                                                                   cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGagt     >  1:2453879/1‑148 (MQ=255)
                                                                                                                                                   cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  1:3101742/1‑150 (MQ=255)
                                                                                                                                                   cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     <  2:2414838/150‑1 (MQ=255)
                                                                                                                                                   cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:1752226/1‑150 (MQ=255)
                                                                                                                                                    ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:4603700/1‑150 (MQ=255)
                                                                                                                                                     cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   <  2:4392023/150‑1 (MQ=255)
                                                                                                                                                      aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  <  2:81189/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  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: ‑