Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,015,583 A→C 6.3% M142L (ATG→CTG) ‡ mocA → CTP:molybdopterin cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,015,5830AC 6.3% 278.7 / 7.2 79M142L (ATG→CTG) ‡mocACTP:molybdopterin cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base C (3/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTT  >  NC_000913/3015322‑3015817
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
aaTGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCa                                                                                                                                                                                                                               <  2:436260/275‑1 (MQ=255)
aaTGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAAATTATGCGCAAGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCa                                                                                                                                                                                                                               <  2:74219/275‑1 (MQ=255)
                        aaCCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCg                                                                                                                                                                                                       >  2:126642/1‑275 (MQ=255)
                             gagCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAgg                                                                                                                                                                                                  <  2:68174/275‑1 (MQ=255)
                                             attCACAACCCAGATTATGCGCAGAGTTTACTAACATCCGTCAAAACCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGCAAAATCTGGTCGTTACGAAATGATGGCGCAATGCTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAGCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGga                                                                                                                                                                                    <  2:432983/273‑1 (MQ=255)
                                                   aaCCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGCGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCACCTGCCATCCGACTTTAGTGTCAAAACCATGCCTGCTGCAGGCAATCCAACGACCCAATCACACCAATATGCGTCAGCCAACCCTATTGGGAGATCattc                                                                                                                                                                            >  1:148846/1‑274 (MQ=255)
                                                        aGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCg                                                                                                                                                                       >  1:39527/1‑275 (MQ=255)
                                                                    gggTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGa                                                                                                                                                                                                                                            >  2:119126/1‑194 (MQ=255)
                                                                               cATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCCCCAATATGCGTCa                                                                                                                                                                                                    >  4:271401/1‑223 (MQ=255)
                                                                                   cGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCa                                                                                                                                                                                                <  2:271405/223‑1 (MQ=255)
                                                                                     tCAAGGCCGCAGTACCCGCAGTGCAAACAGAACAGTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGGTATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCaa                                                                                                                                                                                                                              <  2:305965/191‑1 (MQ=255)
                                                                                       aaGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGa                                                                                                                                        <  2:131767/275‑1 (MQ=255)
                                                                                                    ccGNNGTGCAAACAGAACATTGTTTGCTCACCCACGGTGACATGCCAACCCTCNCCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCa                                                                                                                                                                                                <  4:512617/206‑1 (MQ=255)
                                                                                                         gTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGgac                                                                                                                        <  2:253965/273‑1 (MQ=255)
                                                                                                         gTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGgac                                                                                                                        <  2:480996/273‑1 (MQ=255)
                                                                                                                  gAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAAt                                                                                                                                                                                                                                  >  2:329208/1‑158 (MQ=255)
                                                                                                                   aaCATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATg                                                                                                                                                                                                                                          <  8:286557/149‑1 (MQ=255)
                                                                                                                   aaCATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCAATTTTAGTGTCAAAACCATGCCTGATg                                                                                                                                                                                                                                          >  8:201116/1‑149 (MQ=255)
                                                                                                                      aTTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCaa                                                                                                                                                                                                                              <  4:329204/158‑1 (MQ=255)
                                                                                                                       ttGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAgg                                                                                                                                                                                                                                      >  6:260262/1‑149 (MQ=255)
                                                                                                                          ttttCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCaa                                                                                                                                                                                                                                   <  7:11577/149‑1 (MQ=255)
                                                                                                                               tcACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTTGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCatggttt                                                                                                                                                                                                                              >  5:307085/1‑143 (MQ=255)
                                                                                                                               tcACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCatggttt                                                                                                                                                                                                                              <  6:307085/149‑7 (MQ=255)
                                                                                                                                   ccACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGAc                                                                                                                                                                                                                          <  6:68384/149‑1 (MQ=255)
                                                                                                                                   ccACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGAc                                                                                                                                                                                                                          >  8:141110/1‑149 (MQ=255)
                                                                                                                                             aTGCCGACCCTCACCATCGATATTTTTAGAAAAATCTAGTCGTTACAAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCAc                                                                                                                                                                                                                <  6:88071/149‑1 (MQ=255)
                                                                                                                                             aTGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCAc                                                                                                                                                                                                                <  5:149593/149‑1 (MQ=255)
                                                                                                                                                ccAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCaa                                                                                                                                                                                                             >  8:84324/1‑149 (MQ=255)
                                                                                                                                                   aCCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATGATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAAc                                                                            >  2:106124/1‑275 (MQ=255)
                                                                                                                                                   aCCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAatat                                                                                                                                                                                                          <  7:416088/149‑1 (MQ=255)
                                                                                                                                                   aCCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGga                                                                                                                                                                                    <  2:216622/171‑1 (MQ=255)
                                                                                                                                                     ccTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGc                                                                                                                                                                                                        >  6:76503/1‑149 (MQ=255)
                                                                                                                                                     ccTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGc                                                                                                                                                                                                        >  6:50502/1‑149 (MQ=255)
                                                                                                                                                     ccTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGc                                                                                                                                                                                                        <  7:356066/149‑1 (MQ=255)
                                                                                                                                                       tCACCATCGATATTTTTAGAAAGATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACGAATATGCGt                                                                                                                                                                                                      <  5:408882/149‑1 (MQ=255)
                                                                                                                                                         aCCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCaa                                                                                                                                                                                                                                   >  5:17720/1‑118 (MQ=255)
                                                                                                                                                         aCCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCaa                                                                                                                                                                                                                                   <  6:17720/118‑1 (MQ=255)
                                                                                                                                                              cGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATACTTTTAGATATTGATACTCCGGACGACTTTATTACAGCGa                                                                                                        >  2:381557/1‑236 (MQ=255)
                                                                                                                                                                         gAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCat                                                                                                                                                                                                                                   >  8:291641/1‑101 (MQ=255)
                                                                                                                                                                         gAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCat                                                                                                                                                                                                                                   <  7:291641/102‑2 (MQ=255)
                                                                                                                                                                         gAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCAATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGga                                                                                                                                                                                    <  7:201116/149‑1 (MQ=255)
                                                                                                                                                                          aaaaaTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGgag                                                                                                                                                                                   >  6:220479/1‑149 (MQ=255)
                                                                                                                                                                          aaaaaTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGgag                                                                                                                                                                                   >  6:300702/1‑149 (MQ=255)
                                                                                                                                                                            aaaTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGAt                                                                                                                                                                                 >  8:243753/1‑149 (MQ=255)
                                                                                                                                                                              aTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCa                                                                                                                                                                               >  6:69305/1‑149 (MQ=255)
                                                                                                                                                                                      gTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGACTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCg                                                                                                                                                                       >  5:271706/1‑149 (MQ=255)
                                                                                                                                                                                       ttACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGt                                                                                                                                                                      >  6:421436/1‑149 (MQ=255)
                                                                                                                                                                                        tACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCaa                                                                                                                                                                                                                              >  2:353280/1‑92 (MQ=255)
                                                                                                                                                                                           gAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGaa                                                                                                                                                                  >  7:440144/1‑149 (MQ=255)
                                                                                                                                                                                            aaaTGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGAc                                                                                                                                                                                                                          <  4:353275/92‑1 (MQ=255)
                                                                                                                                                                                                    gcgcAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAAt                                                                                                                                                         <  7:24116/149‑1 (MQ=255)
                                                                                                                                                                                                       cAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCa                                                                                                                                                      <  8:280280/149‑1 (MQ=255)
                                                                                                                                                                                                        aaaacgaaaacgccAAAATTGCATCCCCTNNCATCCAATTTTATTGTCATAACCATGCCTCNTGCAGGCAATCCAACGACCCAATGTCACCAATATACGTCAGGCACTACTTATGGGAGATCATTATTCCGTCGAANTAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGAACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                       <  3:75243/263‑1 (MQ=255)
                                                                                                                                                                                                               ccACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaataat                                                                                                                                              >  5:212938/1‑149 (MQ=255)
                                                                                                                                                                                                               ccACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATAATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATa                 <  4:106123/274‑1 (MQ=255)
                                                                                                                                                                                                                cACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAAtt                                                                                                                                             <  5:69305/149‑1 (MQ=255)
                                                                                                                                                                                                                  cTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAAtttt                                                                                                                                           <  5:300702/149‑1 (MQ=255)
                                                                                                                                                                                                                   tCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTa                                                                                                                                          >  6:221234/1‑149 (MQ=255)
                                                                                                                                                                                                                     cATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGa                                                                                                                                        <  5:260262/149‑1 (MQ=255)
                                                                                                                                                                                                                     cATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGa                                                                                                                                        >  7:157281/1‑149 (MQ=255)
                                                                                                                                                                                                                          tGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGGCCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTg                                                                                                                                   <  6:446521/149‑1 (MQ=255)
                                                                                                                                                                                                                           ggCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGa                                                                                                                                  <  7:243753/149‑1 (MQ=255)
                                                                                                                                                                                                                             cATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGAttt  >  4:217444/1‑275 (MQ=255)
                                                                                                                                                                                                                               tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACt                                                                                                                              >  7:177830/1‑149 (MQ=255)
                                                                                                                                                                                                                               tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatt                                                                                                               >  4:551082/1‑164 (MQ=255)
                                                                                                                                                                                                                               tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa                                                                                                       >  2:491974/1‑172 (MQ=255)
                                                                                                                                                                                                                               tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTAt                           >  1:75243/1‑248 (MQ=255)
                                                                                                                                                                                                                                   cGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAg                                                                                                           <  2:551087/164‑1 (MQ=255)
                                                                                                                                                                                                                                         tCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTAt                                                                                                                                                                                        >  2:534840/1‑81 (MQ=255)
                                                                                                                                                                                                                                             aTTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGga                                                                                                                                                                                    <  4:534835/81‑1 (MQ=255)
                                                                                                                                                                                                                                             aTTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttat                                                                                                                <  5:421436/149‑1 (MQ=255)
                                                                                                                                                                                                                                               tttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatta                                                                                                              <  5:50502/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa                                                                                                       <  8:177830/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAgg                                                                                                     >  7:2108/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTAt                                                                                                                                                                                        >  5:358756/1‑65 (MQ=255)
                                                                                                                                                                                                                                                         aaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTAt                                                                                                                                                                                        <  6:358756/65‑1 (MQ=255)
                                                                                                                                                                                                                                                           aCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaaa                                                                                                  <  8:2108/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGta                                                                                              <  8:157281/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                gCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCa                                                                                                                                                                                                    >  4:518223/1‑46 (MQ=255)
                                                                                                                                                                                                                                                                 ccTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGtata                                                                                            <  6:374529/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCa                                                                                                                                                                                                <  2:518228/46‑1 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
AATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTT  >  NC_000913/3015322‑3015817

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