Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,332 A→T 6.0% N388I (AAT→ATT) 
pseudogene (932/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,3320AT 6.0% 264.8 / 5.2 83N388I (AAT→ATT) 
pseudogene (932/1137 nt)
ymfN
aaaE
e14 prophage; chimeric replication protein/phage terminase YmfN
phage terminase
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base T (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACG  >  NC_000913/1205184‑1205481
                                                                                                                                                    |                                                                                                                                                     
gACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAAt                                                                                                                                                      >  2:650348/1‑150 (MQ=255)
    ccTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCca                                                                                                                                                  >  1:134310/1‑150 (MQ=255)
             cGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCGCCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCt                                                                                                                                         <  2:280379/150‑1 (MQ=255)
                ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGATGGGATTGCCGTCATGatgaaag                                                                                                                                      <  1:1375215/150‑7 (MQ=255)
                ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                      <  2:1090818/150‑1 (MQ=255)
                ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                      >  1:1074873/1‑150 (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:1552180/150‑1 (MQ=255)
                                          gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                            <  1:179769/150‑1 (MQ=255)
                                          gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                            <  1:1146670/150‑1 (MQ=255)
                                          gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                            >  2:936740/1‑150 (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:925269/150‑1 (MQ=255)
                                                                               ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                       >  2:863121/1‑150 (MQ=255)
                                                                               ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                       >  1:1445788/1‑150 (MQ=255)
                                                                               ccAACATGTCCGATCCGATGAAAGACCTGGAAGCGGCGATTGAATCGGAGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACAATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                       >  2:821862/1‑150 (MQ=255)
                                                                                  aCATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGc                                                                    >  2:1481131/1‑150 (MQ=255)
                                                                                   cATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCa                                                                   <  1:572658/150‑1 (MQ=255)
                                                                                     tGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAgg                                                                 <  1:653669/150‑1 (MQ=255)
                                                                                     tGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAgg                                                                 >  2:261457/1‑150 (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                                                         <  1:650681/150‑1 (MQ=255)
                                                                                             ccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAAc                                                         <  2:911887/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:82020/150‑1 (MQ=255)
                                                                                                                                              aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc        <  1:261457/150‑1 (MQ=255)
                                                                                                                                              aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc        <  2:1563600/150‑1 (MQ=255)
                                                                                                                                              aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc        <  1:1355528/150‑1 (MQ=255)
                                                                                                                                              aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc        <  1:1184373/150‑1 (MQ=255)
                                                                                                                                                ggCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTg      <  2:134310/150‑1 (MQ=255)
                                                                                                                                                    aTCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACg  <  1:1379944/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACG  >  NC_000913/1205184‑1205481

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: