Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,517,171 G→A L46F (CTT→TTT)  damX ← cell division protein DamX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,517,1710GA98.5% 209.3 / NA 69L46F (CTT→TTT) damXcell division protein DamX
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/41);  minor base C (1/0);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAATCG  >  NC_000913/3516925‑3517416
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
aGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATGAAATTCa                                                                                                                                                                                                                                                   <  1:293164/251‑1 (MQ=255)
           aGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCaaa                                                                                                                                                                                                                                        <  1:597595/251‑1 (MQ=255)
            gcgcAGAACCGATACCGATGAGCAACAGCAGTAGAACCATAATGACGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTGTTGCGCTGACGCGGACGACTCATAACGCGCTCTTCAGCCACGGTTTCATCTTCGGATTCGACAATTGCTTCTTCGAGTGCCGTTTCCTCAGTGATCTCAATTTGCGCTCGAGTCGGACGGCGATCGTCAGTGGCAGCAAACTCAATACCAGCaaaa                                                                                                                                                                                                                                       <  2:627020/251‑1 (MQ=255)
               cAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTg                                                                                                                                                                                                                                    <  2:301424/251‑1 (MQ=255)
                aGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTAAGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTg                                                                                                                                                                                                                                    >  1:372626/1‑250 (MQ=255)
                  aaCCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATc                                                                                                                                                                                                                                 >  1:403735/1‑251 (MQ=255)
                    ccgataccgatGATCAACAGCAGTAGAACCAGAATGCCTACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAAGATCATCAAAATTGATCTg                                                                                                                                                                                                                               <  1:239417/251‑1 (MQ=255)
                      gatcccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACTCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGGCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCg                                                                                                                                                                                                                             <  2:398115/247‑1 (MQ=255)
                          ccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTc                                                                                                                                                                                                                         <  1:64044/251‑1 (MQ=255)
                           cgatGATCAACAGCAGTAGACCCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGATTACTGGCTGCTTTTTTGCGCTTCCGCGGACGACGCTCTACGCGCTCTTCCTACACGGTTTCATCTTCGGATTCCTCCATTTCTTCTTCTATTTCCGGTTCCTCATTGCTCTCTTTTTGCGCACGAGTCGGCCGGCGCTCGTCCGTGTCCTCCACTTCACTATCATCAAAATTGATCTGCGGTTCg                                                                                                                                                                                                                        >  2:425659/1‑251 (MQ=255)
                            gatgatCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGc                                                                                                                                                                                                                       >  1:4423/1‑251 (MQ=255)
                                           tAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCt                                                                                                                                                                                                        >  1:195134/1‑251 (MQ=255)
                                            aGAACCAGAAGGCCGACGACCATCATAATATACTGACGAGAAGCGGGTATACTGGCTGCTTTTTTGCGCTTACGCGGGCGACGCTCTACGCGCTCTTCCTCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTACGCTCTTTTTGCGCACGAGTCGGACGGCGCTCGTCAGTGTCATCACATTCAATGTCATCAAAATTGATCTGCGGTTCGCCCCTTTGAGTACgttc                                                                                                                                                                                                       >  2:474527/1‑248 (MQ=255)
                                                                catcatcatATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGAc                                                                                                                                                                                   >  1:580699/1‑251 (MQ=255)
                                                                   catcatATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCACTTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGCTCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACg                                                                                                                                                                                  >  2:26306/1‑249 (MQ=255)
                                                                   catcatATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACga                                                                                                                                                                                 <  2:89965/250‑1 (MQ=255)
                                                                    atcatATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACg                                                                                                                                                                                  >  1:118747/1‑248 (MQ=255)
                                                                          taCTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGAAGAGAACGAcc                                                                                                                                                                         >  2:92403/1‑251 (MQ=255)
                                                                              gACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAg                                                                                                                                                                       <  1:600837/249‑1 (MQ=255)
                                                                                 gagaAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTacga                                                                                                                                                                  <  1:593153/251‑1 (MQ=255)
                                                                                      gCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCGTCATCCACGGCGTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGATCAGTACGCTCAGAACGTGCAGAAGATTGACGAGAACGACCAGTACGACGATc                                                                                                                                                             <  2:301522/251‑1 (MQ=255)
                                                                                                    tgctTTTTTTCGATTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGCTCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGtt                                                                                                                                               >  2:630340/1‑251 (MQ=255)
                                                                                                         ttttGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAg                                                                                                                                           <  2:452834/250‑1 (MQ=255)
                                                                                                         ttttGCGCTTAAGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGCCCAGTACGACGATCGCTGGTATCGGGTTTCAGc                                                                                                                                          >  2:460116/1‑251 (MQ=255)
                                                                                                             gcgcTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCg                                                                                                                                       <  1:619558/250‑1 (MQ=255)
                                                                                                                    cgcgGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTgg                                                                                                                               <  1:211596/251‑1 (MQ=255)
                                                                                                                            acgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAAtt                                                                                                                       >  1:484833/1‑251 (MQ=255)
                                                                                                                                 ctACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCAt                                                                                                                    <  1:412115/249‑1 (MQ=255)
                                                                                                                                   gcgcgCTCTTCATCCACGGTTTCATCTGCGTATTCGTCAATTGCGTCTACGAGTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGGCGGACGTCGAGCGTCAGTGTCATCAAATTCAATAGCAACAAAATTGATCGTCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGGCGAGAACGACCAGTACGTCTATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGGTTGAATTCATc                                                                                                                   <  1:323121/247‑1 (MQ=255)
                                                                                                                                   aCGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAAGTTCTTCTTCGATTTCCGGTTCCGCAGTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGGGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATc                                                                                                                   <  1:425659/248‑1 (MQ=255)
                                                                                                                                    cgcgCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCAtt                                                                                                               <  2:559742/251‑1 (MQ=255)
                                                                                                                                    cgcgCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTGCCTCTTCGATGTCCAGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCCGCGGTTCGCCACGGTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGGACGACGATCGCTGGGATCGGGTTTCAGCTCGTCGTCTGGTGTGAAGTCATCCAtt                                                                                                               <  2:277057/251‑1 (MQ=255)
                                                                                                                                     gcgcTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCAttt                                                                                                              >  1:292130/1‑251 (MQ=255)
                                                                                                                                      cgcTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCAttt                                                                                                              >  2:132770/1‑250 (MQ=255)
                                                                                                                                       gcTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTaa                                                                                                            <  2:314071/251‑1 (MQ=255)
                                                                                                                                          ctTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAac                                                                                                           <  2:45924/249‑1 (MQ=255)
                                                                                                                                            tCATCCACTGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGTTTTCCGGTTCCTCATTGCGCTCTTTTTGCTCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAAGATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACTCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAacac                                                                                                         <  1:105579/249‑1 (MQ=255)
                                                                                                                                            tCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCAGCAAAAGTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAaca                                                                                                          <  1:590964/248‑1 (MQ=255)
                                                                                                                                             cATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACAccc                                                                                                       <  2:459576/250‑1 (MQ=255)
                                                                                                                                                   cGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTACCACCCCACTac                                                                                                 >  2:194611/1‑249 (MQ=255)
                                                                                                                                                     gTTTCATCTTCGGATTCGTCTACTTCTACTTCGATTTCCAGTTCCTCATTGCGCTCTTTTTGCGCACTAGTCGGACGGCGATCGTCAGTCTCATCGAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGCCCAGTACGACGATCGCTAGGTGCGGGTTTCAGCTCGTCTTCTGGGTTGAATTCATCCAGTTAACACCCCTCTAAAAg                                                                                              <  2:244006/251‑1 (MQ=255)
                                                                                                                                                             ttCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGc                                                                                      >  1:327391/1‑251 (MQ=255)
                                                                                                                                                              tCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAAt                                                                                        <  1:219030/248‑1 (MQ=255)
                                                                                                                                                                ggATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTa                                                                                   <  1:435824/251‑1 (MQ=255)
                                                                                                                                                                 gATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTAc                                                                                  <  1:452764/251‑1 (MQ=255)
                                                                                                                                                                    tCGTCAAGTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCTCCCGAGTCGGACGGCGCTCGTCAGTGTCATCAAATTCAATATCACCAAAATTGATCTGCGTTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTTCGACGATCGCTGGGATCGGGTTTCATCTCGTCTTCTTGTTTGAATTCATCCATTTATCACCCCTCTAAAAGGTTAATGCTTACCa                                                                                <  2:23166/250‑1 (MQ=255)
                                                                                                                                                                      gTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAGTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACTACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTTGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGt                                                                             <  1:52178/251‑1 (MQ=255)
                                                                                                                                                                       tCAATTTCTTCATCAATTTCCGGTTCCTCATTGCGCTCTTTTTGCGAACGAGTCGGACGTCGATCATCAGTTTCATCACATTCAATATCATAAAAATTGTTCTGCAGTTAGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGCCGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCGTTTCACACCCCACTACAAGGTTAATGCTTActtcgtt                                                                            >  2:559761/1‑245 (MQ=255)
                                                                                                                                                                         aaTTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGc                                                                          <  1:92403/251‑1 (MQ=255)
                                                                                                                                                                                tcttcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAGCGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAAtt                                                                      <  2:527225/248‑1 (MQ=255)
                                                                                                                                                                                tcttcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAAtt                                                                      <  1:132770/248‑1 (MQ=255)
                                                                                                                                                                                tcttcGATTTCCGGCTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGGGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTGCCACGTTGCAATTaa                                                                    <  2:324309/250‑1 (MQ=255)
                                                                                                                                                                                 cttcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCCAATTCAATATCATCAAAATTGCTCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAAcc                                                                  >  1:613086/1‑251 (MQ=255)
                                                                                                                                                                                    cGATGTCCGGTTCCTCATTGCGCTCATCTGGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCTGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACTACTATCGCTGGGATCTGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTACCCTGa                                                               <  2:351993/251‑1 (MQ=255)
                                                                                                                                                                                       tttCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGc                                                            >  1:290006/1‑251 (MQ=255)
                                                                                                                                                                                            ggTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGAACGCTCAGAACGTTCAGAAGATTGACGAGAACGACAAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTACTGCTTACCACGTTGCAATTAACCTGAAGCTAAt                                                        >  1:397219/1‑250 (MQ=255)
                                                                                                                                                                                                 cTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGAccc                                                  >  2:542454/1‑250 (MQ=255)
                                                                                                                                                                                                  tCATTGCGCTCTTTAAGCGCACGCGTCGGACGGCGATCGTCAGTGTCATCAATTTCACTATGATAAAAATTGATCTGGGGTTCGCAACGTTCAGAACGATCAGAACGTTCAGAAGATTGCCGAGACAGAACAGGACGACGATCGCTGGGCTAGGGATTCAGCTAGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAATGGTTAATGCTTAACACGTTGAAATTAACCTGACGCTAATAGAACGc                                                 >  2:422543/1‑251 (MQ=255)
                                                                                                                                                                                                         gcTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATaa                                          <  1:41422/251‑1 (MQ=255)
                                                                                                                                                                                                                    gcACGAGTCGGACGGCGATCGTCAGTCTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGAcc                               <  2:195134/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTAAGTACGCTCAGAACGTTCAGAAGATTGACGTGAACGACCAGTACGGCGATCGCTGGGATCGGGTTTCAGCGCGTCTTATGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAAACTGAAGCTAATAGACCGCTTGATAAGCGGCATGACCtt                             >  1:579602/2‑251 (MQ=255)
                                                                                                                                                                                                                      acgcgTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAATGATCTGCGGTTCGCCACGTTCAGTCCACTCAGAACGTTCAGACGATTAACGAGACCGACCAGTACGGCGATAGATGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGATTTCATCCATTTAACACCCCACTAAAAGGTTGATGCTTACAACGTTGCAATTACCCTGAAGCTAATAGACTGCTTGATAAGCGGCCTGACCtt                             >  1:579598/5‑251 (MQ=255)
                                                                                                                                                                                                                      aCGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCtt                             <  1:299031/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAATGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGTTTAATGCTTACCACGTTGCAATTAACCTGAATCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTa                  <  2:216863/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGc                                                                                      <  1:501414/180‑1 (MQ=255)
                                                                                                                                                                                                                                    gATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGc                                                                                      >  2:501414/1‑180 (MQ=255)
                                                                                                                                                                                                                                    aaTCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGAACAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGCTACGc               >  2:36315/2‑251 (MQ=255)
                                                                                                                                                                                                                                           aGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGttgattga        <  1:482519/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 atcaAATTCACTATCATCAAACTTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCACTTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGCCAATCg  >  2:291715/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
AGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAATCG  >  NC_000913/3516925‑3517416

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