Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,337 A→G 5.7% I390V (ATC→GTC) 
pseudogene (937/1137 nt)
ymfN →
aaaE →
e14 prophage; chimeric replication protein/phage terminase YmfN
phage terminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,205,3370AG 5.7% 265.2 / 2.5 86I390V (ATC→GTC) 
pseudogene (937/1137 nt)
ymfN
aaaE
e14 prophage; chimeric replication protein/phage terminase YmfN
phage terminase
Reads supporting (aligned to +/- strand):  ref base A (38/43);  new base G (3/2);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAA  >  NC_000913/1205188‑1205485
                                                                                                                                                     |                                                                                                                                                    
ccTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCca                                                                                                                                                      >  1:134310/1‑150 (MQ=255)
         cGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCGCCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCt                                                                                                                                             <  2:280379/150‑1 (MQ=255)
            ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGATGGGATTGCCGTCATGatgaaag                                                                                                                                          <  1:1375215/150‑7 (MQ=255)
            ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                          >  1:1074873/1‑150 (MQ=255)
            ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                          <  2:1090818/150‑1 (MQ=255)
            ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                          <  1:1439135/150‑1 (MQ=255)
             ggtgCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgtt                                                                                                                                         <  2:363606/147‑2 (MQ=255)
             cGGGCTGTCACATGACCTTGCTGATGACGACCGGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgtg                                                                                                                                         <  1:1121455/150‑1 (MQ=255)
             cGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgtg                                                                                                                                         <  1:650348/150‑1 (MQ=255)
                   gTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCgg                                                                                                                                   <  2:121940/150‑1 (MQ=255)
                    tCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGc                                                                                                                                  <  2:525518/150‑1 (MQ=255)
                     cacaTGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCAATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCa                                                                                                                                 >  2:1131504/1‑150 (MQ=255)
                      acaTGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCaa                                                                                                                                >  2:1184373/1‑150 (MQ=255)
                        aTGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACg                                                                                                                              >  1:1563600/1‑150 (MQ=255)
                         tGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGt                                                                                                                             <  1:989884/150‑1 (MQ=255)
                            ccTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGt                                                                                                                          >  2:953245/1‑150 (MQ=255)
                            ccTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGt                                                                                                                          >  2:82020/1‑150 (MQ=255)
                               tGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCgg                                                                                                                       >  2:552272/1‑150 (MQ=255)
                                 cTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCa                                                                                                                     <  2:1074873/150‑1 (MQ=255)
                                    atgaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaaa                                                                                                                  >  1:911887/1‑150 (MQ=255)
                                      gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                                >  2:936740/1‑150 (MQ=255)
                                      gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                                <  1:1146670/150‑1 (MQ=255)
                                      gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                                <  1:179769/150‑1 (MQ=255)
                                      gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                                <  2:1552180/150‑1 (MQ=255)
                                           ccTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTcc                                                                                                           <  1:1009933/150‑1 (MQ=255)
                                            cTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCg                                                                                                          >  2:188785/1‑150 (MQ=255)
                                             tGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCgg                                                                                                         >  1:662721/1‑150 (MQ=255)
                                               aaCCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGt                                                                                                       >  2:653669/1‑150 (MQ=255)
                                                 ccccGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTaa                                                                                                     >  1:11326/1‑150 (MQ=255)
                                                   ccGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGAGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACg                                                                                                   >  2:216310/1‑150 (MQ=255)
                                                    cGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACga                                                                                                  <  2:371240/150‑1 (MQ=255)
                                                       cACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACgatga                                                                                               >  1:617183/1‑150 (MQ=255)
                                                         ccATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGatgatg                                                                                             >  2:572658/1‑150 (MQ=255)
                                                         ccATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGatgatg                                                                                             >  2:1379944/1‑150 (MQ=255)
                                                                tCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGaa                                                                                      >  1:1378686/1‑150 (MQ=255)
                                                                 cAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAg                                                                                     >  2:749051/1‑150 (MQ=255)
                                                                 cAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAg                                                                                     >  2:659546/1‑150 (MQ=255)
                                                                       tACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTc                                                                               <  2:1076264/150‑1 (MQ=255)
                                                                        acacCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCa                                                                              >  2:1355528/1‑150 (MQ=255)
                                                                         cgccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAACACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCaa                                                                             >  2:763939/3‑150 (MQ=255)
                                                                           ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                           >  1:1445788/1‑150 (MQ=255)
                                                                           ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                           >  2:863121/1‑150 (MQ=255)
                                                                           ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                           <  1:925269/150‑1 (MQ=255)
                                                                           ccAACATGTCCGATCCGATGAAAGACCTGGAAGCGGCGATTGAATCGGAGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACAATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                           >  2:821862/1‑150 (MQ=255)
                                                                              aCATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGc                                                                        >  2:1481131/1‑150 (MQ=255)
                                                                               cATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCa                                                                       <  1:572658/150‑1 (MQ=255)
                                                                                 tGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAgg                                                                     >  2:261457/1‑150 (MQ=255)
                                                                                 tGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAgg                                                                     <  1:653669/150‑1 (MQ=255)
                                                                                    ccgatccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCgg                                                                  <  1:936740/150‑1 (MQ=255)
                                                                                     cgatccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGa                                                                 <  2:322765/150‑1 (MQ=255)
                                                                                     cgatccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGa                                                                 >  1:275878/1‑150 (MQ=255)
                                                                                        tccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGaaaa                                                              <  2:143988/150‑1 (MQ=255)
                                                                                         ccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAAc                                                             <  2:911887/150‑1 (MQ=255)
                                                                                         ccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAAc                                                             <  1:650681/150‑1 (MQ=255)
                                                                                             tGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACaaaa                                                         <  1:1131504/150‑1 (MQ=255)
                                                                                              gAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAAt                                                        >  1:73245/1‑150 (MQ=255)
                                                                                                aaGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCg                                                      <  1:228232/150‑1 (MQ=255)
                                                                                                aaGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCg                                                      >  2:202618/1‑150 (MQ=255)
                                                                                                  gagCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCAGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGAt                                                    <  2:270229/150‑1 (MQ=255)
                                                                                                   cgCTGGAAGCGGCGATTGAAGCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGCTGTATCGGCAACGTGGTCGGCACAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATg                                                   <  2:517616/149‑1 (MQ=255)
                                                                                                       ggAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGc                                               <  2:834433/150‑1 (MQ=255)
                                                                                                        gAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCa                                              >  1:371052/1‑150 (MQ=255)
                                                                                                         aaGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAg                                             <  2:1336792/150‑1 (MQ=255)
                                                                                                             ggcgATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTgc                                         >  1:249751/1‑150 (MQ=255)
                                                                                                              gcgATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTgcg                                        >  2:169360/1‑150 (MQ=255)
                                                                                                                     aaTCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGaaa                                                               <  1:1064662/120‑1 (MQ=255)
                                                                                                                     aaTCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTAAAGGAGCAGGCGGaaa                                                               >  2:1064662/1‑120 (MQ=255)
                                                                                                                        cGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATgg                              <  2:89617/150‑1 (MQ=255)
                                                                                                                         ggggCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGc                             >  1:977962/1‑150 (MQ=255)
                                                                                                                         ggggCGCTTTCATCATGACGGCAATCCCATCATGACCTGGGGTATCGGCAACGTGGGCGGCAAAACCATTCCGGGTAACGATGAGGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGc                             <  2:1053463/150‑1 (MQ=255)
                                                                                                                           ggCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCgg                           <  1:501600/150‑1 (MQ=255)
                                                                                                                              gcTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTg                        >  1:445927/1‑150 (MQ=255)
                                                                                                                                ccaggtcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGACCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGc                      >  2:1121955/6‑150 (MQ=255)
                                                                                                                                  aggtcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCag                    >  2:965231/4‑150 (MQ=255)
                                                                                                                                    gtcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCagag                  >  2:308260/2‑150 (MQ=255)
                                                                                                                                     tcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGc                 <  1:308260/150‑1 (MQ=255)
                                                                                                                                          aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc            <  1:1184373/150‑1 (MQ=255)
                                                                                                                                          aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc            <  1:1355528/150‑1 (MQ=255)
                                                                                                                                          aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc            <  1:261457/150‑1 (MQ=255)
                                                                                                                                          aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc            <  1:82020/150‑1 (MQ=255)
                                                                                                                                          aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc            <  2:1563600/150‑1 (MQ=255)
                                                                                                                                            ggCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTg          <  2:134310/150‑1 (MQ=255)
                                                                                                                                                aTCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACg      <  1:1379944/150‑1 (MQ=255)
                                                                                                                                                  cccATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGa                                                                          >  2:1093660/1‑80 (MQ=255)
                                                                                                                                                  cccATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGa                                                                          <  1:1093660/80‑1 (MQ=255)
                                                                                                                                                  cccATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACgag    <  1:749051/150‑1 (MQ=255)
                                                                                                                                                   ccATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACgaga   >  2:1584535/1‑150 (MQ=255)
                                                                                                                                                    catcatGACCTGGTGTATCGGAAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGTGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGAAGAGCCATGCTGTACGagaa  <  2:322256/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAA  >  NC_000913/1205188‑1205485

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