Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 515,859 C→G 100% intergenic (‑86/‑60) qmcA ← / → fetA PHB domain‑containing protein/putative iron ABC exporter ATP‑binding subunit FetA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913515,8590CG98.4% 214.3 / ‑6.2 62intergenic (‑86/‑60)qmcA/fetAPHB domain‑containing protein/putative iron ABC exporter ATP‑binding subunit FetA
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base G (32/28);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGT  >  NC_000913/515603‑516122
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                        
gCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGCATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCACCTTATTCCTGCg                                                                                                                                                                                                                                                           >  1:473769/1‑272 (MQ=255)
 cccccAGTGCGCCCCATAAACGGCCCCACCGGACTGGGCCCCGGCTGGAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGa                                                                                                                                                                                                                                                          <  2:473768/269‑1 (MQ=255)
                 tAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGAATACTGAAAGAGAGTTGCCCCTTATTACTGCGAAGGATAAACTGttttt                                                                                                                                                                                                                                          >  1:1004634/1‑272 (MQ=255)
                            accaGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTAATTCCCGCGAAGGAAAAACTGTGTTTAGTAAAAATCa                                                                                                                                                                                                                               >  2:89491/1‑272 (MQ=255)
                                        cccGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAAc                                                                                                                                                                                                                   >  1:131375/1‑272 (MQ=255)
                                                aacccttttGGTGATCGAGACCAAACCGTTCTCCTGTCCACGGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGAAAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGGTAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                                                                                                                                                                                                            <  2:1004633/267‑1 (MQ=255)
                                                                      aaaGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTcc                                                                                                                                                                                                                                                               >  1:1014596/1‑199 (MQ=255)
                                                                      aaaGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACACGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATACAGATAAGCATCGGAAAAACCTACTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCACCTTATTACTGCGAAGGAAAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAAAAGCGATGCAGCAAAATAgaaatatgcatc                                                                                                                                                                                      >  1:174349/1‑261 (MQ=255)
                                                                       aaGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCt                                                                                                                                                                                                                                                              <  2:1014595/199‑1 (MQ=255)
                                                                                                 cTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTCCCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCgg                                                                                                                                                           >  1:1056866/1‑272 (MQ=255)
                                                                                                 cTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCg                                                                                                                                                            >  1:180762/1‑271 (MQ=255)
                                                                                                       tACGATTTTGACTCCCGCTCCGACAATGACCAGCTCGAATAAAAGGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATg                                                                                                                                                     <  1:913892/272‑1 (MQ=255)
                                                                                                                           cGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCt                                                                                                                                                                                                                                                              >  1:661645/1‑147 (MQ=255)
                                                                                                                            gACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTg                                                                                                                                                                                                                                                             <  2:661644/147‑1 (MQ=255)
                                                                                                                                           gACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGc                                                                                                                                                                                                                >  1:985781/1‑177 (MQ=255)
                                                                                                                                           gACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAg                                                                                                                                                                                                 >  1:1064721/1‑192 (MQ=255)
                                                                                                                                           gACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTGGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGATGCGTg                                                                                                                 >  1:513279/1‑272 (MQ=255)
                                                                                                                                            aCAAAAATGCGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCg                                                                                                                                                                                                               <  2:985780/177‑1 (MQ=255)
                                                                                                                                            aCAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGt                                                                                                                                                                                                <  2:1064720/192‑1 (MQ=255)
                                                                                                                                             cAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTTCTGTTGTACCGTCCATAAGCAGCAAAATTTCTTCTTGATTAAACAAATTATACCTGATTACTGAAAGATAGTTGCCTCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGaaaa                                                                                                                                                                                                    <  1:903013/187‑1 (MQ=255)
                                                                                                                                                   ggagaATGGTGATACAGTATAGCTTCGTAACATCCTCCTGTTGTCCCGTCCATAATCAGCAAAATTGTTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGaa                                                                                                         <  2:174349/271‑1 (MQ=255)
                                                                                                                                                         tCGGGATAAAGATAAGCATCGGAAAACCCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTaa                                                                                                    <  1:89491/271‑1 (MQ=255)
                                                                                                                                                                  ataTAAGCATCGGAAAAACCTCCAGTTTTACCGTCCCTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTAcc                                                                                          <  2:131375/270‑1 (MQ=255)
                                                                                                                                                                    aTAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGAt                                                                                                                                                                                                             >  2:124737/1‑155 (MQ=255)
                                                                                                                                                                     tAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                                                                                                                                                                                                            <  1:124737/155‑1 (MQ=255)
                                                                                                                                                                          aTCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATc                                                                                                                                                                 >  1:738155/1‑193 (MQ=255)
                                                                                                                                                                           tCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCt                                                                                                                                                                <  2:738154/193‑1 (MQ=255)
                                                                                                                                                                               aaaaaCCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGt                                                                              >  1:968548/1‑271 (MQ=255)
                                                                                                                                                                                           ttgtACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCAAATTTAAGTTaa                                                                                               >  2:111727/1‑242 (MQ=255)
                                                                                                                                                                                            tgtACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAAt                                                                                              <  1:111727/242‑1 (MQ=255)
                                                                                                                                                                                                     cATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGa                                                                                                                                                                                                                             >  1:70378/1‑106 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTTCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAACACGTAGGATATCTGGCGGGTGATGCGAAGATACTTAATAACATCACTTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACAGGACCTTCTGGTTGTGGCAAAATAACGCTACTaaaa                                                     >  1:1084244/1‑272 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaa                                                                                                                                                                                                                            <  2:70378/106‑1 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAgg                                                                                                                                                                                                        >  1:97544/1‑126 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGc                                                                                                                                                                                   >  1:877337/1‑147 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGAt                                                                                                                                                      >  2:877235/1‑176 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCg                                                                                                                                                   >  2:608530/1‑179 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAAttttt                                                                                                                          >  2:167392/1‑204 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTaaaa                                                     >  2:965658/1‑272 (MQ=255)
                                                                                                                                                                                                      aTAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCAGGTTGTGGGCAAAAgaaagggtcaaaa                                                      >  2:1067195/1‑260 (MQ=255)
                                                                                                                                                                                                       tAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGa                                                                                                                                                                                                       <  2:97544/126‑1 (MQ=255)
                                                                                                                                                                                                       tAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCt                                                                                                                                                                                  <  2:877336/147‑1 (MQ=255)
                                                                                                                                                                                                       tAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATg                                                                                                                                                     <  1:877236/176‑1 (MQ=255)
                                                                                                                                                                                                       tAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGa                                                                                                                                                  <  1:608531/179‑1 (MQ=255)
                                                                                                                                                                                                       tAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTc                                                                                                                         <  1:167392/204‑1 (MQ=255)
                                                                                                                                                                                                        aaTCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTAc                                                                                           >  2:775157/1‑233 (MQ=255)
                                                                                                                                                                                                         aTCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTAcc                                                                                          <  1:775158/233‑1 (MQ=255)
                                                                                                                                                                                                           cAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGt                                                >  1:30184/1‑272 (MQ=255)
                                                                                                                                                                                                                aaaTTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCtt                                           >  2:102463/1‑272 (MQ=255)
                                                                                                                                                                                                                aaaTTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCtt                                           >  1:297247/1‑272 (MQ=255)
                                                                                                                                                                                                                        ttttttttAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATc                                   <  2:968547/266‑1 (MQ=255)
                                                                                                                                                                                                                            gATTAAACAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGcc                               >  1:852006/1‑272 (MQ=255)
                                                                                                                                                                                                                                aaaCAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAc                           <  1:102463/272‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaaaTTTTACCAGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCa                         <  2:297246/270‑1 (MQ=255)
                                                                                                                                                                                                                                  aCAAATTATACCTGATTACTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCa                         <  1:1067196/272‑1 (MQ=255)
                                                                                                                                                                                                                                          tacttttttatgaaaaGGGAGTTCCCCCTTTTTCCTGCGAAGGTTAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGAGTCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATACTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACg                 <  2:1084243/259‑1 (MQ=255)
                                                                                                                                                                                                                                                  taagagttaggagTTGCTCCTTTTTCCGTAGAAGTATAACTTGTTATTAGTAAAAATAAAAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTGCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACGTCAATTTTTCGCTGCGTGCTGGCTAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGtt         <  1:1045683/262‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCTGAAAGAGATTTGCCCCCTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGAt                                                                                                                                                      >  1:302481/1‑131 (MQ=255)
                                                                                                                                                                                                                                                   aCTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCa                                         >  1:392146/1‑239 (MQ=255)
                                                                                                                                                                                                                                                    cTGAAAGAGATTTGCCCCCTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATg                                                                                                                                                     <  2:302480/131‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTGAAAGAGAGTTGCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCAt                                        <  2:392145/239‑1 (MQ=255)
                                                                                                                                                                                                                                                         gtaggtttCCCCTTTTTTCGGGAGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGt  <  2:1056865/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                        
GCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTTCAGCTACAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGT  >  NC_000913/515603‑516122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: