Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,907 A→C 7.9% intergenic (+14/+21) mog → / ← satP molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,9070AC 7.9% 211.2 / 9.2 86intergenic (+14/+21)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base A (37/42);  new base C (4/3);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

GTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTG  >  NC_000913/9762‑10052
                                                                                                                                                  |                                                                                                                                                 
gTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                               >  1:1183135/1‑150 (MQ=255)
 ttGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                                                                              >  2:429819/1‑150 (MQ=255)
     ggTCCACGGTATTTTTGCCAGCGTACCGTACTGCCTTCAGTTGCTGG‑CAGGGCCATACGTTGAACCGGCACCGGAAGTGGTGGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAGATCCCCCCGCTCGGGGGGa                                                                                                                                          <  1:860769/150‑1 (MQ=255)
     ggTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGCGGGGa                                                                                                                                          >  2:161458/1‑150 (MQ=255)
       tACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAACCGGCACCGGAAGGGGTGGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                                                                        <  1:46732/150‑1 (MQ=255)
       tACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                                                                        <  1:999161/150‑1 (MQ=255)
          acGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                                                                     >  2:506540/1‑150 (MQ=255)
           cGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                                                                    <  2:418689/150‑6 (MQ=255)
           cGGTATTTATGCCCGCGTACCGTACTGCAGTCAGTTGCTGG‑ACGGGCCATACGTTCAACCGGCACCGGAAGTGGTTGCAGCATGCAGACCGAAGAGCGCATGACGCGACGTTAGCGAATAAAAAAATCCCGCCGCTCGGGGGGAtttttt                                                                                                                                    <  1:424537/150‑6 (MQ=255)
            ggTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACTGCACCGGAAGTGGTTGCAGAATTCAGAACGAATAGTGCAAGTCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaa                                                                                                                                   >  1:1262294/1‑150 (MQ=255)
              tATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                                                                 <  2:505972/150‑1 (MQ=255)
              tATTTTTGCCAGCGTACCGTACTGCATGCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                                                                 <  1:231661/150‑1 (MQ=255)
                     cccaccttaccgtccTGCATT‑TCTTTCTGGTAAGGGCCAAACGTTGAAACGCCCCCGGAAGGGGGGCCCGCAATCCCCCCCTAGAGGGCAAGCCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGTATCTCAAAACAATTAGt                                                                                                                          <  2:103617/143‑1 (MQ=255)
                      ccAGCGAACCGTACTGCATTCAGTTGATTG‑AAGGGCCATACGTGGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCTTGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                                                                         <  2:820149/150‑1 (MQ=255)
                       ccgCGTACCTTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGCCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACATTTAGTgg                                                                                                                        <  1:171725/148‑1 (MQ=255)
                       cAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGCCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                                                                        <  1:25199/150‑1 (MQ=255)
                        aGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                                                                       >  2:471813/1‑150 (MQ=255)
                             accgtacTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAACTCCCCCCGAGCTGGGGCATCTCAAAACAATTAGTGGGATTCa                                                                                                                  >  1:839512/1‑150 (MQ=255)
                              ccgtacTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                                                                 <  2:1310046/150‑1 (MQ=255)
                               cgtacTGCATTCTGTTGCTGG‑AAGGGCCATACGTTTAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                                                                <  1:1277331/150‑1 (MQ=255)
                               cgtacTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                                                                <  2:406980/150‑1 (MQ=255)
                                  acTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                                                                             >  1:783539/1‑150 (MQ=255)
                                    tGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAATTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                                                                           <  2:1610441/150‑1 (MQ=255)
                                     gCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                                                                          >  1:913023/1‑150 (MQ=255)
                                     gCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                                                                          >  1:442925/1‑150 (MQ=255)
                                     gCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                                                                          >  2:499203/1‑150 (MQ=255)
                                      cATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                         <  1:772496/150‑1 (MQ=255)
                                      cATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                         <  2:1576055/150‑1 (MQ=255)
                                      cATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                         >  1:1635552/1‑150 (MQ=255)
                                       aTTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                                                                        >  1:1036300/1‑150 (MQ=255)
                                        ttCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                                                                       <  2:1366883/150‑1 (MQ=255)
                                         tCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGa                                                                                                      <  2:1232714/150‑1 (MQ=255)
                                                   g‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                                                                            <  2:913023/150‑1 (MQ=255)
                                                          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACt                                                                                      <  1:920912/150‑1 (MQ=255)
                                                               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                                                                 >  2:1485022/1‑150 (MQ=255)
                                                               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                                                                 <  1:506540/150‑1 (MQ=255)
                                                                  tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAATACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                                                              <  1:1068402/150‑1 (MQ=255)
                                                                  tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                                                              <  1:1572988/150‑1 (MQ=255)
                                                                   gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                                                             >  2:407157/1‑150 (MQ=255)
                                                                     aaCGGCAGCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATAAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                           <  1:324203/150‑1 (MQ=255)
                                                                     aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                           <  2:1632558/150‑1 (MQ=255)
                                                                          cACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                                                                      >  1:777103/1‑150 (MQ=255)
                                                                             cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                                                                   <  1:504772/150‑1 (MQ=255)
                                                                              ggAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAACCTGCTCGTTCAGTACTTCAccc                                                                  >  1:723223/1‑150 (MQ=255)
                                                                                 aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                                                               >  1:975313/1‑150 (MQ=255)
                                                                                 aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                                                               >  1:1001939/1‑150 (MQ=255)
                                                                                 aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                                                               <  1:437906/150‑1 (MQ=255)
                                                                                  gTGGGTGCAGCATTCATTCCGAAGAGTTCAAGACGGGTCGTTAGCTAATAAAAAATACCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCGCGTTCAGTACTTCACCCATCg                                                              <  2:257881/150‑1 (MQ=255)
                                                                                    ggTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                                                            <  2:629414/150‑1 (MQ=255)
                                                                                     gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa                                                           >  2:1434444/1‑150 (MQ=255)
                                                                                      ttGCAGCATTCAGACCGAAGAGTGCAAGCCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAGTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGGACTTCCCCCAACGCCag                                                          >  2:295839/1‑150 (MQ=255)
                                                                                        gcagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAATCTGCTCGTTCAGTACTTCACCCATCGCCagat                                                        <  1:711512/150‑1 (MQ=255)
                                                                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata                                                       >  1:711719/1‑150 (MQ=255)
                                                                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata                                                       <  2:1475758/150‑1 (MQ=255)
                                                                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata                                                       >  1:189427/1‑150 (MQ=255)
                                                                                              ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTg                                                  <  1:903504/150‑1 (MQ=255)
                                                                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCGCGCGGGGCTCGCCAACCAAGTAGGGGGAATCACCAATCGGCAGACCGCTGCGACCAAACTGATCGTTCATGACTTCACCCATCGCCAGAAAGATTgcgc                                               >  1:98217/1‑150 (MQ=255)
                                                                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCATATAGCTTgcgc                                               >  1:1561768/1‑150 (MQ=255)
                                                                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcgc                                               >  1:1069377/1‑150 (MQ=255)
                                                                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACGAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCt                                              <  1:429819/150‑1 (MQ=255)
                                                                                                   aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGGCCAAACTGCTCGATCAGTACTTCACCCATCGCCAGATACATTGCGCTg                                             >  2:1015079/1‑150 (MQ=255)
                                                                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGATCGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTgg                                            <  1:864827/150‑1 (MQ=255)
                                                                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCa                                          >  2:1510206/1‑150 (MQ=255)
                                                                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCa                                          <  2:1010038/150‑1 (MQ=255)
                                                                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGACCTTCACCCATCGCCAGATAGATTGCGCTGGcc                                          >  1:113210/1‑149 (MQ=255)
                                                                                                          agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCACAAAACAATTAGTGGGATTCACAAATCGGCAGAACGGTGCGACAAAACTGCTCGTTCAGTACTTCACCCATCGCCAGAAAGATTGCGCTGGCAcccc                                      >  2:274926/1‑148 (MQ=255)
                                                                                                           gTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAGTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAACGCCAGATAGATTGCGCTGGCACCGCa                                     >  1:641663/1‑150 (MQ=255)
                                                                                                            tGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAACCTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAg                                    >  1:139095/1‑150 (MQ=255)
                                                                                                             gCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGa                                   <  1:1592380/150‑1 (MQ=255)
                                                                                                              cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGAt                                  <  2:857752/150‑1 (MQ=255)
                                                                                                                aGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCGAGATAGATTGCGCTGGCACCGCAGATcc                                >  2:1213067/1‑149 (MQ=255)
                                                                                                                aGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTCCGACCAACCTGCTCGTTCAGTACTTCACCCATCTCAAGATAGATTGCGCTGGCACCGCAGATCa                                >  2:1440069/1‑150 (MQ=255)
                                                                                                                         gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATcc                       <  1:499203/150‑1 (MQ=255)
                                                                                                                              cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCg                  <  2:702638/150‑1 (MQ=255)
                                                                                                                                   aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCCCCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCaaa             >  1:607709/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaTCCCACCGAGCGGGGGGATCTCATAACATTTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCACACTGCTCGTTCAGTACTTCACCCATAGCCAGATAGATTGCGCTGGCACCGCAGATCAGCACGATCACACCGGCaaa             >  1:81257/1‑150 (MQ=255)
                                                                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGTCACCGCAGATCAGCCCAATCCAGCCGGCAAAg            <  2:442925/150‑1 (MQ=255)
                                                                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGATGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGCTCAGCCCAATCCAGCCGGCAAAg            >  1:636711/1‑150 (MQ=255)
                                                                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCCGAGCAGCCCAATCCAGCCGGCAAAGTg          >  2:1541422/4‑150 (MQ=255)
                                                                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTg          >  2:1678589/4‑150 (MQ=255)
                                                                                                                                      aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTg          <  2:1036300/150‑1 (MQ=255)
                                                                                                                                      aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTg          <  2:125561/150‑1 (MQ=255)
                                                                                                                                            cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGCCCAAACTGCTCGTTCAGGACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAAGCCAGCCGGCAAAGTGgatgat    >  2:601545/1‑150 (MQ=255)
                                                                                                                                            cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGgatgat    <  2:1274356/150‑1 (MQ=255)
                                                                                                                                             ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGAtt   >  1:356418/1‑150 (MQ=255)
                                                                                                                                              cccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTg  <  1:1463448/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGG‑AAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTG  >  NC_000913/9762‑10052

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