Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,898 T→G 6.2% intergenic (+115/+37) 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,8980TG 6.2% 1802.2 / 8.2 556intergenic (+115/+37)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (201/298);  new base G (52/2);  total (255/301)
Fisher's exact test for biased strand distribution p-value = 4.87e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 566 total aligned reads displayed.

GCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACG  >  NC_000913/3504749‑3505047
                                                                                                                                                     |                                                                                                                                                     
gCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGt                                                                                                                                                       <  1:4383701/150‑1 (MQ=255)
  aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTc                                                                                                                                                     <  1:372807/150‑1 (MQ=255)
  aTCACATTAACTATTTATAATAAAAGACATTAACCGTCGGCCGTATTTGTTGCGACAGCATCGCATCCGTCTTTGATGGCGCGTTGATTTTTACTTCCGGTACGTTCCCCTCACCCTAACCATCTCCCCAAAGGGGCGAGGGGACCGTTc                                                                                                                                                     <  1:2422383/150‑1 (MQ=25)
   tCACCTTAACTATTGTTAATAAAAGACATTAACCGTAGGCCGGATTAGTTGCGCCAGCATCGCATCCGGCGATGCTGTCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCATCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                                    <  1:1226721/150‑1 (MQ=39)
   tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGTATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGAGTTGAATTTTACATCCCGTACGTTCCCCTCACACTAACCCTCTCCCCAAAGGTGCGAGGGGACCGTTCa                                                                                                                                                    <  1:345302/150‑1 (MQ=255)
   tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCAGCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                                    <  1:7131182/150‑1 (MQ=255)
        ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACg                                                                                                                                               <  2:1281526/150‑1 (MQ=255)
        ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCACTACCCCAAGGGGCGAGGGGACCGTTCAGTACg                                                                                                                                               >  1:6554034/1‑150 (MQ=255)
              aTTGATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAttt                                                                                                                                         <  1:8809951/150‑1 (MQ=255)
               ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCACATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTa                                                                                                                                        >  2:9138297/1‑150 (MQ=255)
                tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCCCCCCAACCCCCACCCCCAACGGGCGCGGGGGCCCGTCAGTTCGGGATTTAg                                                                                                                                       >  1:8578242/1‑150 (MQ=25)
                tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                                                       <  1:2313500/150‑1 (MQ=255)
                 gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCAAAACCAAAACCCAAAAGGGCGGGGGGCCCGGTGCGTAAGGAAATTAGc                                                                                                                                      >  1:1405480/1‑150 (MQ=25)
                  ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCt                                                                                                                                     <  1:5192721/150‑1 (MQ=255)
                    aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCCCACCCAAAACCACACCCCCAAGGGGCGCGGGGACCGGTCAGTAAGGAAATTAGCTTg                                                                                                                                   >  2:3787066/1‑150 (MQ=38)
                          aGACATTCACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                                             <  1:7523502/150‑1 (MQ=255)
                           gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGc                                                                                                                            <  2:1215364/150‑1 (MQ=255)
                                 aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                                      <  1:10790003/150‑1 (MQ=255)
                                 aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTGCCCACCCCCCAACCCACACCCCCAAGCGGCGAGGGGGCCCGTAACTACGGAATTTAGCATGACCGTGCTATACa                                                                                                                      >  1:6772925/1‑150 (MQ=25)
                                  aCCGTAGGCCGGATAAGACGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                                                     <  1:12125044/150‑1 (MQ=255)
                                    cggaGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                                   <  2:3286635/147‑1 (MQ=255)
                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                                   <  2:9535160/150‑1 (MQ=255)
                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                                   <  2:2171425/150‑1 (MQ=255)
                                         gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                                              <  2:4891616/150‑1 (MQ=255)
                                             gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGAACGTTGCCCGCACCCAAAACCACAACACAAAGGGGCGGGGGGACCCTTAAGTTCGGAAATTATATTGACGGTGCTATCCCGCCGTAAAAAca                                                                                                          >  2:3479712/1‑149 (MQ=255)
                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                                         <  2:10342969/150‑1 (MQ=255)
                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTAAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                                         <  2:2705048/150‑1 (MQ=255)
                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                                         <  2:1749409/150‑1 (MQ=255)
                                               ttcaacgcgcCAGCATCGCCGGATGCGATGCTGGCGCATTGAATTTTACATCCCGTACGGTGCCCGCACCCAAACCCGCTCCCCAAAGGGGCGCGGGGGCCGGTGAGTACGGGAATTATCGTGGCGGTGCGATCCACCATTAAATAATgg                                                                                                        >  2:2310198/7‑150 (MQ=16)
                                               tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                                        <  1:1890508/150‑1 (MQ=255)
                                               tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                                        <  1:4778759/150‑1 (MQ=255)
                                               tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                                        <  2:6963573/150‑1 (MQ=255)
                                                      cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                                                 <  2:3154054/150‑1 (MQ=255)
                                                       gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAGTTAGCTTGTCGGGGCTATCCAGCATTAAATACTGGCAATActc                                                                                                >  2:8224271/1‑150 (MQ=255)
                                                         cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                                              <  1:3873748/150‑1 (MQ=255)
                                                         cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                                              <  1:2681858/150‑1 (MQ=255)
                                                          aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                                             <  2:8741565/150‑1 (MQ=255)
                                                          aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                                             <  2:6355196/150‑1 (MQ=255)
                                                           gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                                            <  2:9189862/150‑1 (MQ=255)
                                                             atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGc                                                                                          <  2:4699692/150‑1 (MQ=255)
                                                                   tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATc                                                                                    >  2:7451938/1‑150 (MQ=255)
                                                                   tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAAGTTAGCTTGACGGTGCTATCCAGCCTTAAATACGGGCAATACTCACACGCCAAATc                                                                                    >  2:7611445/1‑150 (MQ=38)
                                                                   tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCAAACCCCCAAGGGGCGAGGGGACCGGTAAGTACGGGAATTTTCTTGACGGTGCGATCCACCATTAAATAATGGCAATAAACTAACGCAAAAta                                                                                    >  2:2618650/1‑149 (MQ=255)
                                                                      ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATCTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                                                 <  1:9436092/150‑1 (MQ=255)
                                                                      ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGGGAGGGGACCGTTAAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAATCTCAAGCAAAATCCgc                                                                                 >  2:10965872/1‑150 (MQ=255)
                                                                            gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                                                                           <  2:5255784/150‑1 (MQ=255)
                                                                              tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                                         <  2:7131136/150‑1 (MQ=255)
                                                                                gcgcGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtct                                                                       <  1:4098662/150‑1 (MQ=255)
                                                                                 cgcgTTGAATTGGACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCCCCGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAACATCCGCAGCACCGtctc                                                                      <  2:10191200/150‑1 (MQ=255)
                                                                                   cgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                                                                    <  2:89739/150‑1 (MQ=255)
                                                                                   cgTTGAATTTTACATCCCGTACGGTCCCCTCACCCTCACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTTGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                                                                    <  1:10712441/150‑1 (MQ=255)
                                                                                    gTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTc                                                                   <  1:10260675/150‑1 (MQ=255)
                                                                                     ttGAATTTTACATCCCGTACGTTCCCCTCACCCTAAACCTCTACCCAAAGGGGCGAGGGGACCGGTCCGTACGGAATTTTGCTTTACGGTGCTATCCAGCATTAAATAATGGCAATAATATCACGCAAAATCCGCAGCACCGTCTCTTCa                                                                  >  2:6553017/1‑150 (MQ=255)
                                                                                         aTTTTACATCCCGTACGTTCCCCTCACCCAAAACCACTCCCCCAAGGGGCGAGGGGACCGTTCAGTACCGAAATTAGCTTGAAGGTGCTATACACCCTTAAATAATGGCCAAAATCAAACCCAAAAAACGCACCACCCTCTAATCACCGc                                                              >  2:10747565/1‑150 (MQ=255)
                                                                                         aTTTTACATCCCGTACGTTACCCACACCCAAACCCTCTCCCCAAAGGGGCGAGGGGCCCCGTCAGTACGGAAATTAGGATTACGGTGCGAACCAGCATTAAAAACTGGCAAAAAACTCACGCAAAAACCGCCGCACCGTCTCTTCACCGc                                                              >  2:9389093/1‑150 (MQ=255)
                                                                                          ttttACATCCCGTACGTTCCCCTCACCCAAACCCACACCCCAAATGGGCGAGGGGACCGTTCAGTACGGTATTTAGCTTGACGGTGCTATCCACCAGTAAATACTGGCAATACTCGCGCGCAAAATCCGCAGCACCGTCTCTTCACCGcg                                                             >  2:6836023/1‑149 (MQ=255)
                                                                                           tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                                            <  2:5500327/150‑1 (MQ=255)
                                                                                             tACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAGGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCCTTAAATAAAGGCAATACTCTCACGCAAAAACCGCAGCACCGTCTCTTCACCGCTGCg                                                          >  2:10504667/1‑150 (MQ=38)
                                                                                              acctCCCGCACGTGCCCCTCCCCCGACCCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCGTGACGGTGCTATCCAGCATTAAATACTGGCCATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGa                                                         <  1:10543333/147‑1 (MQ=255)
                                                                                               cATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACCGTGCTATCCAGCATTAAAGACAGGCAATACTCTCACGCAAAAACCGCAGCACCGTCTCGTCACCGCAGCGAt                                                        >  2:8194595/1‑150 (MQ=38)
                                                                                                   ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                                    <  1:4713453/150‑1 (MQ=255)
                                                                                                   ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                                    <  2:2967569/150‑1 (MQ=255)
                                                                                                          ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                                             >  2:433719/1‑150 (MQ=255)
                                                                                                           tCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCAGCAGCACCGTCTCTTCAACGCTGCGATTCGGGTTAAttc                                            >  2:9787435/1‑148 (MQ=255)
                                                                                                              ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                                         >  1:5374444/1‑150 (MQ=255)
                                                                                                               cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                                        >  2:793693/1‑150 (MQ=255)
                                                                                                               cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                                        >  2:5159514/1‑150 (MQ=255)
                                                                                                               cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                                        <  1:11966574/150‑1 (MQ=255)
                                                                                                                tCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCGATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                                       <  2:6278513/150‑1 (MQ=255)
                                                                                                                    ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTCCGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCACTACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                                   >  2:3978337/1‑150 (MQ=255)
                                                                                                                    ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                                   >  1:2734928/1‑150 (MQ=255)
                                                                                                                    ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                                   >  1:6058231/1‑150 (MQ=255)
                                                                                                                    ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                                   >  1:10021883/1‑150 (MQ=255)
                                                                                                                      tAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa                                 >  1:9622847/1‑150 (MQ=255)
                                                                                                                          ccTCTCCCCAAAGGGTCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt                             <  1:1140481/150‑1 (MQ=255)
                                                                                                                          ccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAGTCGGGTTACTGCGAATCGCACaaaggg                             >  1:8205253/1‑146 (MQ=255)
                                                                                                                           ctctCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCGTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc                            <  2:4457711/150‑1 (MQ=255)
                                                                                                                           ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc                            >  2:11598599/1‑150 (MQ=255)
                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGACTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt                           >  2:8533859/1‑150 (MQ=255)
                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt                           >  1:7716612/1‑150 (MQ=255)
                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt                           >  1:6259531/1‑150 (MQ=255)
                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt                           >  2:11921573/1‑150 (MQ=255)
                                                                                                                            tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGCGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt                           >  2:5578697/1‑150 (MQ=255)
                                                                                                                             ctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTg                          >  2:840274/1‑150 (MQ=255)
                                                                                                                             ctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATACGGGTTAATGCGAATCGCACAAAGTTCTg                          >  1:11672832/1‑150 (MQ=255)
                                                                                                                               ccccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCACTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAt                        >  2:6044715/1‑150 (MQ=255)
                                                                                                                                 ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg                      >  2:1122024/1‑150 (MQ=255)
                                                                                                                                 ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg                      >  1:8965688/1‑150 (MQ=255)
                                                                                                                                        ggCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTc               >  1:11080650/1‑150 (MQ=255)
                                                                                                                                         gCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCg              >  1:3318993/1‑150 (MQ=255)
                                                                                                                                         gCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCg              >  1:4931035/1‑150 (MQ=255)
                                                                                                                                           gAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGcc            >  2:7378948/1‑150 (MQ=255)
                                                                                                                                            aGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCt           <  1:2025203/150‑1 (MQ=255)
                                                                                                                                            aGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCt           <  1:5668655/150‑1 (MQ=255)
                                                                                                                                             ggggACCGTTCAGGACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTg          <  1:3570642/150‑1 (MQ=255)
                                                                                                                                              gggACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTgg         >  2:12224660/1‑150 (MQ=255)
                                                                                                                                              gggACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTgg         <  2:7764502/150‑1 (MQ=255)
                                                                                                                                                gACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCg       >  1:1622110/1‑150 (MQ=255)
                                                                                                                                                    gTTCAGTACGGACTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAAc   <  2:10058634/150‑1 (MQ=255)
                                                                                                                                                    gTTCAGGACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAAc   <  2:11539844/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACG  >  NC_000913/3504749‑3505047

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: